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Abstract. We show the existence of global-in-time weak solutions to a general class of coupled 
FENE-type bead-spring chain models that arise from the kinetic theory of dilute solutions 
of polymeric liquids with noninteracting polymer chains. The class of models involves the 
unsteady incompressible Navier— Stokes equations in a bounded domain in R"^, d = 2 or 3, 
for the velocity and the pressure of the fluid, with an elastic extra-stress tensor appearing on 
the right-hand side in the momentum equation. The extra-stress tensor stems from the random 
movement of the polymer chains and is defined by the Kramers expression through the associated 
probability density function that satisfies a Fokker-Planck-type parabolic equation, a crucial 
feature of which is the presence of a centre-of-mass diffusion term. We require no structural 
assumptions on the drag term in the Fokker-Planck equation; in particular, the drag term need 
not be corotational. With a square-integrable and divergence-free initial velocity datum uq 
for the Navier— Stokes equation and a nonnegative initial probability density function ipo for 
the Fokker-Planck equation, which has finite relative entropy with respect to the Maxwellian 
M, we prove, via a limiting procedure on certain regularization parameters, the existence of 
a global-in-time weak solution t i— > {u{t), •ip{t)) to the coupled Navier-Stokes-Fokker— Planck 
system, satisfying the initial condition («(0), ■0(0)) = («o>V'o)i such that t i— > u{t) belongs to 
the classical Leray space and t i— f ^{t) has bounded relative entropy with respect to M and 
t I— > jp{t)/M has integrable Fisher information (w.r.t. the measure d/i := M{q)dq dx) over 
any time interval [0, T], T > 0. If the density of body forces / on the right-hand side of the 
Navicr-Stokes momentum equation vanishes, then a weak solution constructed as above is such 
that t I— (u{t), V'(i)) decays exponentially in time to (0 , M) in the X norm, at a rate that 
is independent of (no, V'o) and of the centre-of-mass diffusion coefficient. 

Keywords: Kinetic polymer models, FENE chain, Navier-Stokes-Fokker— Planck system. 



1. Introduction 

This paper establishes the existence of global-in-time weak solutions to a large class of bead- 
spring chain models with finitely extensible nonlinear elastic (FENE) type spring potentials, — 
a system of nonlinear partial differential equations that arises from the kinetic theory of dilute 
polymer solutions. The solvent is an incompressible, viscous, isothermal Newtonian fluid confined 
to a bounded open Lipschitz domain f2 C M'*, d = 2 or 3, with boundary dfl. For the sake of 
simplicity of presentation, we shall suppose that f2 has a 'solid boundary' dil; the velocity field 
u will then satisfy the no-slip boundary condition m = on dfl. The polymer chains, which 
are suspended in the solvent, are assumed not to interact with each other. The conservation of 
momentum and mass equations for the solvent then have the form of the incompressible Navier- 
Stokes equations in which the elastic extra-stress tensor r (i.e. the polymeric part of the Cauchy 
stress tensor) appears as a source term: 
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Given T e R>o, find u : {x,t) eVLx [0, T] ^ u{x, t) £ R'^ and p : {x,t) enx (0, T] ^ p{x, t) e 
R such that 



du 

(1.1a) -^ + {u-Yx)u-iyA.^u + V.^.p^ f + Vx-T inl]x(0,T], 

(1.1b) W^-u^O inr2x(0,T], 

(1.1c) u = ondnx{0,T], 

(l.ld) u{x,0)^uoix) \/xen. 



It is assumed that each of the equations above has been written in its nondimensional form; u 
denotes a nondimensional velocity, defined as the velocity field scaled by the characteristic flow 
speed Uq; 1/ G K>o is the reciprocal of the Reynolds number, i.e. the ratio of the kinematic 
viscosity coefhcient of the solvent and LqUq, where Lq is a characteristic length-scale of the flow; 
p is the nondimensional pressure and / is the nondimensional density of body forces. 

In a bead-spring chain model, consisting oi K + 1 beads coupled with K elastic springs to 
represent a polymer chain, the extra-stress tensor r is defined by the Kramers expression as 
a weighted average of ipj fhs probability density function of the (random) conformation vector 
q := {qj, . . . ,qJc)'^ & R^'^ of the chain (cf. (II. 7|) below), with q.i representing the d-component 
conformation/orientation vector of the ith spring. The Kolmogorov equation satisfied by ■0 is a 
second-order parabolic equation, the Fokker-Planck equation, whose transport coefficients depend 
on the velocity field u. The domain D of admissible conformation vectors D C R^'* is a K-fold 
Cartesian product Di x ■ ■ ■ x Dk of balanced convex open sets Di C R'', i — 1, . . . , if; the term 
balanced means that qi e Di if, and only if, —qi G Di. Hence, in particular, G Di, i = 1, . . . ,K. 
Typically Di is the whole of R"* or a bounded open d-dimensional ball centred at the origin e R'' 
for each i — 1, . . . ,K . When K — 1, the model is referred to as the dumbbell model. 

Let Oi C [0, oo) denote the image of Di under the mapping qi E Di i-^ ^kiP; and consider the 
spring potential [/^ G (7^(0^; R>o), i — I, . . . ,K. Clearly, G Oi. We shall suppose that Ui{0) = 
and that Ui is monotonic increasing and unbounded on Oi for each i = 1, . . . ,K. The elastic 
spring-force : C R'' — >• R'' of the ith spring in the chain is defined by 

(1.2) F,(g,) = C/;(i|g,|2)g„ t^l,...,K. 

Example 1.1. In the Hookean dumbbell model K — 1, and the spring force is defined by F{q) = q, 
with q E D = R'', corresponding to U{s) ~ s, s E O ^ [0, oo). This model is physically unrealistic 
as it admits an arbitrarily large extension, o 

We shall therefore assume in what follows that Z) is a Cartesian product of K bounded open 
balls Di C R'', centred at the origin E R'^ , i ^ 1, . . . , K, with K > 1. 

We shall further suppose that ioi i — 1, . . . ,K there exist constants Cy > 0, j = 1, 2, 3, 4, and 
7j > 1 such that the (normalized) Maxwellian Mi, defined by 

M.(,0 = ^e-^'^^l?'l^ Z.:./ ^'^m-^'^dq., 

A JDi 

and the associated spring potential Ui satisfy 

(1.3a) c,i [dist(g„ 5A)]^' < ^h{q^) < ca [dist(g„ dD,)y^ Vg, E D„ 

(1.3b) Q3 < [dist(q„ dDi)] U[{\\qJ'') < c,4 Vg, E D,. 

The Maxwellian in the model is then defined by 

K K 

(1.4) M(g) :=[]M,(9,) ^q := {qj , . . . ,ql)^ E D := X D.. 

~ 1=1 ~ ~ ~ ~ 

Observe that, for i ~ 1, . . . , K, 

(1.5) M(g)y,jM(g)]-i = -[M(g)]-iy,.Af(g) = YM^kJ') = ^1(^19.1')?.. 
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Since [[/^(ilgip)]^ = (— logMi{qi) + Const. )-^, it follows from (|1.3a[ b) that (if -fi > 1, as has been 
assumed here,) 

(1.6) / [l + [C/,(i|g,|2)]2 + [[/;( i|5,|2)]2lM,(g,)dg,<oo, t^l,...,K. 



Example 1.2. In the FENE (finitely extensible nonlinear elastic) dumbbell model K — 1 and 
the spring force is given by F{q) = (1 — \q\'^/b)^^ q, q £ D = B{0, b^), corresponding to U{s) = 

— |log (l — ^), s e O = [0, |). Here B{0,b^) is a bounded open bah in M'* centred at the origin 
e R'' and of fixed radius bi , with 6 > 0. Direct calculations show that the Maxwellian M and 
the elastic potential U of the FENE model satisfy the conditions (|1.3a[ b) with K = 1 and 7 := | 
provided that b> 2. Thus, (|1.6p also holds for K ^ 1 and b> 2. 

It is interesting to note that in the (equivalent) stochastic version of the FENE model a solution 
to the system of stochastic differential equations associated with the Fokker-Planck equation exists 
and has trajectorial uniqueness if, and only if, 7 = | > 1; (cf. Jourdain, Lelievre & Le Bris [5S] 
for details). Thus, in the general class of FENE-type bead-spring chain models considered here, 
the assumption 7^ > 1, j = 1, . . . ,K, is the weakest reasonable requirement on the decay-rate of 
M, in (fOal) as dist(gj, dD,) ^0. o 

The governing equations of the general FENE-type bead-spring chain model with centre-of-mass 
diffusion are (|l.lal -d). where the extra-stress tensor r is defined by the Kramers expression: 

(1-7) Zix, t) = k{^J^ V(«, g, t) qj U[ (\\qj^) dg - p{x, t) / j , 

with the density of polymer chains located at x at time t given by 
(1.8) p{x,t) = I ^ix,q,t)dq. 



J D 

The probability density function ip is a, solution of the Fokker-Planck equation 
^ + (u • V,)V' + I] V,. • (a{u) q, ^ 

(1.9) ^eA.V' + ^EE^'^ Y?. -(^^Y?. (^)) inf^xi^x (0,T], 

4=1 j = l ~ \ ~ \ ' / / 

with a{y) = V^jj;, where xU){x,t) e W'''^'^ and {'^xV^ij — The dimensionless constant 

k > featuring in (II. 7p is a constant multiple of the product of the Boltzmann constant kg and 
the absolute temperature T. In (|1.9p . e > is the centre-of-mass diffusion coefficient defined as 
£ ■= i^a / Lq)'^ / {4:{K + 1)A) with £0 := ^fc^T/H signifying the characteristic microscopic length- 
scale and A :— (C/4H)(J7o/io), where C > is a friction coefficient and H > is a spring-constant. 
The dimensionless parameter A G IR>o, called the Weissenberg number (and usually denoted by 
Wi), characterizes the elastic relaxation property of the fiuid, and A = {Aij)fj^i is the symmetric 
positive definite Rouse matrix, or connectivity matrix; for example, A — tridiag [— 1, 2, — 1] in 
the case of a linear chain; see, Nitta 40 . 

Definition 1.1. The collection of equations and structural hypotheses ll.la\ -d)- H1.9^ will be re- 
ferred to throughout the paper as model (Pg), or as the general FENE-type bead-spring chain 
model with centre-of-mass diffusion. 

A noteworthy feature of equation (|1.9p in the model (Pg) compared to classical Fokker-Planck 
equations for bead-spring models in the literature is the presence of the j-dissipative centre-of- 
mass diffusion term e A^t/j on the right-hand side of the Fokker-Planck equation (|1.9p . We refer to 
Barrett & Siili ^ for the derivation of (|1.9|) in the case oi K — 1; see also the article by Schieber 
[45j concerning generalized dumbbell models with centre-of-mass diffusion, and the recent paper 
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of Degond & Liu [T7] for a careful justification of the presence of the centre-of-mass diffusion 
term through asymptotic analysis. In standard derivations of bead-spring models the centre-of- 
mass diffusion term is routinely omitted on the grounds that it is several orders of magnitude 
smaller than the other terms in the equation. Indeed, when the characteristic macroscopic length- 
scale Lq ~ 1, (for example, Lq = diam(il)), Bhave, Armstrong & Brown [M] estimate the ratio 
£q/Lq to be in the range of about 10~^ to 10"''. However, the omission of the term e Axi^ from 
(|1.9p in the case of a heterogeneous solvent velocity u{x, t) is a mathematically counterproductive 
model reduction. When e A^il) is absent, (|1.9p becomes a degenerate parabolic equation exhibiting 
hyperbolic behaviour with respect to {x,t). Since the study of weak solutions to the coupled 
problem requires one to work with velocity fields u that have very limited Sobolev regularity 
(typically u S L°°(0, T; L'^{Q))C]L'^{Q, T; Hl{fl))), one is then forced into the technically unpleasant 
framework of hyperbolically degenerate parabolic equations with rough transport coefhcients (cf. 
Ambrosio [5] and DiPerna &: Lions [T^). The resulting difhculties are further exacerbated by 
the fact that, when D is bounded, a typical spring force F{q) for a finitely extensible model 
(such as FENE) explodes as q approaches dD; see Example 11.21 above. For these reasons, here 
we shall retain the centre-of-mass diffusion term in (jl.9p . In order to emphasize that the positive 
centre-of-mass diffusion coefficient e is not a mathematical artifact but the outcome of the physical 
derivation of the model, in Section [5] and thereafter the variables u and ip have been labelled with 
the subscript e. 

We continue with a brief literature survey. Unless otherwise stated, the centre-of-mass diffusion 
term is absent from the model considered in the cited reference (i.e. e is set to 0); also, in all 
references cited K — 1, i.e. a simple dumbbell model is considered rather than a bead-spring chain 
model. 

An early contribution to the existence and uniqueness of local-in-time solutions to a family of 
dumbbell type polymeric flow models is due to Renardy |33]. While the class of potentials F{q) 
considered by Renardy [44 (cf. hypotheses (F) and (F') on pp. 314-315) does include the case of 
Hookean dumbbells, it excludes the practically relevant case of the FENE dumbbell model (see 
Example 1 1 . 21 above) . More recently, E, Li & Zhang 23 and Li, Zhang & Zhang 321 have revisited 
the question of local existence of solutions for dumbbell models. A further development in this 
direction is the work of Zhang & Zhang (55], where the local existence of regular solutions to 
FENE-type dumbbell models has been shown. All of these papers require high regularity of the 
initial data. Constantin [15] considered the Navier-Stokes equations coupled to nonlinear Fokker- 
Planck equations describing the evolution of the probability distribution of the particles interacting 
with the fluid. Otto & Tzavaras [42^ investigated the Doi model (which is similar to a Hookean 
model (cf. Example 11.11 above ) . except that D = S^) for suspensions of rod-like molecules in the 
dilute regime. Jourdain, Lelievre & Le Bris ^5B^ studied the existence of solutions to the FENE 
dumbbell model in the case of a simple Couette flow. By using tools from the theory of stochastic 
differential equations, they showed the existence of a unique local-in-time solution to the FENE 
dumbbell model for d = 2 when the velocity fleld u is unidirectional and of the particular form 

In the case of Hookean dumbbells {K = 1), and assuming e = 0, the coupled microscopic- 
macroscopic model described above yields, formally, taking the second moment of g i— > ip{q,x,t), 
the fully macroscopic, Oldroyd-B model of viscoelastic flow. Lions & Masmoudi 35! have shown 
the existence of global-in-time weak solutions to the Oldroyd-B model in a simplified corotational 
setting (i.e. with cr(u) = Yx U replaced by ^{Yx y~(yx u)^)) by exploiting the propagation in time 
of the compactness of the solution (i.e. the property that if one takes a sequence of weak solutions 
that converges weakly and such that the corresponding sequence of initial data converges strongly, 
then the weak limit is also a solution) and the DiPerna-Lions theory of renormalized solutions 
to linear hyperbolic equations with nonsmooth transport coefficients. It is not known if an identical 
global existence result for the Oldroyd-B model also holds in the absence of the crucial assumption 
that the drag term is corotational. With e > 0, the coupled microscopic-macroscopic model above 
yields, taking the appropriate moments in the case of Hookean dumbbells, a dissipative version 
of the Oldroyd-B model. In this sense, the Hookean dumbbell model has a macroscopic closure: 
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it is the Oldroyd-B model when e = 0, and a dissipative version of Oldroyd-B when e > (cf. 
Barrett & SiiH [9J). Barrett & Boyaval [7] have proved a global existence result for this dissipative 
Oldroyd-B model in two space dimensions. In contrast, the FENE model is not known to have 
an exact closure at the macroscopic level, though Du, Yu & Liu [20] and Yu, Du & Liu [51] have 
recently considered the analysis of approximate closures of the FENE dumbbell model. Lions & 
Masmoudi [36] proved the global existence of weak solutions for the corotational FENE dumbbell 
model, once again corresponding to the case of e = and K — 1, and the Doi model, also called 
the rod model; see also the work of Masmoudi [37] ■ Recently, Masmoudi [35] has extended this 
analysis to the noncorotational case. 

Previously, El-Kareh & Leal [21] had proposed a steady macroscopic model, with added dissi- 
pation in the equation satisfied by the conformation tensor, defined as 



in order to account for Brownian motion across streamlines; the model can be thought of as an 
approximate macroscopic closure of a FENE-type micro-macro model with centre-of-mass diffu- 
sion. 

Barrett, Schwab & Siili [8] showed the existence of global weak solutions to the coupled 
microscopic- macroscopic model (|l.lal -d). (II. 9p with e — 0, K ~ 1, a.n x- mollified velocity gradient 
in the Fokker-Planck equation and an x-mollified probability density function -0 in the Kramers 
expression, admitting a large class of potentials U (including the Hookean dumbbell model and 
general FENE-type dumbbell models); in addition to these mollifications, u in the x-convective 
term (w • in the Fokker-Planck equation was also mollified. Unlike Lions & Masmoudi [35] . 

the arguments in Barrett, Schwab & Siili [8J did not require that the drag term • (ct(w) q ip) in 
the Fokker-Planck equation was corotational in the FENE case. 

In Barrett & Siili ^ , we derived the coupled Navier-Stokes-Fokker-Planck model with centre- 
of-mass diffusion stated above, in the case oi K = \. We established the existence of global- in- 
time weak solutions to a mollification of the model for a general class of spring-force-potentials 
including in particular the FENE potential. We justified also, through a rigorous limiting process, 
certain classical reductions of this model appearing in the literature that exclude the centre-of- 
mass diffusion term from the Fokker-Planck equation on the grounds that the diffusion coefficient 
is small relative to other coefficients featuring in the equation. In the case of a corotational drag 
term we performed a rigorous passage to the limit as the mollificrs in the Kramers expression and 
the drag term converge to identity operators. 

In Barrett & Siili [TU] we showed the existence of global-in-time weak solutions to the general 
class of noncorotational FENE type dumbbell models (including the standard FENE dumbbell 
model) with centre-of-mass diffusion, in the case of K — \, with microsropic cut-off (cf. (|l.lll) 
and (|1.12l) below) in the drag term 



In this paper we prove the existence of global-in-time weak solutions to the model without cut-off 
or mollification, in the general case oi K > 1. Since the argument is long and technical, we give a 
brief overview of the main steps of the proof here. 

Step 1. Following the approach in Barrett & Siili [TO] and motivated by recent papers of 
Jourdain, Lelievre, Le Bris & Otto [29] and Lin, Liu & Zhang [33] (see also Arnold, Markowich, 
Toscani & Unterreiter [6], and Desvillettes & Villani [18]) concerning the convergence of the 
probability density function tp to its equilibrium value ipoc>ix,q) := M{q) (corresponding to the 
equilibrium value Uca{x) '.= of the velocity field) in the absence of body forces /, we observe 
that if ^p/M is bounded above then, for L e M>o sufficiently large, the drag term (|1.10p is equal 
to 




(1.10) 




(1.11) 
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where G C(R) is a cut-off function defined as 

(1.12) ^^(s) := niin(s,L). 

More generally, in the case of X > 1 , in analogy with (|l.lip , the drag term with cut-off is defined 
by Yqi ■ (sin) 1i ^ (^) ) ■ then follows that, for L ^ 1, any solution ijj of (|1.9p . such 

that ip/M is bounded above, also satisfies 

K 



K K 



(1.13) =£^-V^+2T^^^^ 

1=1 j=l 

We impose the following boundary and initial conditions 



''""-ill 



fix Dx (0,T]. 



M 



K 



a{u)qil3 



(1.14a) 
(1.14b) 



on f2 X dDi X 



K 

X 




X (0,r], for i = 1, 



(1.14c) 



on f2 X Z3, 



\3 = '^,i¥=i 

eWa^i/j-n^O on dn X D X {0,T], 

^(•,-,0) = A/(-)/3^(^o(-,-)A^(-))>0 
where qi is normal to dDi, as Di is a bounded ball centred at the origin, and n is normal to dQ; Tpo 
is nonnegative, defined on 51 x D, with "00(2;, <z) df/ = 1 for a.e. x G J7, and assumed to have finite 
relative entropy with respect to the Maxwelfian M; i.e. /q^u ^o(S' i) log(V'o(S' ^9 dj < 

oo. Clearly, if there exists L > such that < i^o < LM, then M f^^^^po/M) = ipo. Henceforth 
i > 1 is assumed. 



Definition 1.2. The coupled problem \l.la\ -d). |J.7| ), dH^j; il.l3\) . {1.14a\ -c) will he referred 
to as model (P^^i), or as the general FENE-type bead-spring chain model with centre-of-niass 
diffusion and microscopic cut-off, with cut- ojf parameter L > 1. 



In order to highlight the dependence on e and L, in subsequent sections the solution to (jl.l3p . 
(I1.14af c) will be labelled V'e.L- Due to the coupling of (|1.13p to (|l.lap through (|1.7p . the velocity 
and the pressure will also depend on e and L and we shall therefore denote them in subsequent 
sections by Ue,L and Pe.L- 

The cut-off 13^ has a convenient property: the couple {uoojipoo), defined by ^00(5) ■— Q and 
i^<yo{x,q) := M{q), is still an equilibrium solution of pITal -dl with / = 0, dTTTj) . (fTS)) . (fTTTO)) . 
(I1.14a| -c) for all L > 0. Thus, unlike the truncation of the (unbounded) potential proposed in 
El-Kareh & Leal [24] , the introduction of the cut-off function into the Fokker-Planck equation 
(II. 9p does not alter the equilibrium solution (uqo, "iAoo) of the original Navier-Stokes-Fokker-Planck 
system. In addition, the boundary conditions for on 80, x D x {0, T] and il x dD x (0, T] ensure 
that Jjy tpis^ 9; = Jjj ^{x, q, 0) dq for a.e. x £ D. and a.e. t e M>o- 

Step 2. Ideally, one would like to pass to the limit L oo in problem {Pe,L) to deduce the 
existence of solutions to (Pe). Unfortunately, such a direct attack at the problem is (except in 
the special case oi d — 2, or in the absence of convection terms from the model,) fraught with 
technical difficulties. Instead, we shall first (semi)discretize problem (Pc.l) by an implicit Euler 
scheme with respect to t, with step size At; this then results in a time-discrete version (Pf^^) 
of (Pe,L)- By using Schauder's fixed point theorem, we will show in Section [3] the existence of 
solutions to (P^^). In the course of the proof, for technical reasons, a further cut-off, now from 
below, is required, with a cut-off parameter S e (0, 1), which we shall let pass to to complete 
the proof of existence of solutions to (P^£) in the limit of 5 — >■ 0+ (cf. Section [S]). Ultimately, 
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of course, our aim is to show existence of weak solutions to the general FENE-type bead-spring 
chain model with centre-of-mass diffusion, (Pe), and that demands passing to the limits At — > 0+ 
and L — > oo; this then brings us to the next step in our argument. 

Step 3. We shall link the time step At to the cut-off parameter L > 1 by demanding that 
At = o{L~^), as 1/ — >■ oo, so that the only parameter in the problem (P^^) is the cut-off parameter 
(the centre-of-mass diffusion parameter e being fixed). By using special energy estimates, based 
on testing the Fokker-Planck equation in (Pf^^) with the derivative of the relative entropy with 
respect to the Maxwellian of the general FENE-type bead-spring chain model, we show that 
can be bounded, independent of L. Specifically is bounded in the norm of the classical Leray 
space, independent of L; also, the L°° norm in time of the relative entropy of ip^L norm 

in time of the Fisher information of ij^'j^ are bounded, independent of L. We then 

use these L-independent bounds on the relative entropy and the Fisher information to derive L- 
independent bounds on the time-derivatives of and V'i^i in very weak, negative-order Sobolev 
norms. 

Step 4- The collection of L-independent bounds from Step 3 then enables us to extract a 
weakly convergent subsequence of solutions to problem (P^^) as L — >■ oo. We then apply a general 
compactness result in seminormed sets due to Dubinskii |21j (see also |13j). which furnishes strong 
convergence of a subsequence of solutions iuf^L^ii^fi^^) to (P^i) with At = o{L^^) as i — > oo, 
in L2(o,T;L2(rj)) x LP{0,T;L^{n x D)) for any p > I. A crucial observation is that the set 
of functions with finite Fisher information is not a linear space; therefore, typical Aubin-Lions- 
Simon type compactness results (see, for example, Simon |47| ) do not work in our context; however, 
Dubinskii's compactness theorem, which applies to seminormed sets in the sense of Dubinskii, does, 
enabling us to pass to the limit with the microscopic cut-off parameter L in the model (P^£), 
with At = o{L^^), as L — > oo, to finally deduce the existence of a weak solution to the general 
FENE-type bead-spring chain model with centre-of-mass diffusion, (P^). 

The paper is structured as follows. We begin, in Section [21 by stating (Pe^l), the coupled 
Navier-Stokes-Fokker-Planck system with centre-of-mass diffusion and microscopic cut-off for a 
general class of FENE-type spring potentials. In Section |3] we establish the existence of solutions 
to the time-discrete problem (P^^). In Section |4] we derive a set of L-independent bounds on 
in the classical Leray space, together with L-independent bounds on the relative entropy of tp^j^ 
and Fisher information of ip^l^- We then use these L-independent bounds on spatial norms to 

obtain L-independent bounds on very weak norms of time-derivatives of and ■^^l,. Section [S] 
is concerned with the application of Dubinskii's theorem to our problem; and the extraction of a 
strongly convergent subsequence, which we shall then use in Section[5]to pass to the limit with the 
cut-off parameter L in problem (P^|^), with At ~ o(L~^), as L — > oo, to deduce the existence of a 

weak solution {u^, := M tp^) to problem (Pg), the general FENE-type bead-spring chain model 
with centre-of-mass diffusion. Finally, in Section [71 we show using a logarithmic Sobolev inequality 
and the Csiszar-KuUback inequality that, when / = 0, global weak solutions t i— ?> {u^{t),ip^{t)) 
thus constructed decay exponentially in time to (0,M), at a rate that is independent of the 
initial data for the Navier-Stokes and Fokker-Planck equations and of the centre-of-mass diffusion 
coefficient e. We shall operate within Maxwellian-weighted Sobolev spaces, which provide the 
natural functional-analytic framework for the problem. Our proofs require special density and 
embedding results in these spaces, which are proved in Appendix C and Appendix D, respectively. 

For an analogous set of existence and equilibration results for weak solutions of Hookean-type 
bead-spring chain models for dilute polymers, we refer to Part II of the present paper [11) . 



2. The polymer model (Pe,L) 

Let C M'^ be a bounded open set with a Lipschitz-continuous boundary d^, and suppose that 
the set D := Di X • • • X Dk of admissible conformation vectors q := {qj , . . . , qj^)'^ in (|1.9p is such 
that Di, i = 1, . . . , if , is an open ball in R'', d = 2 or 3, centred at the origin with boundary dDi 
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and radius v^, &i > 2; let 

(2.1) ai3:=Uli[aAx(xf=,,,^,7?,)]. 

Collecting (ll.lal ~d). (|1.7p . and (|1.9p . we then consider the following initial-boundary- value prob- 
lem, dependent on the parameter L > I. As has been already emphasized in the Introduction, the 
centre-of-mass diffusion coefficient £ > is a physical parameter and is regarded as being fixed, 
although we systematically highlight its presence in the model through our subscript notation. 



(P,.l) Find ye,L : {x,t) e n x [0,T] ^ Ue,L{x,t) e 
p^^L{x,t) e M such that 



dt 



(2.2a) 
(2.2b) 
(2.2c) 
(2.2d) 



Yx ■ Us,L 
He.L 
y:e.L{x,0) 



where 4's,l '■ {x,q,t) e fl x D x [0,T] i-^ ipe,L{x,q,t) G 
l('0e.L)(a;, <) e M.'^^'^ is the symmetric extra-stress tensor defined as 



^ and pe,L ■■ {x,t) e n X (0,r] ^ 

-- f + Yx- li^ex) 

in 17 X (o,r], 

-- inl7x(0,r], 
-- onar2x(0,T], 
= ^0(5) Va: e fl, 
and t(V'c,l) : {x,t) £ n x (0,T] ^ 



(2.3) 



Here k € ]R>o, I is the unit d x d tensor, 



(2.4a) 
(2.4b) 



C\i,P,,L)ix,t) 
P{i^e,L){x,t) 



'>Pe,L{x,q,t) Ul{^\q^f) qi qj dq, 



'^e,L{x,q,t) dq. 



and 



The Fokker-Planck equation with microscopic cut-off satisfied by V-'e,L is 

K 

+ (^ie,L • V^)^Ae,L + XlY9' 



5^, 



'e,L 



dt 



i=l 



(2.5) 



^ K K 



M 



1=1 j=i 



MV 



M 



inilx Dx (0,T]. 



Here, for a given L > 1, (3^ e C{M.) is defined by (|1.12p . ^(t;) = Yx v, and 
(2.6) A G R^^^ is symmetric positive definite with smallest eigenvalue ag G 
We impose the following boundary and initial conditions: 



M 



1 ^ 

rr X! ^^-J' Y'3 



2A 



M 



a{us^L)qtl3^ 



M 







(2.7a) 

(2.7b) 
(2.7c) 



K 



on X dDj X 



V'e,L(-,-,0)=M(.)/3^(Vo(-,-)/M(-))>0 



X DAx{Q,T], i = l,...,K, 
ondnx Dx (0,r], 



on f2 X _D, 



where ri is the unit outward normal to d^l. The boundary conditions for V-'e,L on d^l x D x (0,T] 
and ft X dD x (0, T] have been chosen so as to ensure that 



(2.8) 



■ipe,Lix,q,t) dq 



^/'e,L(x,q,0)dq y{x,t) G ft X {0,T]. 



EXISTENCE AND EQUILIBRATION OF GLOBAL WEAK SOLUTIONS FOR DILUTE POLYMERS 



9 



Henceforth, we shall write i^^L •= V's-l/^i i'o ■= i'o/M. Thus, for example, (|2.7cp in terms of 
this compact notation becomes: V'e,L(-, •, 0) = /3^('0o(-, ■)) on il x D. 

The notation | • | will be used to signify one of the following. When applied to a real number x, 
\x\ will denote the absolute value of the number x; when applied to a vector v, \v\ will stand for 
the Euclidean norm of the vector v; and, when applied to a square matrix A, \A\ will signify the 
Frobenius norm, [tr(yl'^74)] 5 , of the matrix A, where, for a square matrix B, tx{B) denotes the 
trace of B. 



3. Existence of a solution to the discrete-in-time problem 

Let 

(3.1) H ■.^{weL^{n):V,j:-w^O} and F {w G : • u; = 0}, 

where the divergence operator V^^- is to be understood in the sense of distributions on 51. Let V' 
be the dual of V. Let S : V' V he such that 5 u is the unique solution to the Helmholtz-Stokes 
problem 

(3.2) / S V • wdx + / Yx {§ v) : Yx wdx — {v,w)v Vu; € V, 

~ ~ ~ J^^~ ~ ~ ~ ~ ~ 

where (•, •)y denotes the duality pairing between V' and V. We note that 

(3.3) {v,Sv}y^\\Sv\\l,^,,^ yv&Y', 

and \\S ■ ||_ffi(si) is a norm on V' . More generally, let denote the closure of the set of all 
divergence-free C'^{ft) functions in the norm of HQ{il)r\H°'{fl), cr > 1, equipped with the Hilbert 
space norm, denoted by || • || y^, inherited from H"'{n), and let V'^ signify the dual space of Va, with 
duality pairing (•, ■)v^- As £7 is a bounded Lipschitz domain, we have that Vi = V (cf. Temam 
Hn], Ch. 1, Thm. 1.6). Similarly, (•, •)_f/i(n) 'v^iU denote the duality pairing between (ffj(fi))' and 
i?J(r2). The norm on {^^{il))' will be that induced from taking WVx ■ \\L^(n) to be the norm on 

Him- 

For later purposes, we recall the following well-known Gagliardo-Nirenberg inequality. Let 
r £ [2,oo) if d = 2, and r G [2,6] if c? = 3 and 6 = (| — 7)- Then, there is a constant 
C = Cin,r,d), such that, for all 77 e H^ifl): 

(3.4) ||r/|U.(o) <C^||r/||ii(%||ry||?,,(^). 

Let J" e C(M>o) be defined by J"(s) := s (logs - 1) -I- 1, s > 0. As lims^o+ -^(s) 1, the 
function T can be considered to be defined and continuous on [0, cx)), where it is a nonnegative, 
strictly convex function with J^(l) = 0. Wc assume the following: 

dn e C"'^ uoeH; V'o ^ > a.e. on fl x D with 



J^ii^o) X D) and / M{q)Mx,q) 

Jd ~ ~ ~ 



dg = 1 for a.e. x £ H.; 



(3.5) 7i>l, i = l,...,K in (HSib); and / e ^^(0, T; (i?J(17))'). 

Here, L^j{^ x D), for p E [l,oo), denotes the Maxwellian- weighted space over il x D with 
norm 



^hljinxD) \ M\(p\Pdqdx 



nxD 



Similarly, we introduce L\.j{D), the Maxwellian- weighted space over D. Letting 



dq dx 



(3.6) mHl,inxD) - (/ [i^i' + iy.^i' + iYg^i" 

we then set 

(3.7) X = Hlj{nxD) := e L1,,(0 x i^) : ||^||^_L^(f,,^) < 00} 
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It is shown in Appendix C that 

(3.8) C°°(f7 X D) is dense in X. 
We have from Sobolev embedding that 

(3.9) H\n;Ll{D)) ^ L^{n-Ll{D)), 

where s G [1, oo) if d = 2 or s G [1,6] if d = 3. Similarly to p.4p we have, with r and 6 as there, 
that there is a constant C, depending only on fi, r and d, such that 

(3.10) ML^in;Ll,iD)) < ^^ll^llil(0;Ll,(D)) Mm (n-.Ll.iD)) y^eH\n;LliD)). 
In addition, we note that the embeddings 

(3.11a) HIj{D) ^ Ll,{D), 

(3.11b) Hli{n X D) = L^{n;Hlj{D))r\H\n;Lli{D)) ^ LIj{VL x D) = {Q; LIj (D)) 

are compact if 7^ > 1, i = 1, . . . , if, in (ll.3al b'): see Appendix D. 

Let X' be the dual space of X with Lj^ji^ x D) being the pivot space. Then, similarly to p.2p . 
lei Q : X' ^ X be such that Qrj\s the unique solution of 



M 



nxD 



(3.12) 
where (M • 
(3.13) 



dq dx 

) is the duality pairing between X' and X. Then, as in 



2 

J? 



and 11^7 • II is a norm on X' . 

We recall the Aubin-Lions-Simon compactness theorem, see, e.g., Temam ^35" and Simon ^47] . 
Let Bq, B and Bi be Banach spaces, Bi, i = 0, 1, reflexive, with a compact embedding Bq ^ B 
and a continuous embedding B ^ Bi. Then, for > 1, i = 0, 1, the embedding 

(3.14) {r,€ L"o(0,r;So) : |f e (0, T; 61) } L"°(0,T;e) 

is compact. 

Throughout we will assume that p.5p hold, so that (II. 6p and p.llai b) hold. We note for future 
reference that (|2.4ap and (II. 6p yield that, for G L\i{fl x D), 



|C,(M(^)|2 dx 



(3.15) 



< 



< C 



MipU^q,qJdq 



A./({/;)2|g,|4 dq 



dx 



D 



M\(p\^dq dx 



flxD 



M\ip\^dq dx 



flxD 



1 = 1,. ..,K, 



where C is a positive constant. 

We establish a simple integration-by-parts formula. 

Lemma 3.1. Let ip e H\,j{D) and suppose that B E 
then, 



idxd 



is a square matrix such that tz{B) = 0; 



(3.16) 



/ M V(Bq,) •y,,(?d(z = / M 



K 



V ? ^q^qIU^{^\q.\^)■■Bdq. 



1=1 



Proof. By Theorem C.l in Appendix C, the set C°°{D) is dense in H\,j{D); hence, for any (p e 
H\,j[D) there exists a sequence {<^„}„>o C C°°[D) converging to (p in H\f{D). As M g C^{D) 
and vanishes on dD, the same is true of each of the functions Mipn, n > 1. On replacing If by 
(pn on both sides of (|3.16p , the resulting identity is easily verified by using the classical divergence 
theorem for smooth functions, noting (|1.5p . that M(pn vanishes on dD and that ix{B) = 0. Then, 
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p.l6|) itself follows by letting n oo, recalling the definition of the norm in H\j{D) and hypothesis 

(ira . □ 

We now formulate our discrete-in-time approximation of problem {Pe,L) for fixed parameters 
£ G (0, 1] and L > 1. For any T > and > 1, let Ai = T and i„ = n At, n = 0, . . . , iV. 
To prove existence of a solution under minimal smoothness requirements on the initial datum uq 
(recall P-Sp *). we introduce — y°(At) € V such that 



(3.17) 
and so 
(3.18) 



• w -)- At w" : Va; w 



dx — uq ■ vdx yv £ V] 
'n~ ~ ~ ~ ~ 



u"\^ + At\Y^u"r]dx< / \uo\^dx<C. 



In addition, we have that converges to uq weakly in H in the limit of At — !■ 0+. For p G [1, oo), 
let 



Zp := {if G (fi xD):if>0 a.e. on fl x D and 



(3.19) 



D 



M{q) (p{x, q)dq <1 for a.e. x G VI}. 



Analogously to defining iP for a given initial velocity field uq, we shall assign a certain 'smoothed' 
initial datum, ij^^ — 'ip'^{At), to the initial datum ipQ. The definition of ip'^ is delicate; it will be 
given in Section [6] All we need to know for now is that there exists a ip^, independent of the 
cut-off parameter L, such that: 



(3.20) 



G Zi; 



Tiip°)eL\,inx D); 



V° G Hij{n X D); JnxD 



M dqdx< / M J"(^o) dq dx. 



The proofs of these properties will be given in Lemma 16.21 in Section |6l It follows from (|3.20p and 
that /3^(^A°) G ^2; m fact, /3^(??°) e L°°{nx D)n HIj{VL x D). 
Our discrete-in-time approximation of (P£,l) is then defined as follows. 

(P^i) Let uIj^ := yO g F and -0^^ := /3^(<) G Z2. Then, for n = 1, . . . , A^, given 
(M"i\ e F X ^2, find [u^j^, 4^1^) eV X {XnZ2) such that 



,n-l 



At 



• It; da; + / V^; w"^ : V^: w dec 



(3.21a) 



^ f 



AI — — : — $ dq dx 

At ^ I ~ 



K 

JnxD „_i 



1 ^ ^ 

— ^ A, V,^ V"l - [ ^(<l) 9. ] /S"^ (Cl) 



i=i 



Vq.^dg dec 



(3.21b) 



M 



nxD 



where, for t G [t„-_i, t„), and n = 1, . . . , iV, 



(3.22) 



/^*^+(-,t) = /"(•) 



Va;!^ dq dx = \f(f G X; 



/(•,t) dt G (ill(l}))' C F'. 
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It follows from and ((X^ that 

(3.23) /'^*^+ ^ / strongly in L'^{0, T; as At 0+. 

Note that as the test function w in p.21al) is chosen to be divergence-free, the term containing 
the density p in the definition of r (cf. p.3p ) is eliminated from p.21al) . 

In order to prove the existence of a solution to (P^£), we require the following convex regular- 
ization e C^'^(M) of J" defined, for any S € (0, 1) and L > 1, by 



(3.24) 

Hence, 

(3.25a) 

(3.25b) 

We note that 
(3.26) 



^^ + s(log(5-l) + l fors<<5, 
T^-is) := { T{s) = s (logs - 1) + 1 for S < s < L, 
+ s (log L - 1) + 1 for L < s. 



His) > 



I + log (5 - 1 for s < S, 
logs for S < s < L, 

f + log L - 1 for L < s, 



- -1 



for s < S, 
ior S < s < L, 
for L < s. 



25 

TL-C[L) 



for s < 0, 
for s > 0; 



and that [^s]"{s) is bounded below by 1/L for all s e K. Finally, we set 



(3.27) 



Pns) m]")-\s) = max{/3^(s),<5}. 



and observe that /3|'(s) is bounded above by L and bounded below by 6 for all s e M. Note also 
that both and /3|' are Lipschitz continuous on M, with Lipschitz constants equal to 1. 

3.1. Existence of a solution to (P^^)- It is convenient to rewrite p.21ap as 



(3.28) 



where, for all Wi G i/J(il), i — 1,2, 



(3.29a) 



b{wi,W2): 



wi + At (u 



■W2 dx + At / Va; wi : Wx W2 dx. 



and, for all w <E ffKQ) and ip G L\[{Q x D), 
(3.29b) ib{ip){w) ■.^At{r\w)Hi^n)' 



K 



Atk^C^{M(f) : Vx 



w 



dx. 



We note that, for all u € y and all wi,W2 G H^{il), we have that 



(3.30) 



{v ■ Vx)wi 



W2 dx = — 



{V ■ Vx)w2 



wi dx 



and hence b{-,-) is a continuous nonsymmetric coercive bilinear functional on Hq(JI) x i/J(r2). 
In addition, thanks to (|3.15p . ib{'f){-) is a continuous hnear functional on Hl{Q) for any £ 
Llj{n X D). 
For r > d, let 



(3.31) 



y € L''{n) ■■ J^v ■ "^xwdx = Vw e wi'^(n)| 
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It is also convenient to rewrite p.21b|) as 

(3.32) a(^"^,^) = V^e X, 



where, for all <^i, <^2 G X, 
a(<?i,^2) 



K K 



(3.33a) +^EE^*^ 

i=i j=i 

and, for all v G H^fl), rj £ L°^{VL x D) and 0eX, 

K 

nxD ' .-1 ~ ~ ~ 



(3.33b) 



dq dx. 



It follows from (P3T|) and that, for r >d, 



(3.34) 



Mvip ■Yx$dqdx = Vu G y^ £ X. 



flxD 



Hence a(-, •) is a continuous coercive bilinear functional on X x X. In addition, we have that, for 
all V e H^i^), 77 e L°°{n X D) and X, 

Ka(^',^)(<?)l < IIC,I^llL|,(nxD) II^IIl5^(OxD) 



(3.35) 



At 



n-1 



Af|g| dq\ M\L--(ny.D)\\'^ xV\\L^(^n)\\'^ q$\\Ll^(nycD)- 



Therefore, by noting that t/;"^ G .^2 and recalling (|1.4p . it follows that (!-a{v,vj){-) is a continuous 

linear functional on X for all u G H'^{^) and 77 G L°°{VL x D). 

In order to prove existence of a solution to p.21al b). i.e. p.28p and p.32p . we consider a 
regularized system for a given 5 G (0, 1): 

Find {uIl,s^$Il,s) e F X ^ such that 

(3.36a) b{ulL,5, w) = ib{re,L,5){w) Vw G V, 

(3.36b) a{$lLj,^) = iaiul\,s,l3H$lL,s)m G X. 

The existence of a solution to p.36al b) will be proved by using a fixed-point argument. Given 

-0 G Lli{n X D), let {u*, ip*) e V X X he such that 



(3.37a) 
(3.37b) 



b{u* ,w) — £b{^)iw) Vw G 



The Lax-Milgram theorem yields the existence of a unique solution to p.37a[ b). and so the overall 
procedure p.37al b') is well defined. 

Lemma 3.2. Let G : i^f(^ x D) ^ X C L\,f{n x D) denote the nonlinear map that takes the 
function ip to ip* — G{^) via the procedure (|3.37a[ b). Then G has a fixed point. Hence there exists 
a solution {u^^l,S''$6,l,s) ^ Y x X to p.36a[ b). 

Proof. Clearly, a fixed point of G yields a solution of p.36al b'). In order to show that G has 
a fixed point, we apply Schauder's fixed-point theorem; that is, we need to show that: (i) G : 
L\f{^l X D) ^ X D) is continuous; (ii) G is compact; and (iii) there exists a C* G IR>o such 

that 

(3.38) mLl,inxD) < a 

for every ^ G Lj^{Q x D) and k G (0, 1] satisfying ^ = k G{'ip). 
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Let {V'^^''}p>o be such that 
(3.39) ^ strongly in L\^{n x D) as p — >■ oo. 

It follows immediately from (13.271) and (|3.15p that 

(3.40a) M5 /3|'(7/^(P)) ^ Af3 /3|'(V^) strongly in L'' (17 x £>) as p ^ cx), 

for all r e [1, oo) and, for i = 1, . . . , /C, 

(3.40b) CiiM^^P'^) ^ C^{M^) strongly in as p ^ cx). 

In order to prove (i) above, we need to show that 

(3.41) 77<P) := G(V^(P)) -> G(V') strongly in L\^{n x D) as p ^ oo. 

We have from the definition of G (see (|3.37al b')') that, for all p > 0, 
(3.42a) a{ffP\(p)^ea{vJP\pfi^^P^))m G X, 

where v^p^ £ V satisfies 

(3.42b) b{v^P'>,w)^ebi$^P'>)iw) ywev. 

Choosing ip — rfp") in p.42ap yields, on noting the simple identity 
(3.43) 2{si-S2)si^sl + {si-S2f -si Vsi,S2e]R, 

dsn, (j^TMIl and (p:?7)l that, for all p > 0, 



M 
nxD 



dg da; 



(3.44) < I M\i^^l^\^dqdx + C{L,\,a^^)At [ \V ^ v^p'^ \^ dx . 

JnxD ' ^ ~ Jn ~ ~ ~ 

Choosing w = v^P'^ in (|3.42bp . and noting ^A3^ . (jS^O]) . ([XTS]) and (IXSOj) yields, for all p > 0, 
that 



(3.45) 



< f K-'Pd.T + CAi||r||^ +CAt / M|^Wpdgdx<C. 



Combining p.44p and p.45p , we have for all p > that 

(3.46) 1157(f) II ^ + \\v(P^\\min) < C{L, {At)-') . 

It follows from (|3.46l) . (|3.9p and the compactness of the embedding (j3.11bp that there exists a 
subsequence {{rf P''\ v'^P''^)}p^>o and functions rj & X and v &V such that, as pk 00, 



(3.47a) 




^ 9j 




weakly in L" (Q.; LIj(D)) , 


(3.47b) 


M5 V^^'-P"^ - 


-> 




weakly in L^(ri x D), 


(3.47c) 


M5 VqrfP"^ - 


-> 




weakly in L^(f2 x _D), 


(3.47d) 




-> ry 




strongly in i|/(f7 x D), 


(3.47e) 


y(Pk) _ 


-> u 




weakly in H'{n); 



where s e [1, 00) if d = 2 or s g [1,6] if d = 3. We deduce from (|3.42bl) . (|32Sib), (ICTTel) and 
(I3.40b[) that the functions v &V and e X satisfy 

(3.48) Hu^w) =4(V^)(w) VweF. 
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It follows from p.42ap . p.33a[ b). (|3.47af -e) and p.40al) that ry, -0 e X and u e satisfy 
(3.49) a(77,^)=4(«,/3i^W)(^) V^eC°°(IHrD). 

Then, noting that a(-,-) is a bounded bilinear functional on X x X, that £a{v, f3g{tp))(-) is a 
bounded linear functional on X, and recalling p.Sp . we deduce that p.49p holds for all ^ € X. 
Combining this X version of (|3.49l) and (I3.48p . we have that fj = G{'iIj) £ X . Therefore the whole 
sequence 

^^(P) = G(?A(p)) ^ G(^) strongly in Llf{n x £»), 

as p — ?> oo, and so (i) holds. 

Since the embedding X ^ -^mI^^ x -D) is compact, it follows that (ii) holds. It therefore remains 
to show that (iii) holds. 

As regards (iii), "0 = kG(V') implies that {v,ip} eV x X satisfies 

(3.50a) b{v,w) = IbWiw) Vw G V, 

(3.50b) aiip,ip) = Kea{v,(3f{^)){(p) V^eX. 

Choosing u; = u in (I3.50ap yields, similarly to (|3.45l) . that 

dx + Ati' / iW.^vl'^dx 



\vI + \v-U:;^'\'-\u:r'\' 



(3.51) 



At 



Choosing if = [-^i"]'(V') in (|3.50bp and noting the convexity of Tg, p.27p and that v is divergence- 
free, yield 



dg dx 



At ^ ^ 

4 = 1 j = l 



SlxD 

a: 



M V,^7A-V,,([J-f]'(^^))dg dx 



< 



(3.52) 



/vAt / M (7{v)qi -^J q^%l)dqdx 

^ r 

AvAt V / Gj(MV') : cr(u)dx. 



where in the transition to the final inequality we applied (j3.16p with B := q;{v) (on account of 
it being independent of the variable q), together with the fact that tv{q{v)) ~ Wx ■ = 0, and 

recalled (|2.4al) . Next, on noting p. 271) and that u"^^ e V, it follows that 



Mu^2V-y,([j-i^]'(v))dg dx 



-1 ^ 



nxD 



Yx'ip dx 



My^^i •y,(G|'(7/^))dg dx = 0. 



(3.53) 

where Gf G C^'^R) is defined by 



(3.54) 



J_„2 I 5^ 
2(5* ^ 2 

* 2 

2L^ 



if S < (5, 

if s e [6,L], 
if s > L; 
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and so [Gf]'{s) = s//3|'(s). Combining (P3T|) and and noting (15351) . (l3^25bl) and 

yields that 



(3.55) 



dx + nAtv / iV^wl^da; + A: / M J^l" (tp) dq dx 
'n ~ ~ ~ JnxD ~ ~ 



kL^^ At 



M 



nxD 



dq dx 



<KAt{r,v)H.(^) + - / |<-^|^dx + fc / M^^i>^ijl^-^')dqdx 



n 

HiAtv f ,„ ,9 , K At ,, ^„,,o 



(Hi(i^))' 



It is easy to show that -Fsis) is nonnegative for all s S M, with J^g{l) — 0. Furthermore, for 
any k € (0, 1], T^iKs) < T^{s) if s < or 1 < k s, and also T^iKs) < T^{0) <lifO<KS<l. 
Thus we deduce that 



(3.56) 



Hence, the bounds (|3.55p and p.56p . on noting p.26p . give rise to the desired bound (|3.38p with 
C* dependent only on S, L, k, ag, f and ^^^^ ■ Therefore (iii) holds, and so G has a fixed point, 
proving existence of a solution to p.SGai bV □ 



Choosing w = g in p.36ap and (p = [^s]'('4'^ ^ s) (|3.36b[) . and combining, then yields, 
as in (|3.55p . with C(L) a positive constant, independent of 5 and At, 



(3.57) 



\uIl,5? + \uIl,8-u::^? 



dx + k M J'^{^PlL g)dqdx 



+ At 

k L^^ flo 



^ / \V.ulL^,\'dx + kL-'e 



MlV.^l^J^dqdx 



2A 



nxD 



M\V,^Pl\s\^dqdx 



< C{L). 



-Td.T + fc / MFf{re:L )dqdx 



nxD 



We are now ready to pass to the limit 5 — !• 0+, to deduce the existence of a solution {(w"^,, V'"i,)}nLi 
to (Pf i), with G V and V^"l € ^ H Z2, n = 1, . . . , iV. 



Lemma 3.3. There exists a subsequence (not indicated) of {{u^^ g,ip" g)}s>o, and functions 



ie,L 



g V and -ip^^ e X n Z2, n G {1, . . . , N}, such that, as S ^ 0+, 



(3.58a) 



weakly in V, 



(3.58b) 



strongly in L'^(il), 
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where r e [1, oo) if d — 2 and r G [1,6) if d ~ 3; and 



(3.59a) $Il.s weakly in L^{n x D), 

(3.59b) M5 VgV^"i.5 ^ Af5 V^V'^.L weakly in {Q x D) , 

(3.59c) M^W^i^^L s^ M^W^i^^L weakly in L^{nx D), 

(3.59d) $Il s ipi^ strongly m L^{VLx D), 

(3.59e) M^/3,^(V^^^i_,)^MM^C,L) strongly m x D) , 

for all s € [2, oo) and, for i = 1, . . . , K , 

(3.59f) C,{M^Il^s) CriMi^lL) strongly in L^{n). 



Further, (u"^, V'"i) solves i 3. 21 al b) torn — 1, . . . ,N. Hence there exists a solution {(u"^, '(/'"l)}^ 
to (Pf l), with yJ^j^ e V and ij^^^^ E X D Z2 for all n = 1, . . . , N . 

Proof. The weak convergence results p.58al) . p.59ap and that > a.e. on 51 x fohow 

immediately from the first two bounds on the left-hand side of p.57p . on noting p.26p . The 
strong convergence (j3.58bp for y"^ g follows from p.58ap . on noting that V C i?o(51) is compactly 
embedded in L^{^) for the stated values of r. 

It follows immediately from the bound on the fifth term on the left-hand side of p. 571) that 
p.59bp holds for some limit g G L^{flxD), which we need to identify. However, for any r] € ($7 x 

D), it follows from (jl.5p and the compact support of 77 on D that [V^ • (Af 5 ry) ]/M^ g L^(ri x D), 
and hence the above convergence implies, noting p.59ap . that 

r /• , ^ V, • (Mi 77) 

/ g-r]dqdx^- M ^ ^jj^ ^ g ~ dq dx 

JnxD- ~ ~ ~ JnxD M2 ^ ~ 

r V, -(Afir,) ^ 

(3.60) -> - / 2 d(7 da; = - / V, • (M^ 77) dq dx 

JnxD Af2 ^ ~ JnxD ' ~ ~ ~ ~ 

as 6 — >■ 0+. Equivalently, on dividing and multiplying by Afs under the integral sign in the 
left-most term in p.60p . we have that 

/ M-^g ■ A'Hrjdqdx^ ~ I i^^^Y g ■ (M^) dq dx VrjeCHrixD). 
JnxD ~ ~ ~ JnxD ' ~ ~ 

Observe that 77 e Cl{ft x D) 1-^ M^rj £ Cl{fl x £)) is a bijection of CK^l x D) onto itself; thus, 
the equality above is equivalent to 

/ M-ig-xdqdx^ - f ^^^^ , ■ x) dq dx Vx £ Cj(r! x i?). 
JnxD ~ ~ ~ JnxD 

Since Cff (OxD) dCliflxD), the last identity also holds for all rj £ C^iflxD). As £ L°°{D) 
andAf"5 e L^^{D), it follows that Af-^g e Lf^^{nxD) and^A^^^ G L^^^{nxD). By identification 
of a locally integrable function with a distribution we deduce that M^^g is the distributional 
gradient of "0"^ 'w.r.t. q: 

M-^g ^Yq^e.L inV'inxD). 
As M-ig e Li^^i^ X D), whereby also Y qi^lb *^ Jd\oA^ x D), it follows that 

g = M'^Y,i^l^£Ll^{SlxD). 
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However, the left-hand side belongs to L^{^ x D), which then implies that the right-hand side 
also belongs to L^{^ x D). Thus we have shown that 

(3.61) g = M3y,V'"L e L^{^ X D), 

and hence the desired result (j3.59bl) . as required. A similar argument proves (|3.59c[) on noting 
(I3.59ap . and the fourth bound in ((XFT)) . 

The strong convergence result (j3.59dl) for "iA^L 5 follows directly from (I3.59a| -c) and (I3.11bp . 
Finally, (IXsgel f) follow from (|339dl) . (IX?7l) . and (1X151) . 

It follows from (I535al b). (l339bl -f). (IXMal b). (IH^Hal b). (1535]) and dM]) that we may pass to 
the limit (5 ^ 0+ in (|3.36a[ b') to obtain that (u^x^'^'Il) ^ F x ^ with -0^'^ > a.e. on r2 x Z? 
solves (IX^ . (IS32]); i.e. it solves (l3.21al bV 

Next we prove the integral constraint on "0"^. First, for m = n — 1, n, let 

(3.62) p^^{x):^ j M{q)i^^^{x,q)dq, x & ^. 

J D ~ ~ ~ ~ 

For TO — n — 1, n, as € X, we deduce from the Cauchy-Schwarz inequality and Fubini's 
theorem that p"^ e i7^(f2) and p"^^ € L^(f2). We introduce also the following closed linear 
subspace of X = HIj{^} x _D): 

(3.63) H\n)(g)l{D):^l^$£Hlj{nxD) : i^{-,q*)^lf{-,q**) Vg*, g** G i:>| . 

Then, on choosing (p — ip £ H^{il) (E) !(£*) in (I3.21b[) . we deduce from p.62p and Fubini's theorem 
that, for all ip £ H^i^l), 

with 

(3.65) 0<P°l:= [ M0°^dg= / M/3^(0")dg< / Af0°dg<l, a.e. on 17. 



By introducing the function z™^ := 1 — p™^, to = n — 1, n, and noting that z"^ € i/^(ri) and 
z"^^ e L^(ri), we deduce from p.64p . and as u"^^ is divergence-free on f2 with zero trace on d^, 
that 

(3.66) / ^"-^ . ^"'^ ifdx+ [ [e V.zIl - f-j^ zlj\ • V,^ dj = 0, 

for all if £ H^{il). Let us now define by [x]± :— ^ {x i: \x\) the positive and negative parts, [a;]-|- 
and of a real number x, respectively. As "0"^^ £ Z2, we then have that [z^^]- = a.e. on 
fi. Taking ip = [z"^]- as a test function in p.66p . noting that this is a legitimate choice since 
€ H^{^), decomposing z^j^, m = n — l,n, into their positive and negative parts, and 
noting that Jt"^^ is divergence-free on fl and has zero trace on dQ, we deduce that 

\\[zl^U\' + Ats\\VA4\U\'^0, 

where || • || denotes the L'^{n) norm. Hence, [z^ i]- = a.e. on ft. In other words, z"^ > a.e. 
on ri, which then gives that p"^ < 1 a.e. on fi, i.e. 0"^ £ Z2 as required. As (u"^,V'e,L) & 
V X Z2, performing the above existence proof at each time level tn, n = 1, . . . , N , yields a solution 

{fe",L:^?,L)}n=lt0(Pf^i). □ 
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4. Entropy estimates 

Next, we derive bounds on the solution of (Pf^^), independent of L. Our starting point is 
Lemma 13.31 concerning the existence of a solution to the problem (P^^). The model (P^^) 
includes 'microscopic cut-off' in the drag term of the Fokker-Planck equation, where i > 1 is a 
(fixed, but otherwise arbitrary,) cut-off parameter. Our ultimate objective is to pass to the limits 
L — > oo and At — > 0+ in the model (Pf^^), with L and At linked by the condition At = o{L~^), as 
L — >■ oo. To that end, we need to develop various bounds on sequences of weak solutions of (P^^) 
that are uniform in the cut-off parameter L and thus permit the extraction of weakly convergent 
subsequences, as L — >■ oo, through the use of a weak compactness argument. The derivation of 
such bounds, based on the use of the relative entropy associated with the Maxwellian M, is our 
main task in this section. 

Let us introduce the following definitions, in line with p. 221) : 



(4.1a) 
(4.1b) 



.At 



t tn 



At,+ 
'■e,L 



At 

(•,t) :=^"(.), 



At,- 



tr, 



t 



At 



(•), t e {tn-l,t„ 



n — 1, 



1,...,A^, 
..,iV. 



At(,±) 



We shall adopt ^ 



as a collective symbol for u 



At 



At,± 



The corresponding notations ip 



At(,±) 



-0^1^, and L defined analogously; recall p.lSp and (I3.20p . 



We note for future reference that 



(4.2) 



At 



At,± /, ,±N "Hex 



dt 



t e (i„„i,t„), n^l,...,N, 



where '■=tn and t^ :~ t„_i, with an analogous relationship in the case of 

Using the above notation, p.21a|) summed for n = 1, . . 
reminiscent of a weak formulation of (jl.lal dl: find w^^*''^'' 



. , N can be restated in a form that is 
e y, t G (0,r], such that 



Jn 



dt 



; da; dt 



Jn 



(4.3) 



At, 

''e,L 

K 



V7 N Af,+ 



V7 At,+ 



Va; W 



dx dt 



{f^''+,w)H^^^n) -kJ2 f C.,{Mi^f'£+) : ^^wdx 



dt. 



for ah w e L^(0, T; V), subject to the initial condition w^^(-, 0) = e V. 

Analogously, (after a minor re-ordering of terms on the left-hand side for presentational reasons,) 
(I3.21bp summed through n = 1, . . . , N can be restated as follows: find V'^^ '^^(t) G Z2, t G (0, T], 
such that 



M 



T 



dt 



(fi dq dx dt 



1 

2A 



/ M \e\7oci^f*£~^ - uf]^ -V^ipdqdxdt 
JnxD ~ ' ~ ' ' J ~ ^ ~ 

/ / Vg, V^f L + • V,, ^ dq dx dt 



(4.4) 



Jo JnxD 



i=i i=i 

K r 



/ At, + N 



/3^(V^fl'+) • V,,^dgda:dt = 0, 



for ah e L^{0,T;X), subject to the initial condition ■0^[(-,-,O) = /3^(-0"(-, •)) ^ Z2. We 
emphasize that (|4.3p and (|4.4p are an equivalent restatement of problem (P^^), for which existence 
of a solution has been established (cf. Lemma [373]) . 

1, . . . , can be 



Similarly, with analogous notation for {p"^}^^0' P-641) summed for n 

Af(, 



restated as follows: Given uf^*^'^\t) eY,t€ (0,T], solving (gS]), find pff^\t) e IC :^ {t] e 
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H^{n) : e [0, 1] a.e. on Q.}, t e (0,r], such that 

Jo Jn ci^ ~ Jo Jn ~ ' ~ ' ' J ~ 
(4.5) '^^eL\0,T;H\n)), 

subject to the initial condition pf^i(-,0) = J^M{q)l3^{^(;°{-,q))dq e JC: cf. ([^^ and recall that 

■0°^ = ^■'^{■ip^). Once again, on noting (|3.62p and (|3.64p . we have established the existence of a 
solution to (|4.5p and that 



(4.6) ptLi^^t) = / M{q)i,'^'^^^\x,q,t)dq for a.e. {x,t) e 17 x (0,r). 

jd ~ ' ~ ~ ~ 

In conjunction with /3^, defined by (|f .f 2p . we consider the following cut-off version of the 
entropy function F : s ^ R>o H> T[s) ~ s(log s — 1) + 1 € M>o: 



(4.7) 
Note that 

(4.8) (J-^)'(.) 
and 



s(logs-f) + f, 0<s<L, 
' s(logL-l) + l, L< s. 



2L 



(4.9) {F^)"{s) 



log s, < s < L, 

f+logi-1, L<s, 

' i, < s < L, 



i L<s. 



(4.11) (■^'^)"(s) > -F"(s) = s-\ s e 



Hence, 

(4.10) /3^(s) =min(s,i) = [(J-^)"(s)]-\ s e M>o, 

with the convention l/oo := when s = 0, and 

. „ — 1 „ ^ TTD 

''>0- 

We shall also require the following inequality, relating J^^ to T: 
(4.12) F^{s)>F{s), seR>o. 

For < s < 1, (|4.12l) trivially holds, with equality. For s > 1, it follows from (|4.11[) . with s replaced 
by a dummy variable cr, after integrating twice over a e [l,s], and noting that (J"^)'(l) P [I) 
and (J"^)(l) = J"(l). 

4.1. L- independent bounds on the spatial derivatives. We are now ready to embark on 
the derivation of the required bounds, uniform in the cut-off parameter L, on norms of u^^£^ 

and V'^L^- As far as u^^'^ is concerned, this is a relatively straightforward exercise. We select 

W ~ X[o,i\ u'^l'^ test function in (j4.3p . with t chosen as tn, n S {1, . . . ,N}, and X[o,t] denoting 
the characteristic function of the interval [0,^]. We then deduce, with t = i„ and noting p.lSp . 
that 



Jo 

<\\iiof+i rii,/^*^+(.s)ii^^.(,,)),d. 

(4.13) -2k [ [ Miq)Tq,qJU^{^\q,\^)^fl:+ -.V.uf'l+dqdxds, 

Jo JnxD ~ i=i ~ ~ ~ ' ~ 

where, again, || • || denotes the norm over Q. We intentionally did not bound the final term 
on the right-hand side of ()4.13p . As we shall see in what follows, this simple trick will prove 
helpful: our bounds on ijj^ ^'^ below will furnish an identical term with the opposite sign, so then 
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by combining the bounds on u^^^ and "0^^'^ this pair of, otherwise dangerous, terms wih be 
removed. This fortuitous cancellation reflects the balance of total energy in the system. 

Having dealt with Jt^^'^, we now embark on the less straightforward task of deriving bounds 

on norms of "0^1' that are uniform in the cut-off parameter L. The appropriate choice of test 

function in (|4.4p for this purpose is ^ = X[o.t\ (•^^)'('0^l^) with t = tn, n e {1, . . . , N}; this can 
be seen by noting that with such a (p, at least formally, the final term on the left-hand side of (|4.4I) 
can be manipulated to become identical to the final term in (I4.13p . but with the opposite sign. 
While Lemma [3.31 guarantees that 0^ ^'^(•, •, i) belongs to Z2 for all t € [0,T], and is therefore 
nonnegative a.e. on il x D x [0,T], there is unfortunately no reason why "0^^!'^ should be strictly 
positive on £7 X L) X [0, T], and therefore the expression (J^^)' {'tp'^*£~^) may in general be undefined; 
the same is true of (J"^)"(0^^['+), which also appears in the algebraic manipulations. We shall 
circumvent this problem by working with (T^)' {ipf'*£^ +a) instead of {J^'")' {ipf*£^), where a > 0; 
since '0^1'''' is known to be nonnegative from Lemma (J^^)'(0^^''^ -fa) and (J-^)" {^/j'^'i^^ + a) 
are well-defined. After deriving the relevant bounds, which will involve J^^(V'f^|^'^ -I- a) only, we 
shall pass to the limit a 0+, noting that, unlike {J-^)' {tpfj^'^) and (J^^)"(0^^'^), the function 
is well-defined for any nonnegative ■0^l^- Thus, we take any a e (0,1), whereby 

< a < 1 < L, and we choose (p = X[o,t] {^^)' {''I'f]]^ + ct), with t = tn, n G {1, . . . , N}, as test 
function in (|4.4p . As the calculations are quite involved, for the sake of clarity of exposition we 
shall manipulate the terms in (|4.4p one at a time and will then merge the resulting bounds on the 
individual terms with (|4.3p to obtain a single energy inequality for the pair {u^ ip^ l^)- 

We start by considering the first term in (|4.4I) . Clearly -I- a) is twice continuously differ- 
entiable on the interval (—a, 00) for any a > 0. Thus, by Taylor series expansion with remainder 
of the function 

s e [0,00) ^ F^{s + a) e [0,00), 
we have, for any c € [0,oo), that 

(s - c) (J-^)'(s + a) - J-^(.s + a) - J-^(c + a) + i(s - cf (F^^y^es + (1 - 6)c + a), 
with e e (0, 1). Hence, on noting that t e [0,T] 1-^ ^^X'^'' ''0 ^ ^ 

is piecewise linear relative to 

the partition {Q = to,ti, . . . ,tN — T} of the interval [0, T], 



ds 



nxD ' ' 



JnxD 



Af^(^^i + a) (J-^)'(^^i'+ + a) dq dx ds 



nxD ' JnxD 



MF^ (V^f {t) + a)dqdx~ / MF^ {(3^ + a) dq dx 



1 



2Ai 



JnxD 



M{F^ri9$t;^+ + (1 - e)$^l- + a) (V^^l^+ - ip^l- f dq dx ds. 



Noting from (|4.9p that {F^)"{s + a) > l/L for all s G [0,oo) and all a > 0, this then implies, 
with t = t„, n e {1, . . . , N}, that 

Ti > / MF^{i)fl+{t)+a)dqdx- I Af J'^(/3^(?A") + a) dg dx 
JnxD ~ JnxD 

ft 



1 



At,+ 7At, 



(4.14) +2A^ / / (Cr- err dqd^ds 



JnxD 



The denominator in the prefactor of the last integral motivates us to link At to L so that At L — 
0(1) as At— >-0+ (or, equivalently. At — o{L~^) as L ^ 00), in order to drive the integral multiplied 
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by the prefactor to in the hmit of L — > oo, once the product of the two has been bounded above 
by a constant, independent of L. 

Next we consider the second term in (|4.4I) . using repeatedly that V^; • uf*'^' — and that uf^'^ 
has zero trace on dH. 

T2 



M 



JnxD 
t 



At,- 



At,- TAt.- 



= e 



JnxD 



MYAi^^L^ + a) • y.(^^)'(v5^^i'+ + a) dq dxds 



'0 JnxD 

where in the last line we subtracted in the form of 



Myfl;-(^f['+ + a) ■ y,(^^)'(^fi'+ + a) dq dx ds, 



a 

Iq JnxD 
Hence, similarly to p. 531) . 



Mufl'' ■ YAJ^'^Yii'Tl'^ + a) dg dx ds = 0. 



-L\//7:At,+ 



To 



A/(J-^)"(Vfi'+ + a)\YA^t.L^ + dq dx ds 



JnxD 



JnxD 



M(J-^)"(V?i'+ + a)ly.Cr I' dg dxds 



At,+ |2 , 



JnxD 



JnxD 



JnxD 



[ 1/i 

A/(^^)"(V^fi'+ + a)|y,^fi'+p dg dxds 



At.H 
At.H 



a < i 
a > i 



dg dx ds 



Muf'£^ ■ YxiG'^ii'fL^ + a)] dq dx ds 



JnxD 



where denotes the (locally Lipschitz continuous) function defined on R by s — if s < L and 
Y^s^ otherwise. On noting that the integral involving vanishes, (|4.1ip yields the lower bound 

(4.15) T2>e f f M(^f['++a)-i|y,(V^fi'++a)pdgdxds. 

Jo JnxD 

Next, we consider the third term in (|4.4p . Thanks to (|2.6I) we have, again with t — tn and 
n e {1,...,A^}: 



1 

Ts := — 



/ / E E y-j.^s'^L^ • y«,x[o,t](-^^)'(v^^'L+ + «) dg dxds 
/* / Af (j-^)"(v^^i'+ + ") E E y?.^^e''I''' • y.-v^s'^r dg dxds 

Jo JnxD 

>Txf I M{F^ni.f^^ + a)j^\Yjf;^+\^ dgdxds 
^-^ Jo JnxD ~{ 



1 

2A 



(4.16) 



ao 
2A 



JnxD 



M (-F^)"(V^fi'+ + a) |y,^^fl^+|' dg dxds. 
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We are now ready to consider the final term in (|4.4p . with t — tn, n E {1, . . . , N}: 



T4 := - 



(4.17) 



/ / MX:[£(MeT)9.]/3^(^e,L 
Jo JnxD j^]^ 



Jo JnxD j^-^ 



xUe,L We,£ > 



{TL)"{^f'i^ + a) dq dx ds 



Jo JnxD j^j^ 



(V'e,L ) Y qiip^i,'^ dq dx ds 



a 



f I M V[(y,Mfl'+)gd ■y,,v^fi'+ dgdxds 

Jo JnxD j^j^ 

K 



JnxD 



4=1 

if 



1 



y^i V'^L^ dq dx ds 



Jo JnxD ~ ~ ~ 



a: 



JnxD 



1=1 



Yqii^^L^ dq da; ds, 



where in the transition to the final equality we applied p.l6p with B Wx £ (on account of 
it being independent of the variable q), together with the fact that tr (y^; uf^£^) = Yx ■ y^i'^ = 0. 
Summing (|4.14l) - (|4.16p and (|4.17p yields, with t = tn and n G {1, . . . , A^}, the following inequality: 

/ MJ-^(v7fi'+(i)+a)dQdx+ -1- /* / Mii^f'^-i^f'j^fdqdxds 
JnxD ~ ^^i^Jo JnxD 



M 



f I ' ' '?At,+ , 

JnxD ip^j] +a 

Af (J-^)"(V^^i'+ + a) |y,^^l^+P dgdxds 



dq dx ds 



a-o 
2A 



< 



Jo JnxD 

7L I rL 1 2o 



JnxD 



M F^{(3^($^)+a)dq dx 



- f I MEg.g^c/K^kn^fr :y.yfrdgdxds 

Jo JnxD .-I ~ ~ 



1=1 
K 



Jo JnxD j^]^ '- 



1 - 



• yijiV'f^L^ dxdi 



(4.18) 



Comparing (|4.18p with (|4.13p we see that after multiplying (|4.18p by 2k and adding the resulting 
inequality to (|4.13p the final term in (j4.13p is cancelled by 2fc times the second term on the 
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right-hand side of (I4.18p . Hence, for any t — tn, with n e {1, . . . , N}, we deduce that 



(4.19) 



At,+ At,-||2 



ds + iy f ||y,yfi'+(s)f ds 

JQ 



AtL 



JnxD 



+ 2k£ 



M 



At. 



"At,+ |2 



JnxD ip T + a 



■ dq dx ds 



[ Mi^ni^t:L^ + a)\YJt^^\'dqdxds 
Jo JnxD 



<\\uor + - ||r'-+(.s)|| ^,(f,)),d,s + 2fc / MT^{(3^{r)+a)dqdx 

^ Jg ~ y " JnxD 



2k f f mY^V^^V^Tl')'! 
Jo JnxD j^]^ 



Yqi^-'til'^ d(7 dx ds. 



It remains to bound the last term on the right-hand side of (j4.19p . Noting that is Lipschitz 
continuous, with Lipschitz constant equal to 1, and f3^ {s + a) > a for s > 0, we have that 



< 1 



< 



'e,L ) ^ \ ^/a when "0^^^ < i, 
- I when ^f^j]^ > L. 



With this bound we now focus our attention on the last term in the inequality (|4.19p . Let 
b := . . . , bn) and b := \b\i := bi + ■ ■ ■ + bx; as \qi\ < y/hi, i — 1, . . . ,K, we have that |g| < Vb 
for all q E D. For t — tn, n E {1, . . . , N}, we then have that 



-2fc 



K 



JnxD 



z=l 



<2k / M\V,U^l^\\q\ 
Jo JnxD 



e,L 



y g. dq dx ds 



< 2k 



MlY.y^rikl 



JnxD 



At.+ ii„i J V" when ■0^^^' < L, 
when i'f*£~^ > L 



M ly, yfl'+l ly^V^fl'+l dq dx ds 



< 2kVab 

lo JnxD 

2kV^b f I ly.ufri ( / ly.cn ) dxds 

Jo Jn ~ \Jd 



< 2kVab I I |y, ufi'+l ( / Af(J-^)"(d'+ + a) |y,^fi^+|^ dq ] dxds 
Jn ' 



-L\///7At,+ 
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JnxD 



ao k 



(4.20) +a 



Jo Jn 



Substitution of (|4:20| into (|4l9l) and use of (|41T|) to bound + a) from below by 

P'i^f^j]^ + a) = + and Bl^ to bound + a) by + a) from below 

finally yields, for alH = n e {1, . . . , iV}, that — 

\\nfr{t)f + ^j\\uf:i'-uf^frds + .j\\V.u^^^ 

+ 2k f M ^(V^f i'+ (t) +a)dqdx+-^ f [ M {i^f'j^ - )2 dg dx ds 

Jo JnxD ijj^ £ + a ~ Jo JnxD -tp^ £ + a 

<hor + - r||/^*'-*'(5)||^^i(o)vd.s + 2fc / MJ-^(/3^(V^«)+a)d9dj 
^ Ja ~ ^ " JnxD 

t 



(4.21) / lly.MfrWiPd. 

flo Jo ~ 



The only restriction we have imposed on a so far is that it belongs to the open interval (0, 1); let 
us now restrict the range of a further by demanding that, in fact, 

(4.22) 0<«<min(l,^). 

Then, the last term on the right-hand side of (|4.2ip can be absorbed into the third term on the 
left-hand side, giving, for t — tn and n G {1, . . . , N}, 

+ 2k f M J-(V^f i'+(t) + a)dqdx + -^ f f M (^f - dq dx ds 

JnxD ~ ^t^Ja JnxD 

M%|^d,d,d.+ f|^/7 M^|i^d,d,d, 

Jo JnxD tp^ j] + a ~ 2A Jo JnxD ip^_£ +a 

(4.23) <\\uof + l f \\f^t'+{s)\\^ , ds + 2k f Af J-^(/3^(^0)+a)d(zd.T. 

Jo " \ ov " JnxD 

We now focus our attention on the final integral on the right-hand side of (|4.23p : 



T5(a) := / A/J"^(/3^(V/') + a)dgd;E= / M F'^ [P'^ {il)") + a) dq dx, 

JnxD ~ J'2I,l,c'J<3l,c 

where 

Sli^a {ix,qJenxD : 0<p'^iip'^ix,q))<L-a}, 

Si.a := {ix,q)eflxD : L-a<p^i$°ix,q))<L}. 
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We begin by noting that 

MF^{l3^{i>°) + a) dq Ax 
For the integral over *8l.q we have 

MJ-^(/3^(V^*') + a)dg dx 
"(/3^(V^0)+a)2-L2 



Si 



Si 



< 



M 



M 



+ (/3^(V/)+a)(logL-l) + l 



a (l 



2L 



dq dx 



+ (/3^(^°) + aHlogiP'^ir) + a) - 1) + 1 



Lf20\ 



2L 



dq dx 



2LJ 



'Sl 



M dq dx 



MJ"(/3^«) + a)dg dx 



AfJ"(^^(V'") + a)dg dx. 



Thus we have shown that 



(4.24) 



TsH < -a|r!| 



MT{/3^{^P°) + a) dq dx. 



Now, there are two possibihties: 

Case 1. If + a < 1, then < < 1 - a. Since L > 1 it follows that < /3^(s) < 1 

if, and only if, /3^(s) — s. Thus we deduce that in this case — ip'^, and therefore 

< F{I3^{$°) + a) = Tii}^ + a). 

Case 2. Alternatively, if 1 < /3-^('i/;°) + a, then, on noting that /3^(s) < s for all s G [0, oo), it follows 
that 1 < (■;/''') + a < ^/^^ + a. However the function F is strictly monotonic increasing 
on the interval [1, oo), which then implies that = < F{I3^{%Ij^) + a) < J"(?A° + a). 

The conclusion we draw is that, either way, < J^(/3^('0°) + a) < + a). Hence, 



(4.25) 



Tz{a)<-a\Q.\+ I M + a) dq dx. 



nxD 



Substituting (|4.25l) into (|4.23l) thus yields, for t = t„ and n e {1, . . . ,N}, 



At,+ 



{t)f 



At Jo 



1 At,+ At.-||2 1 



2\bk 
ao 



iiy.u^l;+(s)rds 



-2fc 



M F{-4)f^£^{t) + a) dq dx + 



SlxD 



-2fc£ 

/o JnxD 

aok 
"2A" 



AtL 
dq dx ds 



JnxD 



Mii^^}+ -ipf}:')'dqdxds 



JnxD 



M ,' dq dx ds 



At,+ 



+ a 



< 



1 



iiyoii^+^ iir'^(.)ii^H„^ 



(4.26) 



3afc|17|+2fc/ MT{i^° + a)dqdx. 



nxD 



The key observation at this point is that the right-hand side of (|4.26p is completely independent 
of the cut-off parameter L. 
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We shall tidy up the bound (|4.26l) by passing to the limit a — > 0+. The first a-dependent 
term on the right-hand side of (|4.26p trivially converges to as a — )■ 0+; concerning the second 
Qf-dependent term, Lebesgue's dominated convergence theorem implies that 

lim / M + a)Aqdx= I M dq dx. 

°'^°+JnxD ~ ~ JnxD 

Similarly, we can easily pass to the limit on the left-hand side of (|4.26l) . By applying Fatou's 
lemma to the fourth, sixth and seventh term on the left-hand side of (|4.26l) we get, for t — i„, 
nG {1,...,A^}, that 



lim inia^Q 



lim infc(_s.o , 



M 



(IxD 



M^(^f['+(0+a) > / MT{^P^l^{t))dq dx, 



flxD 
t 



I I T'At.+ , 

10 JnxD ip^j^ + a 



dq dx ds > 



M 



JnxD 



dq dx ds 



lim infc,^04 



M dgdxds > 



10 JnxD £ + a 



JnxD 



M 



M|yWV^fI'+|'d(?dxds, 



JnxD 



TAt, 



^ — dq dx ds 



JnxD 



Thus, after passage to the limit a 0_|_, on recalling p.20p . we have, for alH — n G {!,..., N}, 
that 



At.+ At.-||2 
r - r 



ds + v I \\^xuf'{^{s)fds 
Jo ~ ~ ' 



'2k 



Af J'(V'fi'+(t))d(7 dx + 



nxD 



AtL 



M{^f'£+ -^f'£-fdq dxds 



JnxD 



10 JnxD 
2aofc /"* 



A^ly-V^e^L^l'dg Axds 



JnxD 



(4.27) < llyoll' + f- / ||/''*'+(s)||^„i(j,)),ds-H2fc / MF{^'^)dqdx 



A/lVg-y/V-f^'+l^ dgdxds 



(Hhm' 



nxD 



< 



hor + - \\f{s)\\'u^.yds + 2k MTii^o)dqdx^:[B{uoJ,i^o)] 

-In ~ JnxD 



(4.28) 



where, in the last line, we used (|3:20)) to bound the third term in g^, and that t e [0,r] 
together with the definition (|3.22p of f^^'^ to bound the second term. 

We select tp = X[o,t] pf^l^ as test function in (|4.5p . with t chosen as t„ and n e {!,..., N}. 
Then, similarly to (|4.13l) . we deduce, with t — tn, that 



ii/'.''S:"(»)-f.';rwii'<i»+2s / iivrf^rwr d. 



At, 



At,+ / 



(4.29) 



< 



/3^(^0)dg 



D 



< \n\ 



where we have noted ([3331), (13301) and that I3^{tp°) E Z; 
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4.2. L- independent bounds on the time-derivatives. Next, we derive L-independent bounds 

itl ^eX and pfX- 



on the time-derivatives of tire functions u^^, i^f]^ and pf\^- We begin by bounding the time- 
derivative of '(/'f^^ using (|4.28l) : we shall then bound the time-derivatives of and uf\^ in a 



similar manner. 



4.2.1. L -independent hound on the time- derivative of ijj^^^. It follows from (|4.4p that 



(4.30) 



, at 

JnxD oi 



(pdq dx dt 







< 






JO Jn 



M ■ ^x^dq dx dt 



J(lxD 
T 



A'l uf'^£ V'f^i'^ ■ Va;<^ dq dxdt 



-t- 



/ / E E ^^i' ^ • V,. ^ dg dx di 

^'^Jo JnxD i=i j=i ~ ~ ^ ~ 

/ / 

Jo JnxD ._i 



/3^(V^fi'+) • Vg.(?dqda;dt 



S1 + S2 + S3 + S4 yipeL\0,T;X) 



We proceed to bound each of the terms Si, ... , S4, bearing in mind (cf. the last sentence in the 
statement of Lemma 13.31) that 



(4.31a) 



> a.e. on r2 X D X [0, T], 



M{q)dq = l, 



D 



(4.31b) 0</ M{q)il,f*j^^{x,q,t)dq <1 ior a.e. {x,t) e Q x D. 

Jd ~ ' ~ ~ 

We shall use throughout the rest of this section test functions ip such that 
(4.32) £ L^{0,T;H^{n;L°°{D))n L^{n;W^^°°{D))). 

We begin by considering Si; noting (I4.31al b') and (I4.28p . we have that 



Si = 2e 



< 2e 



JnxD 

T 



M V^Ti'^ y Vv^fi^ • y^-^dg dxdt 



Jn 



M 



D 



2 Y 

dg WYx'fh'^iD] 



dx dt 



M 



JnxD 



dqdxdtj (j^Jjy^^\ 



L=°(D)dg; dt\ 



Hence, by with t^t^ =T 

(4.33) 



Si < J:^B(yo,/,^o) ( / / |iy,^|li..p)dxdt 



T 

Jn 
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Next, we consider the term S2: 



JQ Jn \Jd ~/ 

(4.34) < Cp{n)(^J^ ^ly.yf^l'-pdxdt^ ^ ^J|y.^||ioo(^)dxdtj , 

where Cp{Q) denotes the (positive) constant appearing in the Poincare inequahty < 
Cp{n) WYxufL^'W on n, with ufl~ eV d Hl{n). On recaUing g^H]), the definitions of uf^['+ 
and £~ from (|4.1bp . and noting p.lSp . we have that 



~ ' Jo ~ ' 

Ai||y,y°||2+ / iiy.yf^+fdi 

Jo ~ ' 



At,- ||2 jj. 

r dt 



"'O 



rT-At 



< \\uof+ r WY.uf^i^W'dt 
Jo 

(4.35) < (l + i)[B(wo,/>o)]'. 

Therefore, 



(4.36) S2 <Cp(r!)(l + i)*B(yo,/,^^o) ^ y^lly,^|!ioop)dxdtj . 
Alternatively, directly from the second line of (|4.34|) . we have that 

(4.37) S2 < Vr(^ess.snp,g[o,T]^l?i.^ri'da;y ^||y.^||ioo(^)dxdt 
Similarly as above, 

ess.suptg[o_r] / |M^^^~pdx = ess.suptg[o,T]lly^L"WlP 
Jo 

= max (\\uof, ess.suptg (0 1 1 uf^£+ {t)\\^'j 
< max(||uo|P,ess.suptg[o_r]||M^^+(Olp) 

(4.38) < [B(yo,/,V^o)]'. 

Combining (gSH), (^37)) and we have that 

(4.39) S2 <min(Cp(f]) (l + i) ^ Vt) B(uo, /, V^o) U ^ ||y.^||io.p) dxdt j . 
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We are ready to consider S3; we have that 

T . K K 



1 

2A 



< 



/ / M2y^X^X^A,,y,^y^^.y,.^dgdxdi 

Jo Jny.D 



12. 



JnxD 



K.L 



E 



.\A\ 



/Q JnxD 
i-T 



K 

E 



At,+ 



if 



ly^,^!^ dqdxdt 



At,+ 
9j V '^e.-t 



/Q JnxD 
I-T 



TAt.- 



?At,- 



|^^|y,,(^Pj dqdxdt 
lYq'fl di? dx dt 



Jn \JD 



Mi,f'l+dq 



M 



D 



< 



\A\ /2aofc 



Thus, by (1123), 



dg da: dt 



dg \\Yq^\\L--(D}dxdt 



Jsi 



dx dt] 



(4.40) 



S, < 



1^1 



v2aofcA 



B(wo,/, V'o) 



g'f IIL~(D) dg; dt 



10 Jn 

Finahy, for term S4, recahing the notatfon b := \b\i (cf. the paragraph before (|4.20p ) together 
with the inequahty (3'"{s) < s for s E K>o, we have that 

K 



f f M\a{uf^l+)\P^i^;f}+) \q^\ \YqM dqdxdt 
JO JnxD 



< Vb 



M\<j{uf'£ + )\P'^{^f'j:^)\Yq^\ dqdxdt 



< Vb 



10 JnxD 

I-T 



< Vb 



10 Jn 

I-T 



k{uf!l^)\ [ I Af/3^(^,-r)dg) \\Yqlp\\L-iD)dxdt 



< Vb 



< Vb 



10 Jn 

I-T 



D 



\,raLf;^At^ 



I dx dt 



/ / k{ufL^)\\\Yq'P\\L°°{D)<ixdt 

Jo Jn 

1 

T \ / T \ 

I f \aiuf'^)fdxdt] if [ \\V,^\\l^^^^dxdt 
yo Jn J \Jq Jn j 

V? vio In l^^^^-^^'^l'^^;^^) Jj^Ml-iD)dxdt 



Hence, by (14281) . 
(4.41) 



S4 < W- B(uo,/,-0o) 



Jn 



\Yq'?\\L--(D) dx dt 



Upon substituting the bounds on the terms Si to S4 into (|430l) . with ip € L'^{0, T; H^{il] L°°{D))n 
L'^{Q;W^'°°{D))), and noting that the latter space is contained in L^{0,T;X) we deduce from 
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(P:^ that 



T 

JnxD 



M — TT— ^ dq dx dt 



dt 



(4.42) 



< aB(yo,/,V^o) 



T 

Jn 



dx dt 



for any (p G L'^{0,T; H^{n; L°°{D)) (1 L'^{n;W'^'°^ (D))), where C* denotes a positive constant 
(that can be computed by tracking the constants in (|4.33l) - (|4.4ip ). which depends solely on e, i^, 
Cp{n),T,\A\,ao,k,K,X,Keindb. 

We now consider the time-derivative of pf^. It follows from (|4.29p . (|4.31bp and (|4.38p 

that 



Jn 



< 



dt 



ip dx dt 



< 



T 

Jn 



T-, At,+ At,- At,+ 



l^xpf^fdt 



-eSS.SUPtg[0^T]l|P^L^IlL-(O) 



da; dt 



(4.43) 



< 



B(uo,/,'0o) 



X ^ ||V,(p||2dij 

\v,ipfdt] y^eL''{o,T-H\n)). 



4.2.2. L-independent bound on the time-derivative of u^\^. In this section we shall derive an L- 
independent bound on the time-derivative of m^^. Our starting point is (|4.3p . from which we 
deduce that 



At 



JQ. 



< 



dt 



Jn 

T 



w dx dt 



/ At.- v7 \ At,+ 



• wdxdt 



+ V 



T 

Jn 



(4.44) 



(/ ■+,«^)h„Mo) dt 
Ui + U2 + U3 U4 







E 



4=1 



V^: u^i^ : V^^ w dx dt 



C^iMi>f*£^) : V^wdxdt 
Vw e L\0,T;V). 



Jn 



On recalling from the discussion following p. 31) the definition of V^o-, we shall assume henceforth 
that 

w e L^{0,T;V„), cr>^d, a > 1. 

Clearly, L'^{0,T; V^) C ^^0,^; V). By Lemma 4.1 in Ch. 3 of Temam [49 and using (|438)) and 
(|4.28p . we have, with cr > ^d, a > 1, and c(ri, d) a constant that only depends on and d, that 



Ui < c(n,d)ess.suptg[o^T]lly^L 



(4.45) 



< c(f7,d)\/-[B(yo,/,V'o)] 



Iwll?/ dt . 
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For term U2 we have, 



At,+ ||2 



a; iie,L 



dt 



< y^B(yo,/,^o) / l|y,u;||"dt 



(4.46) 



Concerning the term U3, on noting the definition of the norm || • ||(H,5(f2))' s-^id that thanks to 
P-2^P we have ||/'^*'+||L2(o,T;(ffi(f^))') - U\\L^o,T;(Hl(n)y), it follows that 



(4.47) 



U3 < y^B{uoJ,^o) / WYxwfdt 







Before we embark on the estimation of the term U4 we observe that 

K 



U4 



k 



(4.48) 



/ / / MV;,^l'+5]Q.gTf,,/i|^^|2\ v.z^dg 

T r ^ 

Jo Jn \ J D 



da; d< 



where we used the integration-by-parts formula p.l6p to transform the expression in the square 
brackets in the first line into the expression in the square brackets in the second line. Thus we 
have that 



U4 < fc 



< k 



r [ M\Y.w\ f V|g.||y,,^fl'+| ) dqdxdt 

Jo JnxD ~ \i=l ~ ' y ~ 



M\Yxw\\q\\YJtl:^\dq dxdt 



JnxD 

T 



< kVb [ [ \Yxw\ ( I M \V yip f*£+\dq] dxdt 
Jo Jn ~ ~ \Jd ~ ^ ~J ~ 

i-T 



T 



(4.49) 



< 2k\/b 



Jn 



ly.u'l ( / M\^yJi^f'£-^\^dq\ dxdt, 



D 



where in the transition to the last line we used the Cauchy-Schwarz inequality and (j4.31bp . Hence, 
by 6251), 



U4 < 2kVh\ 



JnxD 
2Xbk ( 2aak /"^ 



Mly^VV^fi'^rdg dxdt] 



2 / „ J" 



iVx w| dxdt 



nxD 



M\Wq\l'4)fl'^? dq d^dt 



Jo 

2 / 



Wxwfdt] 







(4.50) 



< 



2\hk 



ao 



B(yo,/,V'o) / ■ 
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Collecting the bounds on the terms Ui to U4 and inserting them into (|4.44l) yields 

< C maxf[B(Mo,/>o)]',BK,/,^^o)) ( / llu'lly dt) 



■wdxdt 



Jn 



dt 







(4.51) 

for any w G L^(0, T; 1/^), cr > 5^, cr > 1, where C** denotes a positive constant (that can be 
computed by tracking the constants in (|4.45p - (|4.50l) '). which depends solely on f2, d, v, k, K, A, 
oo and b. 

5. DUBINSKIl'S COMPACTNESS THEOREM 

Having developed a collection of L- independent bounds in Sections 14.11 and we now describe 
the theoretical tool that will be used to set up a weak compactness argument using these bounds: 
Dubinskii's compactness theorem in seminormed sets. 

Let A he a linear space over the field M of real numbers, and suppose that is a subset of A 
such that 

(5.1) (V(^ e X) (Vc e M>o) opeM. 

In other words, whenever tp is contained in M, the ray through tp from the origin of the linear 
space A is also contained in M . Note in particular that while any set A4 with property (|5.ip must 
contain the zero element of the linear space A, the set M need not be closed under summation. 
The linear space A will be referred to as the ambient space for Ai. 

Suppose further that each element of a set A4 with property (|5.ip is assigned a certain real 
number, denoted [(^]ai, such that: 

(i) [p]m > 0; and [ip]M — if, and only if, tp = 0; and 

(ii) (Vc e M>o) [cp]m = c[ip]M- 

We shall then say that is a seminormed set. 

A subset S of a seminormed set Ai is said to be bounded if there exists a positive constant Kq 
such that [(p]m < Kq for all p £ B. 

A seminormed set Ai contained in a normed linear space A with norm || • is said to be 
embedded in A, and we write A4 ^ A, ii the inclusion map i : ip £ At 1-^ i{ip) — (p £ A (which is, 
by definition, injective and positively 1-homogeneous, i.e. i{cp) = ci(p) for all c £ M>o and all 
ip £ Ai) is a bounded operator, i.e. 

[BKo £ M>o) (V(^ £ M) ||i((^)|U < Ko[p>]m- 

The symbol i( ) is usually omitted from the notation i{p), and (/s G is simply identified with 
p £ A. Thus, a bounded subset of a seminormed set is also a bounded subset of the ambient 
normed linear space the seminormed set is embedded in. 

The embedding of a seminormed set Ai into a normed linear space A is said to be compact if 
from any infinite, bounded set of elements of Ai one can extract a subsequence that converges in 
A; we shall write Ai ^-H- A to denote that Ai is compactly embedded in A. 

Suppose that T is a positive real number, ip maps the nonempty closed interval [0, T] into 
a seminormed set M, and p £ M., p > 1. We denote by LP{Q,T; Ai) the set of all functions 
p} : t £[0,T]^ tp{t) £ M such that 

LP(0,T; A^) is then a seminormed set in the ambient linear space LP(0,T;^), with 

We denote by L°°(0,r;A^) and ['p\l'=°(o.t;M) the usual modifications of these definitions when 
p = 00. 

For two normed linear spaces, Ao and Ai^ we shall continue to denote by Aq ^ Ai that Aq is 
(continuously) embedded in Ai. The following theorem is due to Dubinskii [21]; see also [13]. 



(/J 
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Theorem 5.1 (Dubinskii |21)). Suppose that Aq and Ai are Banach spaces, Aq ^ Ai, and A4 

is a seminormed subset of Aq such that A4 ^-H- Aq. Consider the set 



dip 



dt 



< oo 



LPl(0,T;^l) 



where 1 < p < oo, 1 < pi < oo, \\ ■ is the norm of Ai, and dip/dt is understood in the sense 
of Ai-valued distributions on the open interval (0,T). Then, y, with 

dip 

LPi{a,T;Ai) 

IS a seminormed set m LP{0,T;Ao) n W^'P^{0,T;Ai), and y ^ LP{0,T;Aq). 



']y ['P]lp{0,T-M) + 



We note that in Dubinskii [21] the author writes M instead of our M>o in l|5.ip and property 
(ii). The proof of Thm. 1 in Dubinskii's work, stated as Theorem 15.11 above . reveals however that 
the result remains valid with our weaker notion of seminormed set, as (5.1) and property (ii) are 
only ever used in the proof with c > 0; we refer to our paper [13) for details. In the next section, 
we shall apply Dubinskii's theorem by selecting 



and 



Aq — L\,j{^ X D) with norm 



M = <ip e Aq ■■ 'P>0 with 



mlAo - 



Miq) \'p{x,q)\dxdq 



and, for p £ Ai, we define 



Miq) nV.,^pix,q)\ + \Vg^p{x,q) 



dxdq < oo 



da; dq. 



Note that is a seminormed subset of the ambient space Aq. Finally, we put 

Ai := Ar^H^in X D))' := {ip : Mip G (i7"(f7 x D))'}, 

equipped with the norm H^H^i := \\AI(p\\(^H'{nxD))' ^ and take s > 1 + ^{K + l)d. Our choice of 
Ai is motivated by the fact that, thanks to the Sobolev embedding theorem onflxDc K'^^^'' = 
■^(K+i)d^ the final factor on the right-hand side of (|4.42p can be further bounded from above by a 
constant multiple of \\'p\\L^{o,T;H'(nxD))7 with s > 1 + ^{K + l)d. For such s it then follows, again 
from the Sobolev embedding theorem that, for any p £ Aq, 



miAi 



mip,x)\ ^ 

sup - — < sup 

X€H={nxD)\\x\\H={nxD) x<SH={nxD) 



m\Ao\\X\\L--{nxD) 
\\x\\H'^{nxD) 



< KqWpI 



Ao 



where Kq is any positive constant that is greater than or equal to the constant Ks, the norm of 
the continuous linear operator corresponding to the Sobolev embedding {H'^{n x D) ^)H'^~^ {fl x 
D) ^ L°°{nx D), s > 1 + \{K + l)d. Hence, Aq ^ Ai. 

Trivially, M ^ Aq. We shall show that in fact A4 ^-H- Aq. Suppose, to this end, that B is 
an infinite, bounded subset of Ai. We can assume without loss of generality that B is the infinite 
sequence {(pn}n>i C Ai with [(^„]x < Kq for all n > 1, where Kq is a fixed positive constant. 
We define p„ :— and note that (0„ > and /?„ e Ii\i{VL x D) for all n > 1, with 

\pA?R\^(SlxD-) = \Vn\M<KQ Vn > 1. 

Since Hlj{il x D) is compactly embedded in L\.j{n x D) (see Appendix D for a proof of this), we 
deduce that the sequence {pn}n>i has a subsequence {pnk}k>i that is convergent in L\,j(fl x D); 
denote the limit of this subsequence by p; p £ L\.i{V, x D). Then, since a subsequence of the 
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sequence {j5„j.}fc>i also converges to p a.e. on Q x D and each p„j. is nonnegative onQ x D, the 
same is true of p. Now, define (p :— p^, and note that cp E L\f{^l x D). Clearly, 



^n^~0\\L\,{nxD)= / M {pn,+p)\pn^,-p\dxdq 



< 



< 



II Prik + P\\Llj(nxD) II Prik - P IIl|j(OxD) 

(ll Prik WLfjinxD) + II P IIl|j(J2xD)) II Prik ~ P llL|j(Oxr>)- 



As {pnfc}fc>i converges to pin Llj{n x D), and is therefore also a bounded sequence in L\i{n x D), 
it follows from the last inequality that {ipnk}k>i converges to ip in L\j{il x D) — Aq. This implies 
that M. is compactly embedded in hence the triple M. ^-K> Ao ^ Ai satisfies the conditions 
of Theorem 15.11 

Remark 5.1. In fact, there is a deep connection between Ai and the set of functions with finite 
relative entropy on D; this can he seen by noting the logarithmic Sobolev inequality; 

(5.2) M{q)\p{q)?^og-^ dq<- M{q)\V,m\^dq ^p e HI,{D), 

Jd ~ ~ \\P\\lI,(d) ~ f^JD ~ 



with a constant k > 0; the inequality ()5.2p is known to hold whenever M satisfies the Bakry-Emery 
condition.- Hess(— log Af(g)) > k Id on D, asserting the logarithmic concavity of the Maxwellian on 
D, with the last inequality understood in the sense of symmetric Kd x Kd matrices. The inequality 
()5.2p follows from inequality (1.3) in Arnold et al. |4j, with the Maxwellian M extended by to 
the whole o/M^'* to define a probability measure on R^"* supported on D = Di x ■ ■ ■ x Dk- 

The validity of the Bakry-Emery condition for the FENE Maxwellian, for example, is an easy 
consequence of the fact that 

Hess(-logAf(g)) = Hess [jy, (i|gi|')^ 

(5.3) = diag (Hess (t/idlgiP)) , ■ • ■ , Hess (uK{\\qK?))) , 

for all q :— {qj , . ■ . ^qJc)'^ G -^i x • • • x T)k = D, and the following lower bounds (cf. Knezevic 
& sail jSOj . Sec. 2.1) on the d x d Hessian matrices that are the diagonal blocks of the Kd x Kd 
Hessian matrix Hess(— log A/(q)).- 

a (Hess (c/.(i|g,n)) i^ > (1 - \^^\Vbr%^f > ICd', 

for all qi G Di and all £_i G M.'^, i = 1, . . . , K . Hence, 

f Hess(-logM(g)){> 
for all q E D and all ^ G M^'*, yielding 

Hess(-logM(q)) > Id Vq G -D; 

i.e. K = 1. 

More generally, we see from (j5.3p that if qi G -Di i— > [/i(i|gjp) is strongly convex on Di for 

each i = 1, . . . , K , then M satisfies the Bakry-Emery condition on D. 
On writing 'fi{q) := |p(g)p ( > 0) in (|5.2p . we then have that 



<2{q) 2 f 

(5.4) / M{q)^{q)\0g---^ dq<- M(q)V,J^{q) 

~ ~ m\L\,{D) - - " - 



2 

dq, 



for all (p such that ip > on D and G H\j{D). Taking (p = (p/M where ip is a probability 
density function on D, we have that ||^||li (d) = ||</'||li(_d) = 1/ thus, on denoting by pi the Gibbs 
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measure defined by d/i — M{q)dq, the left-hand side of (|5.4p becomes 



^(^|M)-/^(log^)d,, 



referred to as the relative entropy of (p with respect to M. The expression appearing on the right- 
hand side of (j5.4p is 1/(2k) times the Fisher information, of (p: 



Vglog if{q) 



yglog^(g) ^(g)d/i = 4 



dfj., 



where, E is the expectation with respect to the Gibbs measure fi defined above. o 

Lemma 5.1. Suppose that a sequence {ipn}'^=i converges in L^{0,T; L\j{fl x D)) to a function 
(p G L^{Q,T; L\j[Vl x D)), and is bounded in L°° (Q^T] L\^(VL x D)), i.e. there exists Kq > such 
that ||<Pn||L~(o,T;Lij(OxD)) < ^0 for all n>l. Then, if G LP{0,T; L\,j{n x D)) for all p G [l,oo), 
and the sequence {ipn}n>i converges to p in LP{Q,T] L\j{Vl x D)) for all p G [l,oo). 

Proof. Since {^ri}n>i converges in L^(Q,T; L\j{^ x D)), it follows that it is a Cauchy sequence 
in L^{Q,T;L\i{Q. x £))); thus, for any p G [l,oo), there exists no = no(e,p) G N such that for all 
m,n> no{e,p) we have 



II - ^m\\L\, (Ox D) ^t < J^J^^ 



Hence, for all m,n> no(e,p). 



i/p 



< eSS.SUPtg[o^T]ll^« - ^mlliynxi)) 



< e. 



i/p 



W^n - 'Pm\\Ll^(nxD) 



This in turn implies that {(^,i}n>i is a Cauchy sequence in the function space i^(0, T; L\j{fl x D)), 
for each p G [1, oo). 

Since LP{0,T; L\f{Q x D)) is complete, {pn}n>i converges in LP{0,T; L\j{fl x D)) to a limit, 
which we denote by P(p), say. As, by assumption, {^„}„>i converges in L^(0, T; L\,j{il x D)), and 

LP{0,T;Lii{n X i?)) c L\0,T;L\j{n x Z?)) 



for each p G [1, oo), it follows by uniqueness of the limit that ^(p) = ip for all p G [1, oo). Hence, 
also, ip G LP{0, T; x £))) for all p e [1, oo). That completes the proof. □ 



6. Passage to the limit L ^ oo: existence of weak solutions to the FENE chain 

MODEL with CENTRE-OF-MASS DIFFUSION 

The bounds (|4.28p . (I4.42p and (|4.5ip imply the existence of a constant > 0, depending only 
on B(uo, /, V'o) and the constants C* and C««, which in turn depend only on e, v, Cp{Vl), T, \A\, 
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Oq, fc, K , A, f2, d and h, but not on L or At, such that: 

i-T 



ess.suptg[o,T]l!y, 
+ ess.suptg[o,T] 



1 
At 



At,+ 







At,- 



ly.yf^+(s)|pds 



1 



AtL 



JnxD 



JnxD 



At.+ |2 
e.L I 



dqdxds+ [ I M\VqJ^f]^\^ dqAxAs 
Jo JnxD 

2 



(6.1) 
where 



dt 



At 



M 



dt 



At < a, 



(H={QxD))' 



denotes the norm of the dual space V'^ of with a > ^d, a > 1 (cf. the paragraph 
foUowing (|4.44l) ): and || • \\(H={nxD)Y is the norm of the dual space x D))' of H^{i} x D), 

with s > 1 + ^{K + l)d. The bounds in (|6.ip on the time-derivatives follow from ()4.5ip . and from 
(|4.42p using the Sobolev embedding theorem. 

By virtue of (j4.38p . (|4.35p . the definitions (|4.1a[ b'). and with an argument completely analogous 
to (|4.35p on noting p.Sp in the case of the fourth term in (|6.ip . and using (I4.10p . (I3.20p and recalling 
that L > 1 we have (with a possible adjustment of the constant C^,, if necessary,) that 

T 



ess.suptg[o_-r] 



e,L 

At(,±) 



1 

A< 



At,+ 
ie,L 



At,- 
ie.i 



(s)||^ds + ess.suptg[o^ 



UxD 



■As 



MT{i>ff^\t))AqAx 



ML 



T 

JflxD 



V'^i' )'^ d(j d;5 As 



+ 



(6.2) 



10 JQxD 

Jo 



'^sX I 



Aq Ax As 



dt 



At 



M 



JnxD 

2 



dt 



At < C\. 



At,+ |2 



Aq Ax As 



(H^inxD))' 



On noting (|4.31ai b). (|4Ja[ b). (13^20)) and ([331), we also have that 



(6.3) 

and 

(6.4) 



At(,±) 
e,L 



> a.e. onnx D X [0,T] 



M{q)i^f'j^''^\x,q,t)Aq < 1 for a.e. (x,t) G 1) x [0,T]. 



Henceforth, we shall assume that 
(6.5) At = o{L^'^) asL^oo. 

Requiring, for example, that < At < Co / {L logL), L > I, with an arbitrary (but fixed) constant 



Co will suffice to ensure that (16. 5p holds. The sequences {; 



At(,±) 
'e,L 



}l>i and {V'^^*''^^}l>i as weU 
as all sequences of spatial and temporal derivatives of the entries of these two sequences will thus 
be, indirectly, indexed by L alone, although for reasons of consistency with our previous notation 
we shall not introduce new, compressed, notation with At omitted from the superscripts. Instead, 
whenever L — oo in the rest of this section, it will be understood that At tends to according 
to (|6.5I) . We are now almost ready to pass to the limit with i — >• oo. Before doing so, however. 
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we first need to state the definition of the function Vj" that obeys (|3.20p , for a given -00 satisfying 
(I3.5p . We emphasize that up to this point we simply accepted without proof the existence of a 
function obeying (I3.20p for a given 'ipQ. The reason we have been evading to state the precise 
choice of was for the sake of clarity of exposition. The definition of and the verification of 
the properties listed under p.20p rely on mathematical tools that were not in place at the start of 
Section [3] where the notation 'ip'^ was introduced, but were developed later, in the last two sections. 
The details of 'lifting' ipo into a 'smoother' function ip'^ are technical; they are discussed in the 
next subsection. 

A second remark is in order. One might wonder whether one could simply choose as ipo'^ 
indeed, with such a choice all of the properties listed in (|3.5I) would be automatically satisfied, bar 
one: there is no guarantee that [-^o]^^^ S Hlj{flx D). Although the property [-0"]^/^ e H\j{Q.xD) 
has not yet been used, it will play a crucial role in our passage to the limit with L — >■ oo in Section 
16.21 In fact, in the light of the logarithmic Sobolev inequality (15. 4p . on comparing the requirements 
on 00 in (|3.5|) with those on ■0° in (|3.20p . one can clearly see that the role of the condition 
|-^oji/2 g HIj{Q X D) in p.20p is to 'lift' the initial datum 0o with finite relative entropy into a 
'smoother' initial datum ip'^ that also has finite Fisher information, in analogy with the process 
of 'lifting' the initial velocity uq from H into u° in V. That the choice of 'ipo as ip^ is not a good 
one can be seen by noting the mismatch between the third term in ()6.2p arising from the Navier- 
Stokes equation on the one hand, and the sixth and seventh term in ()6.2p that stem from the 
Fokker-Planck equation. The absence of bounds at this stage on and tp'^* in those terms 

is entirely due to the fact that, to derive (|6.2I) . we did not use that [-0'^]^/^ e ^lii^ ^ D). This 
shortcoming of (|6.2p will be rectified as soon as we have defined -0" and shown that it possesses 
all of the properties listed in p.20p . 



6.1. The definition of V-""- Given -00 satisfying the conditions in (|3.5I) and A > 1, we consider 



the following discrete-in-time problem in weak form: find C"^'^ G H\,i{^l x D) such that 



^A.l _ /-A.O 

(6.6) / M- ^ — (pdqAx+ I Af V^C^^^ • V^^-l- VgC^^i • Vq^ dgdx = 



n-xD ~ ~ JnxD 

for all if e Hliin X D), with C^'O :== /3'^(0o) e Llf{n x D). Here fi^ is defined by (frT2l) . with 
L replaced by A. The function J^^, which we shall encounter below, is defined by (j4.7p . with L 
replaced by A. 

The existence of a unique solution (^'^ g H\,j{D. x D) to (|6.6p . for each Ai > and A > 1, 
follows immediately by applying the Lax-Milgram theorem. The parameter A plays an analogous 
role to the cut-off parameter L; however since we shall let A ^ oo in this subsection while, for 
the moment at least, the parameter L is kept fixed, we had to use a symbol other than L in 
(|6.6p in order to avoid confusion; we chose the letter A for this purpose in order to emphasize the 
connection with L. 

Lemma 6.1. Let C,^'^ he defined by (j6.6p . and consider j^'^ defined by 
(6.7) 7^'"(S) := / M(9)C^'"(s,9)d«, n = 0,l. 



JD 

Then, (^^'^ is nonnegative a.e. on nx D, andO< 7^^^ < 1 a.e. on fl. 

Proof. The proof of nonnegativity of C^'^ is straightforward (cf. the discussion following (|3.66p ). 
Indeed, we have that = a.e. on n X D, thanks to p.Sp and the definition of /3 ; 

we then take ip — [C^'^]- as a test function in (|6.6I) . noting that this is a legitimate choice 
since C^'^ € H\j{Vl x D) and therefore [C"^^^]- e H\.i{Vl x D) also (cf. Lemma 3.3 in Barrett, 
Schwab & Siih i8j). On decomposing C,^-^ = \C,^'^]+ + [C^'^]-, and using that [C^^i]+ [C^^i]_ = 0, 
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Yx[C^'^]+ ■ YxiC'^'^]- = and YqiC^'^h ■ YqK^'^]' = a.e. on 17 x D, we deduce that 

= j-J MC^'°[C^'Vdgds = i- / M[C^-0]+[C^^Vdgdx<0, 

«/ O X D "/Ox D 

where || • || denotes the L^{n x D) norm. This then imphes that \\M^ [C^'^]-\\'^ ^ 0. Hence, 
[C'^'^]- — a.e. on X 13. In other words, ^^'^ > a.e. on x 13, as claimed. 

In order to prove the upper bound in the statement of the lemma, we proceed as follows. 
With 7^'"' as defined in (|6.7p . we deduce from the definition of C^'^ and Fubini's theorem that 
£ H^{il). Furthermore, on selecting (p ~ ip <^ H^{n) ® 1{D) in (|6.6I) . recall p.63p . we have 

that 

(6.8) / j^,o ^^^^ r ^^^^^^ ^^^^^^^ y^eH^n). 

Jn ~ Jn~ ~ ~ 

As C^^° = P^{-ipo), and < /3^(s) < s for all s £ R>o, we also have by (|33]) that 

(6.9) 0<7'^'°= / MP^{$o)dq< I M$odq = l on n. 

Jd ~ J D 

Consider z^-" 1 - 7^-", n = 0, 1. On substituting 7^'" ^ 1 - .z'^^", n = 0, 1, into (HH), we 
have that 

(6.10) / ipAx+ / V^z^^^ ■ V:,ipdx^Q y^cH^{VL). 

Also, by ()6.9p . we have that < z^''^ < 1. By using an identical procedure to the one in the first 
part of the proof, we then deduce that — a.e. on n. Thus, z^'i > a.e. on Q, which 

then implies that 7'^'^ < 1 a.e. on fl, as claimed. □ 

Next, we shall pass to the limit A — > 00; as we shall see in the final part of Lemma l6.2l below. this 
will require the use of smoother test functions in problem (|6.6p . as otherwise the term involving 
^A,o _ ^Aj-^p-j jg j^Q^ defined in the limit. In any case, our objective is to use the limit of the 
sequence {C^'^}a>i, once it has been shown to exist, as our definition of the function ■0°. We shall 
then show that thus defined has all the properties listed in p. 201) . 

To this end, we need to derive A-independent bounds on norms of C,^'^, very similar to the L- 
independent bounds discussed in Sectional Since the argument is almost identical to (but simpler 
than) the one there (viz. (|6.6p can be viewed as a special case of (j3.21b[) . with /", l 
taken to be identically zero, A = i,e=l,iV = l, and A chosen as the K x K identity matrix), 

we shall not include the details here. It suffices to say that, on testing (|6.6I) with F'{C,^'^ + a) and 
passing to the limit a — >■ 0+, analogously as in the proof of ()4.27p in Section |4j we obtain that 

[ MF{Q^'^)Aqdx+ AM [ A/jy^rVC^I^dqdx 
JnxD ~ ~ JnxD ~ ~ ~ 

(6.11) +AAt [ M dqdx < [ MT{ipo)dqdx. 

JnxD ~ JnxD 

Our passage to the limit A — >■ 00 in (j6.6l) is based on a weak compactness argument, using (|6.1ip . 
and is discussed below. 

We have from Lemma l6.ll that {[C^'^]^}a>i is a bounded sequence in L\f{D, x D). Using 
this in conjunction with the second and third bound in (I6.1ip we deduce that, for At > fixed, 
{[C^'^]^}a>i is a bounded sequence in H\.^{Vl x D). Thanks to the compact embedding of H\j{Vl x 
D) into L\,i{Vt X D) (cf. Appendix D), we deduce that {[C'^'^]^}a>i has a strongly convergent 
subsequence in L\.i{Q, x _D), whose limit we label by Z, and we then let C,^ := 1? . For future 
reference we note that, upon extraction of a subsequence (not indicated), C,^'^ then converges to 
C, ^ a.e. on n x D; and C^'^(g;, •) converges to C ^(Sj ') a.e. on D, for a.e. a; G fi. 
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By definition, we have that C ^ ^ 0; furthermore, thanks to the upper bound on 7^'^ stated in 
Lemma 16.11 the remark in the last sentence of the previous paragraph, and Fatou's lemma, we 
also have that 

(6.12) / M{q)C^{x,q)dq <l for a.e. x e fl. 

Jd ~ ~ ~ ~ 

Further, again as a direct consequence of the definition of ^, we have that 



(6.13) y C^^.i ^ V ^ strongly in Llj{fl x D). 

Application of the factorization ci — C2 = (Vci ^ V^J) (\/^+ V^J) with ci, C2 G R>0; the Cauchy- 
Schwarz inequality and (j6.13l) then yields that 

(6.14) C^'i ^ strongly in Llj{n x D). 
Finally, we define 

(6.15) V'°:=C^- 
It follows from the nonnegativity of C ^ and (|6.12p that 

(6.16) <>0 a.e.onQxD and 0< / M{q) ip°{x,q) dq < I ior a.e. x £ Q. 

Jd ~ ~ ~ ~ 

Further, from the bound on the first term in (|6.11l) and Fatou's lemma, together with the fact 
that, thanks to the continuity of J^, (a subsequence, not indicated, of) {^(C'^'^)}a>o converges to 
J^(C ^) = a.e. on X D, we also have that 



(6.17) / M F{ip'^)Aqdx< A/ J"(V'o)dgdx. 
Next, we note that from (|6.13p we have that, as A — > 00, 

(6.18) Af5 ^ strongly in L'^{n x D). 

We shall use (I6.18P to deduce weak convergence of the sequences of x and q gradients of C,^'^ . 
We proceed as in the proof of Lemma 13.31 The bound on the third term on the left-hand side of 
(|6.1ip implies the existence of a subsequence (not indicated) and an element g E L^{il x D), such 
that 



(6.19) M^Vq\/c^'^->-g weakly in L2(f| X D). 

Proceeding as in (|3.60p - (|3.61l) in the proof of Lemma [3.31 with tp" g, tp" l S ^ Q+ replaced 
by v C^'^, y and A — > 00, respectively, we obtain the weak convergence result: 



.20a) Mi Vg\/C^^i -> Aft Vg\/C 1 weakly in L^{n x D), 

and similarly for the a; gradient 



(6.20b) Af5 V^yC^^'i Af3 yC 1 weakly in L^{n x D), 

as A — > 00. Then, inequality (|6.1ip . (|6.20al b) and the weak lower-semicontinuity of the L^(f2 x D) 
norm imply that 



(6.21) AAt f M + dqdx< f MT{ipQ)dqdx. 

After these preparations, we are now ready to state the central result of this subsection. Before 
we do so, a comment is in order. Strictly speaking, we should have written instead of in 
our definition (|6.15p . as "0" depends on the choice of At. For notational simplicity, we prefer the 
more compact notation, with the dependence of on At implicitly understood; we shall only 
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write V'^j, when it is necessary to emphasize the dependence of At. Of course, i/jq is independent 
of At. 

We shaU show that, with our definition of the properties under (|3.20p hold, together with 
some additional properties that we extract from (|6.15p . 

Lemma 6.2. The function tp^ — ip'^^ defined by (|6.15p has the following properties: 
O < e Zi; 



& 



JnxD ~ ~ JnxD 



inxD 
© 4At 



M 



flxD 



0|2 



dq dx < 



nxD ' ^ 

limAt-i-o+ V'" = "00, weakly in L\,j(fl x D); 

© limAt-i.o+ = "00; weakly in L\,j{n x D). 



M J'(V'o)d(? d^; 



nxD 



Proof. 



O This is an immediate consequence of (|6.16p and the definition p.l9p of Zi. 
@ This property was established in (|6.17p above. 

© The inequality follows by using (|6.15l) in the left-hand side of (|6.21l) . 
We begin by noting that an argument, completely analogous to (but simpler than) the 
one in Section l4.2.1l that resulted in (I4.42p . applied to (|6.6p now, yields 



^A,l _ ^A,0 

M — (p dq dx 



nxD 



At 



< 2 



M 



nxD 



dq dx 



L°°(_D) + llYg'/'llL~(D) 



da; 



for ah H^{n;L°°{D))nL^{n;W^'°°{D)). By noting (lOTl) we deduce that 



M{C^'^ -C''")vdq dx 



A,0\ 



nxD 



nxD 



M T{^Jo) dq dx 



(6.22) 



|V:r^|li=o(c) + llVq^iiio 



dx 



for all if e H^{n;L°°{D)) n L^{n;W^'°°{D)). As the right-hand side of is inde- 

pendent of A, we can pass to the limit A — > oo on both sides of (I6.22p . using the strong 
convergence of C^'^ to -0" in Llj{n X D) as A -> cx) (see (|5T^ and the definition of ^J^) 
together with the strong convergence of (^'^ — /^^('i/'o) to "00 in ^m(^ x Z?) as A — > oo, 
with At kept fixed. We deduce that 



M (-0° --0o)^dq da; 



QxD 



MF{^o)dq dx 



(IxD 



(6.23) 



2 

L~(I5) 



L'^{D) 



dx 



for all (p G H'^{n;L°°{D)) n L'^{VL]W'^^°°{D)) and therefore in particular for all (p € 
H^{n X D) with s>l + \{K + l)d. 

As the last two factors on the right-hand side of (I6.23P are independent of At, we can 
pass to the limit At 0+ on both sides of (|6.23p to deduce that -0° = ip^^ converges to 
V^o weakly in M-^{H''{VL x D))' , s>l + \{K + l)d, as At ^ 0+. 

Noting (|6.17p and the fact that F{r)/r — )■ cx) as r — cxd, we deduce from de la Vallee- 
Poussin's theorem that the family {V'At}At>o is uniformly integrable in L\i {Q.xD). Hence, 
by the Dunford-Pettis theorem, the family At>Q is weakly relatively compact in 

L\,j{Vl X D). Consequently, one can extract a subsequence {V'Atfcife^i that converges 
weakly in L\j{Q. x D)\ however the uniqueness of the weak limit together with the weak 
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.24) 



convergence of the (entire) sequence -^^ = ^/j^^ to -i/jq in M ^{H'^{n x D))' , s > 1 + ^{K + 
l)d, as At 0+, established in the previous paragraph, then impUes that the (entire) 
sequence ip'^ — -0^^ converges to i/'o weakly in L]^{U x I?), as Ai — >■ 0+, on noting that 
Llj{n X D) is (continuously) embedded in M'^H^n x D))' for s > 1 + i(i<: + l)d (cf. 
the discussion following Theorem IS.ip . 
It follows from -0" € Zi and (fri2|) that 







< / MLdqdx< [ A//3^(i/^°)dgdx < / M$°dqdx<\n 
J^,«>L ~ ~ JnxD ~ ~ JnxD ~ ~ 



.25) 
.26) 



On noting that J- is nonnegative and monotonically increasing on [1, cx)), and that J-{s) € 
[0, 1] for s e [0, 1], we deduce that 

MT{[ip° - L]+)dq dx 

MJ-([0"-L]+)dqdx+ / MT{[il^" - L] + )dqdx 

< [ Mdqdx+ [ MF{^°)dqdx<C. 
JnxD ~ ~ JnxD ~ ~ 

Let us recall the logarithmic Young's inequality 

r s < r log r — r + e* for all r, s e M>o. 
This follows from the Fenchel- Young inequality: 

r s < g*{r) + g{s) for all r, s G R, 

involving the convex function g : s G M i— > g(.s) G (— oo, +oo] and its convex conjugate g* , 
with g{s) = e" and 

{+c» if r < 0, 

if r = 0, 

r (log r — 1) if r > 0; 



.27) 



with the resulting inequality then restricted to R>o- It immediately follows from (|6.26p 
that rs< T{r) + e" for ah r, s G M>o. 

Applying the last inequality with r = [-0" — L]+ and s — logL, we have that 

[$°-L]+{\ogL)<J^{[$'>~L] + )+L. 



The bounds (|6.24p , (|6.25p (noting that the integrand of the left- most integral in (|6.25p is 
nonnegative) and (I6.27P then imply 



/ M[$° - L]+dqdx = I M[i)^ ~ L\+dqdx 

JnxD ~ ~ J$o>L ~ ~ 



(6.28) 



< 



logL 



ipo>L 

MT{[ip° - L] + )dqdx+ / MLdqdx 

'4i'>>L ~ ~ J^°>L ~ ~. 



< 



C 
logL' 



Hence, for any ip G L°°(fl x D) we have from (|6.28l) on recalling the relationship 
At = o(L^i) that $° = V'At satisfies 



lim 

At->-0+ 



M« - /3^(Vr))^dg dx 



nxD 



lim 



M[0°-L]+(^d(7 dx 



nxD 



.29) 



< ( lim / M[0O-L]+d<zda; ||^|U»(oxD) 



Therefore, similarly to (j6.23p . we have that the sequence {'ip^ — /3^(V'At)}At>o converges 
to zero weakly in M'^{H''[n x D))' for s > \{K + l)d, as At 0+. 
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Noting (|6.25l) and the fact that J^(r)/r oo as r — ?► oo, we deduce from de le Vallee 
Poussin's theorem that the family 

- /?^(^L)}At>0 ^ Wit - L] + }At>0 

is uniformly integrable in L\j{il x D). Hence, we can proceed as for the sequence 
{ipAt}'^t>o in the proof of O to show that the (entire) sequence 

V^o - = ifAt - P^i$%t) ^ weakly in L\j{n x L>), as At ^ 0+, 

on noting that L\,i{Q, x D) is (continuously) embedded in M~^{H^{Q, x D))' for s > 
^{K + l)d (cf. the discussion following Theorem 15. ip . Hence, we have proved the desired 
result. 

□ 



Noting item 3 in Lemma 16.21 we can now return to the inequality (j6.2[) , and supplement it with 
additional bounds, in the sixth and seventh term on the left-hand side. The first additional bound 
can be seen as the analogue of (|4.35p : 



M 



(6.30) 



JnxD 
= AAt 



M 

T-At 



< 4At 





M 

T-At 



< 



+ 4 



flxD 



+ 4 



dq dx dt 

iy.V/?^(v^")i' + iy,\//3^(v>)i' 

M 

dq dx 

M 



dq dx 



nxD 



iy.V'A°P + iy,V'/'°i' 



nxD 
M T{$o) dq dx 

M 



dq dx dt 



dq dx dt 



JnxD 



dq dxdt < C^, 



where in the last inequality we used p.20p and the bounds on the sixth and seventh term in 
()6.2p : here and henceforth C* signifies a generic positive constant, independent of L and At. On 
combining ()6.30p with our previous bounds on the sixth and seventh term in (j6.2p . we deduce that 



.31) 



M 



JnxD 



dq dxdt < C^. 



It remains to derive an analogous bound on "0^1^. To this end, let n £ {!,... ,N} and consider 
t e {tn-i,tn)', we recall that 



(6.32) 



At 



For ease of exposition we shall write 

t — tn~l 



7+ 



At 



and 



7- 



At 



in the argument that follows, noting that 7+ + 7_ = 1 and both 7+ and 7^- are positive. The 



functions t S (tn-ij^n) 4'^L^i',','t) are constant in time and 'ijj^^£'"^'{x,q,t) > over the set 



■At,(±) 



44 



JOHN W. BARRETT AND ENDRE SULI 



il X D X {tn-i,tn), n £ {1, . . . , N} . For any a £ (0, 1) we have that 



At 12 

e,L\ 



7+(^^l'+ + a)+7-(^^l'-+a) 



At -|2 



At 



a 



< 2 



< 2 



Hence, on boundmg 7* by 1, we deduce that 
(6.33) 



At 



for all {x,q,t) £ ilx D X (t„-i, t„), n = 1, . . . , iV, and ah a e (0, 1). On multiplying (|633)) by M, 
integrating over n x D x (tn-i, tn), summing over n — 1, . . . , N , and passing to the limit a — > 0+ 
using the monotone convergence theorem, we deduce that 



At,+ 



At- 



JnxD 



Af y^i/Vfil dgdxdi < 2 



JnxD 

T 



JnxD 



Analogously, 



JnxD 



M\Yq\^^tx\ dqdxdt ^ 2 



JnxD 

T 



M\VgJiif*£+\^dq dxdt 



M\Yq\/TptL'\'^dq dxdt 



JnxD 

Summing the last two inequalities and recalling (j6.3ip . we then deduce that 

i-T 



(6.34) 



M 



^At 1 2 



■■e.Ll 



'At 1 2 



dq dx dt < Ci,, 



/o JnxD 

where, again, denotes a generic positive constant independent of L and At. 

Finally, on combining (|6.31l) and (|6.34p with (|6.2I) we arrive at the following bound, which 
represents the starting point for the convergence analysis that will be developed in the next 
subsection. 

With a > ^d, (T > 1 and s > 1 + ^{K + l)d, we have that: 

1 rT 

ess.supjgfo^rjiiu^^^' " 



At,+ , 



\u-z'-{t)^uf'j:-{t)rdt 



Yx ufi: {t)fdt + ess.supt 



e[o,T] 



nxD 



M T{i^ff^\t))dq dx 



AtL Jq JnxD ^'^ 



M{$^j'+ -$^}-)'dqdxdt 



JnxD 



^At(,±)|2 



dq dx dt 



JnxD 



rAf(,±)|2 



dq dx dt 



(6.35) 



dt 



it) 



dt^ 



M 



dt 



it) 



dt < a. 



(H={nxD))' 
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Similarly, 



(6.36) 



At 



dt 



it) 



dt < C*. 



Here, the bound on the first term on the left-hand side follows from (|4.6I) . (|6.3p and (|6.4p : the 
bound on the second term comes from (|4.29p . and the bound on the last term from (|4.43p . The 
bound on the third term on the left-hand side of (|6.36p is obtained by applying Vx to both sides 
of (|4.6p with the integrand 



rewritten as 



Miq) 



n 2 



then exchanging the order of Vj^ and the integral over D on the right-hand side of the resulting 
identity, applying V^; to the integrand using the chain rule, followed by taking the modulus on both 
sides and applying the Cauchy-Schwarz inequality to the integral over D on the right, integrating 
the square of the resulting inequality over [0, T] x $7 and, finally, recalling again the definition 
()4.6|) and using the bound on the first term in ()6.36p and the bound on the sixth term in ()6.35p . 
In fact, in the case of p^^l^ the stated bound on the third term on the left-hand side of (|6.36p 
follows directly from (|4.29p . 

6.2. Passage to the limit L — > oo. We are now ready to prove the central result of the paper. 

Theorem 6.1. Suppose that the assumptions p.Sp and the condition ()6.5p . relating At to L, 
hold. Then, there exists a subsequence o/{(y^i,'0^i)}L>i (not indicated) with At — o{L^^), and 
a pair oj Junctions {ue,4'e) such that 

Ue e L°"iO,T; L^{n)) n L'^iO,T;V) n H\0,T;V'^), a > id, a > 1, 

and 

e L\0,T;Llj{n X D)) n H\0,T; M-\H%n x £>))'), s>l + ^{K + l)d, 
with > a.e. on n x D x [0,T], 



(6.37) pe{x,t):= M{q)'iljAx,q,t)dq = l for a.e. {x,t) £ n x [0,T], 

Jd ~ ~ ~ ~ 

whereby ip^ G L°°{0,T; L\,j(fl x D)); and finite relative entropy and Fisher information, with 

(6.38) L^{'^,T;L\,i{Vtx D)) and \f^, £ L'^{0,T; Hl^fl x D)), 
such that, as L oo (and thereby At — > 0+ ), 



(6.39a) 


At(,±) 


— ?> Ue 


weak* in L°°{0,T; L^{n)), 


(6.39b) 


At(,±) 


— )• Ue 


weakly in L^(0,T;V), 


(6.39c) 


At{,±) 


-> Me 


strongly in L^{Q,T-U'{n)), 


(6.39d) 


du^X 
dt 


dUe 


weakly in L^(0, T; V'^), 
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where r e [1, oo) if d ~ 2 and r £ [1,6) if d — 3; and 
(6.40a) 
(6.40b) 

(6.40c) 
(6.40d) 



Mi "^x\Ji^f^t^^ Miv^y^)^ weakly in L^{0,T;L^{n x D)), 
Vq\f$^^ Vg^/V-e weakly m L^{0,T;L^{n x £>)), 



M 



Q^At 



dt 



dt 



7At(,±) ^ 7 



t/;ea% in L^{0,T; {H^in x £>))'), 
strongly in LP{0,T; Llj{n x D)), 



for all p [1, oo); and, 



K 



K 



(6.40e) V, • ^C,(AfV^^['+) ^ V, • ^C,(M^e) wea^ L2(0,T;y;) 

i=i ~ ~ i=i ~ 

The pair {u^^tpi,) is a global weak solution to problem (Pe), in the sense that 



Jfi 



dw 

Ue ■ dx dt - 
~ ot ~ 



T 

Jn 



K 



dx dt 



(6.41) 

and 



uq{x) ■ w{x,Q)dx + / (/, w)jji(o) - fc / CiiMxj)^) -."sj ^wdx 

~ ~ ~ ~ ~ Jo \_ ~ ~ i=i-^^~ 



dt 



JilxD 
T 



^ dip ''^ 
M ip^ — — dq dx dt 
ot ^ ~ 



M 



JnxD 



jx [ I ^ ^ ■ ^ '^'^ "^f 



Va;<^ d<7 dx dt 



K K 



i=i j=i 



V'e ■ dq dx dt 



'ijjo{x,q)ip{x,q,0) dq dx 



nxD 



(6.42) 



V(^e W'^'\0,T;H''{n X D)) s.i. (^(^^r) ^ 0. 



/n addition, the function is weakly continuous as a mapping from [0, T] to H , and ip^ is weakly 
continuous as a mapping from [0,T] to L\i{Q. x D). The weak solution {ue,ipe) satisfies the 
following energy inequality for a.e. t € [0,T]: 

\\ue{t)f + iy[ \\YxUeis)fds+ k [ MT{i^e{t))dqdx 
Jo ~ ~ JnxD 



+ 4ke 



JnxD 
t 



M\^xy '(l^e\ dq dx ds 



a^k 



JnxD 



M\^q\^tpe\ dqdxds 



(6.43) 



< 



\Hr + - ll/(5)||fffi(o))'ds + fc/ MF{i:o)dqdx<[E{uo,f,^o)] 

>^Jr, ~ \ ow; JnxD 



with J-{s) ~ s(log.s — 1) + 1, s > 0, and [B(uo, /, -^o)]^ as defined in (|4.28l) . 

Proof. Since the proof is long, we have broken it up into a number of steps. 

Step 1. On recaUing the weak* compactness of bounded bahs in the Banach space L°°{Q, T; L^(r2)) 
and noting the bound on the first term on the left-hand side of (I6.35p . upon three successive ex- 
tractions of subsequences we deduce the existence of an unbounded index set C C (1, oo) such that 
each of the three sequences {u^^^'^^jiec converges to its respective weak* limit in L°°{0, T; i^(n)) 
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as i -> OO with L e C. Thanks to (|4.1al b). 

(6.44) r WufXi^) - v!^%^{s)fAs = \ r \\uf^l+{s) - y,^i'-(=s)|p d. < ^C.At, 
Jo Jo 

where the last inequality is a consequence of the second bound in (j6.35p . On passing to the 
limit L — ^ cx) with L ^ C and using (I6.5|) we thus deduce that the weak* limits of the se- 
quences {u^*^'^'^}Lec coincide. We label this common limit by u^; by construction then, e 
L°°(0,r; L'^[n)). Thus we have shown (|6.39ap . 

Upon further successive extraction of subsequences from {y^i '^''}Le£ and noting the bounds 
on the third and eighth term on the left-hand side of (|6.35l) the limits (j6.39bl d') follow directly 
from the weak compactness of bounded balls in the Hilbert spaces L^(0,T; V) and L'^{0,T; V'^) 
and (I6.39ap thanks to the uniqueness of limits of sequences in the weak topology of L^(0,T; V) 
and L^(0, T; Fjj), respectively. 

By the Aubin-Lions-Simon compactness theorem (cf. p.l4[) ). we then deduce (|6.39c|) in the 
case of on noting the compact embedding of V into L^{^) H H, with the values of r as in 
the statement of the theorem. In particular, with r = 2, {t/^^}ig£ converges to u^, strongly 
in L^(0, T; i^(57)) as L — > oo. Then, by the bound on the left-most term in (|6.44p . we deduce 
that {uf]^^}Lec also converges to u^, strongly in L^{0,T; L^{V,)) as L — )■ oo (and thereby At — > 
0+). Further, by the bound on the middle term in (|6.44p we have that the same is true of 
{u'^*£ }l<£C- Thus we have shown that the three sequences {u'^^^'^^}li£C all converge to Ue, 
strongly in L^(0, T; L^(^l)). Since the sequences {yfi'^^jiGz: are bounded in L^(0, T; H^{ft)) (cf. 
the bound on the third term in (j6.35p ) and strongly convergent in L'^{0, T; L'^{n)), we deduce from 
(13. 4p that (|6.39cp holds, with the values of r as in the statement of the theorem. Thus we have 
proved (j6.39al -d). 

Step 2. Dubinskii's theorem, with ^o, Ai and A4 as in the discussion following the statement 
of Theorem 15. II and selecting p = 1 and pi = 2, implies that 

(f : [0,T]^ M : [ip]mo,T;M) + 

^ L\0,T;Ao) = L\0,T;Llj{n x D)). 

Using this compact embedding, together with the bounds on the sixth, the seventh and the last 
term on the left-hand side of (|6.35l) , in conjunction with (16. 3p and (j6.4p , we deduce (upon extraction 
of a subsequence) strong convergence of {V'^^}l>i in L^iO,T; L\.i{n x D)) to an element ip^ € 
L\0,T;Llj{nx D)), as oo. 

Thanks to the bound on the fifth term in (j6.35p . by the Cauchy-Schwarz inequality and an 
argument identical to the one in (|6.44l) . we have that 

/7 MlV^^^i-V^fridgdxdA <^/7 M(^fr-V?r)'d<zd:.dt 

Jo JilxD ~ ~ J Jo JflxD ~ ~ 

(6.45) <lC^T\n\AtL. 

On recalling (|6.5p . and using the triangle inequality in the L^{0,T; L\j{^ x D)) norm, together 
with (|6.45p and the strong convergence of {V'^|^}l>i to in L^{0,T; L\,j{Q x D)), we deduce, 
as L oo, strong convergence of {ipfi,~^}L>i in L^{0,T; L\j{n x D)) to the same element e 
L^(0, T; L\.f{n X D)). The inequality (|6.45p then implies strong convergence of {V'^^'~}l>i to ip^, 
also. This completes the proof of (j6.40dp for p = 1. 
From (|6.3p and (|6.4p we have that 

(6.46) IIV^^l''^^(t)llLi,(nxi.) < M 



dip 
'dt 



< oo 



L2(0,T;^l) 
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for a.e. t in [0,T] and all L > 1. In other words, the sequences {V'^^*''^''}l>i are bounded in 
L°°{0,T;Lli{nxD)). By LcmmaEHl the strong convergence of these to -^e in (0, T; L]^^ (f7 x Z))) , 
shown above, then implies strong convergence in Lp{Q,T; L\.j{n x D)) to the same limit for all 
values of p S [l,oo). That now completes the proof of (|6.40d|) . 

Since strong convergence in £''(0, T; L\j{il x -D)), p > 1, implies convergence almost everywhere 
onQxDx[0, T] of a subsequence, it follows from (I6.3P that > on fix D x[0,T]. Furthermore, 
by Fubini's theorem, strong convergence of {V'f^L '*''}l>i to ip^ in L^{0,T; L\,i{fl x D)) implies 
that 

/ M{q)\%l)f^^'^\x,q,t) — il}^{x,q,t)\(lq ^ Q as i ^ oo 

J D ~ ' ~ ~ ~ ~ 

for a.e. {x,t) e n x [0,T]. 

Hence we have that Jj^ M{q) ^pf*j^''^\x,q,t) dq converges to M{q) -tjj^lx, q,t) dq, as L — ^ oo, 
for a.e. {x,t) e fl x [0,T], and then ()6.4p implies that 

(6.47) / M{q)^^^{x,q,t)dq <1 for a.e. (a;, i) G x [0, T]. 
Jd ~ ~ ~ 

We will show later that the inequality here can in fact be sharpened to an equality. 

As the sequences {V'^l '*''}l>i converge to ipe strongly in L^(0, T; L\j{il x D)), it follows (upon 
extraction of suitable subsequences) that they converge to tpe a.e. on $7 x x [0,T]. This then, 
in turn, implies that the sequences {-^('0^^ ''''^)}l>i converge to T{ipe) a.e. on x _D x [0,T]; in 
particular, for a.e. t € [0,T], the sequences {T{ipf*j^''^\-,-,t))}Lyi converge to J"(V'£(-, ■, i)) a.e. 
on n X D. Since T is nonnegative, Fatou's lemma then implies that, for a.e. t E [0, T], 

[ M{q)F{^e{x,q,t))dxdq 
JnxD ~ ~ ~ 

(6.48) < lim infi^^oo / M{q) J'{i^fl^-^\x, q,t)) dxdq < C,, 

JnxD ~ ' ~ 

where the second inequality in (|6.48p stems from the bound on the fourth term on the left-hand 
side of (|6.35p . As the integrand in the expression on the left-hand side of (|6.48p is nonnegative, we 
deduce that J-{^e) belongs to L°°(0, T; L\.f{fl x D)), as asserted in the statement of the theorem. 

We observe in passing that since | ^/c2 \ < \/|ci — C2I for any two nonnegative real numbers 
ci and C2, the strong convergence (I6.40d[) directly implies that, as L 00 (and thereby At 0+), 

(6.49) y^JF^^y^ strongly in LP{0,T;Llj{n X D)) Vpe[l,oo), 
and therefore, as L 00 (and At — > 0+), 

(6.50) M5 y/^^^p^ ^ Af ^ strongly in iP(0,r;i2(f7 X £>)) Vpe[l,oo). 

By proceeding in exactly the same way as in the previous subsection, between equations (j6.18p 
and (j6.20bp . with C^'^ and C ^ = t/j" replaced by V'^^*''^'' and ipe, respectively, but now using the 
sixth and the seventh bound in (|6.35p . and (|6.4p . we deduce that ()6.40al b) hold. 

The convergence result (I6.40cp follows from the bound on the last term on the left-hand side 
of (|6.35l) and the weak compactness of bounded balls in the Hilbert space L^{0, T; {H''{il x D))'), 
s>l + ^{K + l )d. 

The proof of (|6.40el) is considerably more complicated, and will be given below. 

After all these technical preparations we are now ready to return to (|4.3p and (14. 4p and pass to 
the limit L — > 00 (and thereby also At — > 0+); we shall also prove (|6.40ep . Since there are quite 
a few terms to deal with, we shall discuss them one at a time, starting with equation ()4.4p . and 
followed by equation (|4.3p . 

Step 3. We begin by passing to the limit L 00 (and At — > 0+) on equation (|4.4p . In what 
follows, we shah take test functions G C^{[0,T];C'^{n x D)) such that ^(^sT) = 0. Note 
that, for any s > 0, the set of all such test functions ^ is a dense linear subspace of the linear 
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space of functions in W^'^{0,T; H'^{fl x D)) vanishing at t ~ T. As each of the terms in (|4.4I) 
has been shown to be a continuous hnear functional with respect to (p on L^(0, T; H^{^l x D)) for 
s > 1 + i(K + l)d, and therefore also on W^^^iO^T; H'^in x D)) for s > 1 + ^{K + l)d, which 
is (continuously) embedded in L^(0, T; iJ*(r2 x D)) for s > 1 + ^{K + l)d, the use of such test 
functions for the purposes of the argument below is fully justified. 

Step 3.1. Integration by parts with respect to t in the first term in (|4.4|) gives 

r I ^ dxdt^~ r f M^jfl ^ Aq dxAt 

Jo JnxD cit ^ ~ Jq Jqxd at ~ 

(6.51) - / M{q)(3^{^\x,q))^{x,q,Q)dqdx 

JnxD ~ ~ ~ ~ ~ ~ ~ 



for all $ e C\[0,T];C°°{n x D)) such that (p{-, ■,T) ^ 0. Using (|6.40dp and noting point © of 
Lemma 16.21 we immediately have that, as L — >■ oo (and At 0+), the first term on the right-hand 
side of (|6.5ip converges to the first term on the left-hand side of (|6.42p and the second term on the 
right-hand side of ()6.5ip converges to — j^y^jj ij^oix, q) lp{x, q, 0) dq dx, resulting in the first term 
on the right-hand side of (j6.42p . That completes Step 3.1. 

Step 3.2. The second term in (|4.4[) will be dealt with by decomposing it into two further terms, 
the first of which tends to 0, while the second converges to the expected limiting value. We proceed 
as follows: 



T 

JnxD 



M V^?^^['+ • V^^ dq dx dt 



2e / / M f v/^fl,'+ - J^e\ V,7^f['+ • V,^ dqdxdt 
Jo JnxD \ y~V' ^ ^ ^ 

+ 2e [ [ M \/^,V^\f$f^ ■ V^ip dqdxdt 



Jo JnxD ■ ~ • ~ ^ ~ 

=: Vi +V2. 

We shall show that Vi converges to and that V2 converges to the expected limit. First, note 
that 

T 

"'SI 



|Vi|<2e / / ( / MK/V^,^r- V^ePdg 



y^A/lV^y^g+pd^y ||V,^|Uo.(^) dxdt 

<2e/ ([ M|7^g+- V^pdqda; 

\JnxD ~ ~ 

f A/lV^J^^i'+pdgdxV ||V,^|Uoo(nxD)di 

JnxD ~ V . ^ ^ 



X \\M2 "^x\l'$fj^^\\L'^{nxD) ||Va;^||L~(J7x_D) dt 

I yM^g+ - ^M^eVmnxD) A \ fwY-M S(nxD) 
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were r e (2, oo). Using the bound on the sixth term in (j6.2p together with the Sobolev embedding 
theorem, we then have (with C* now denoting a possibly different constant than in (|6.2[) . but one 
that is still independent of L and Ai) that 



|Vi| < 2ch\\^M^f^l+ - VM^|U.(o,T;L2(nxD)) W^M 



L'^ (a,T;L°°(nxD)) 



- lire,L ^''"L^(0,T;Lij(nx_D)) 'L ^^"L^(0,T;L°°(nxD))' 



where we also used the elementary inequality l-^ci — ^/cJI < \/\ci — C2I with ci,C2 G M>o- The 
norm of the difference in the last displayed line is known to converge to as L — !• 00 (and At 0+) 
by (j6.40dp . This then implies that the term Vi converges to as L — 00 (and At —> 0+). 
Concerning the term V2, we have that 

V2 = 2e/ / M^V.J ■ J M 4 V,(? dq dx dt. 
Jq JnxD ~ * ~ ^ ~ 



Once we have verified that y M 'ipeYx'P belongs to L^(0,T;i^(r2 x £))), the weak convergence 
result (|6.40ap will imply that 



V2^2e [ I MiV^\J^]^- ^ M^^V^ip dqdxdt 
Jo JnxD 

= e I ( MVxi)s- Vx^M<^xdt 
Jo JnxD ~ ~ ~ ~ 

as L — > 00 (and At —>■ 0+), and we will have completed Step 3.2. Let us therefore show that 

\J M ips Yx'f belongs to L^{0, T; L^{^1 x D)); the justification is quite straightforward: using (|6.47p 
we have that 

pT I- I i-T 



f f lyM^V^^pd-Tdt = f f M^i,\V^^\^dqdxdt 
Jo JnxD ~ ~ Jo JnxD ~ ~ ~ 



< 



I I |lV,^||i..p)da;dt 
Jo Jn ~ 

= llYx</?|li2(o,T;L2(0;L°°(D))) < 

That now completes Step 3.2. 

Step 3.3. The third term in (|4.4I) is dealt with as follows: 



\ \ Mu^^l V'^L^ ' ^ x'f' dq dx dt = — / / AIu^^p^- Si'x'^dqdxdt 
Jo JnxD ~ ' ' ~ ^ ~ JnxD ~ ~ ~ ~ 

+ / M (u^ - uf*£^)il)fl^'^ ■ V^ip dqdxdt 
Jo JnxD ~ ~ ' ' ~ ^ ~ 

T 



JnxD 



M Ue ($e - V^i^i'^) • V^^ dq dx dt. 



We label the last two terms by V3 and V4 and we show that each of them converges to as i — > 
(and At — >■ 0+). We start with term V3; below, we apply Holder's inequality with r e (l,oo) in 
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the case of c? = 2 and with r E (1, 6) when d — 3 



IV3 



< 



Jn 

T 



( At.--, 



Jn 

T 



At.- 



Vx^pj dq dxdt 
I|V2:^||l^(d) dx dt 



< / lue-uf^' \\\Vx(p\\L'^{D)dxdt 
Jo Jn ~ ~ 



< 



< 



At,- 



dx 



< 11^ - ^^L"llL^(0,T;L-(n)) \\Y-^h-iO.T;L^in;L^iD))y 

where in the transition from the second hne to the third hne we made use of (j4.31bp . Thanks to 
()6.39cp the first factor in the last hne converges to 0, and hence V3 converges to also, as L ^ 00 
(and At ^ 0+). 

For V4, we have, by using Fubini's theorem, the factorization 



L°°(D) 

dt 



dx 



dt 



(6.52) 



M (4 - ^tr) 



= M3 



M5 



together with the Cauchy-Schwarz inequality, (j4.31bl) . (|6.47l) and the inequality \y/ci — -^/cj I — 
\/|ci — C2I with ci,C2 e M>o, that 



|V4|< 



M \u,\ - |V,^| dq dx dt 

M l^^e - dg I ||Va;^||L~(D) dxdt 



< 



JnxD 

T 



Jn 

T 



< 2 



< 2 



< 2 



Jn 

T 



D 



\U,\[I M\^^ij,~jipf}: + \^dq] \\Vxf>\\L^iD)dxdt 



lueP da; 



"e||L2(0) 



nxD 



X ||V:r^||L~(OxD) dt 



nxD 



M\iPe-$tL^\dq dx 



X \\Vx$\\L'^{nxD) dt 



By (I6.39ap the first factor in the last line is finite while, according to (j6.40dp (with p = 1), the 
middle factor converges to as L — >■ 00 (and At — 0+). This proves that V4 converges to as 
L —i' 00 (and At 0+), also. That completes Step 3.3. 

Step 3.4. Thanks to (|6.40bp . as L ^ 00 (and At 0+), 



M3 q\J i:f^^'^^ A/5 VgxJ^^, weakly in L^{Q,T;L^{n x D)). 

This, in turn, implies that, componentwise, as L — 00 (and At — > 0+), 

^1, \J^^f!t^ ^'^^ weakly in ^^(o, T; L^{^ x D)), 
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for each j = 1, . . . ,K, whereby also, 



K 



K 



3 = i 



for each i = 1, . . . ,K. That places us in a very similar position as in the case of Step 3.2, and we 
can argue in an identical manner as there to show that 



K K 



kl I E E ^^i'^ ■ ^ dx di 

^ / / E E ^9.^= ■ ^s."? dxAt 

"^■^ Jo JQxD ^^-^ ~ ~ ~ ~ 



as L — > cx) and At 0+, for all ip G L'^{0,T; W'^-°°{n x £>)), r e (2,cx)), and in particular for 
all if e C\[0,T];C°°{n x D)). That completes Step 3.4. 

A'tep 5.5. The final term in (|4.4p . the drag term, is the one in the equation that is the most 
difficult to deal with. We shall break it up into four subterms, three of which will be shown to 
converge to in the limit of L ^ oo (and At 0+), leaving the fourth term as the (expected) 
limiting value: 



.53) 



r [ M V[a(7.fi'+)<z,]/3^(V^fi'+) • V,,^dQda;dt 

Jo JflxD j^-^ ~ ~ ' ~ ' ~ ~ ~ 

Jo JnxD j^j^ 

-- / ME 

Jo JnxD ,_i 



V7 At,+ \ 



Jo JnxD 



i=l 
K 



V7 At,+ \ 



V7 1 



i=l 
K r 



Jo JflxD -1 



(/3^(V^S'^)-/3''(V^e)) • V,.^d<7 
(^^(^e)-^s) • V5,(?dqdxdt 
T/jg • Vq . 1^ dg da; dt 
ipg ■ Vq.(p dq dx dt. 



dx dt 



Strictly speaking, we should have written "L — >■ oo, with L g as in Step 1 above, instead 
of "L — >■ oo" ; for the sake of brevity we chose to use the latter, compressed notation. The same 
notational convention applies below. 

We label the first three terms on the right-hand side by V5, Ve, V7, respectively. We shall show 
that each of the three terms converges to 0, leaving the fourth term as the limit of the left-most 
expression in the chain, as i — > 00 (and At — > 0+). 

We begin by bounding the term V5, noting that is Lipschitz continuous, with Lipschitz 
constant 1, writing, as before b := \b\i, and using the factorization (|6.52p together with (j4.31bp 
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and (|6.47p . and then proceeding as in the case of term V4 in Step 3.3: 



IV5I <Vb [ [ M|V 



At,+ 



At,+ | I 7At,+ 



< \/6 



-(/^el I Vg(^| dq dx dt 



M 1 7^, 



At,+ 
e.L 



■06 1 dg I da; 



|V„<?| 



L=°(f2x_D) 



dt 



< 2V6||V.uf*'+| 



X llYq^lU°=((0,T)xnxD))- 

On noting the bound on the third term on the left-hand side of (|6.35p and the convergence result 
(|6.40dp that was proved in Step 2, we deduce that term V5 converges to as L — >• 00 (and 
At 0+). 

We move on to term Ve, using an identical argument as in the case of term V5: 

i-T 



iVel <^/b [ f M\V^: 
Jo JnxD ~ ' 



ufl+\\P'^{i>e)-^e\\y,^\ dqdxdt 



< Vb 



IV7 At,+ | 



I / M|/?^(V'e)-^e| dq] dx 



D 



Vg(^||L~(nx_D) dt 



< 2Vb\\y,ufl+\\mo,T;LHn)) Wl^'^i^e) - $e\\l,^o,T;Ll,{nxD)) 

^ ll'^9<?llL°°((0,T)xnxD))- 

Observe that < tpe — f^^i^e) < V'e and that — f3^{'ip^) converges to almost everywhere on 
n X D X {0,T) as L 00. Note further that, thanks to (I6.40dp with p = 1, V^^s G L^iS^,T\L\i{n x 
D)). Thus, Lebesgue's dominated convergence theorem implies that, as L — >■ 00, the middle factor 
in the last displayed line converges to 0. Hence, recalling the bound on the third term on the 
left-hand side of (|6.35p . we thus deduce that Vg converges to as L 00 (and At — > 0+). 
Finally, we consider the term V7: 



V7 := - 



M 



JnxD 



K 

E 

=1 



(y 



A*,-|- 



V^rUj qi 



ipe ■ '^qi'f dq dxdt. 



We observe that, before starting to bound V7, we should perform an integration by parts in order 



to transfer the x-gradients from the difference V^; uf'^£~^ — onto the other factors under the 



integral sign, as we only have weak, but not strong, convergence of Va, y^^^ — Va, to 0, (cf. 



(|6.39bp ) whereas the difference 



At,+ 
ie,L 



Ue converges to strongly by virtue of (I6.39cp . 



We note is this respect that the function x G fl 1-^ i^s{x, q,t) £ M>o has a well-defined trace on 
dfl for a.e. {q,t) £ D x (0,r), since, thanks to (|6.40ap . 



,P,i;q,t)eH\n), 



and therefore 



for a.e. {q,t) e D x (0,T), implying that ^ij,{-,q,t) \on G L^P{dfl) for a.e. {q,t) e D x (0,T), 

with 2p £ [1,00) when d = 2 and 2p £ [1,4] when d = 3, whereby ipeldn £ i^'(9i7) for a.e. 
{q,t) £ D X (0,T), with p £ [l,oo) when d = 2 and p £ [1,2] when d = 3. As the functions 

and uf'*^^ have zero trace on dV,, the boundary integral that arises in the course of integration by 
parts is correctly defined and, in fact, vanishes. With these preliminary remarks in mind, we first 
write 

»T ^ K d 



V7 



JnxD 



.^ ^E E 



d 



i—l m,n— 1 



dXji 



-ipe {'^q,'P)n dq dxdt. 



Here, {uf'^£~^)n and {ue)n denote the nth among the d components of the vectors uf'^£^ and Wg, 
1 < n < d, respectively, and {YqiV)n denotes the nth among the d components of the vector Yqi'P^ 
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1 < n < d, for each i e Similarly, {qi)m denotes the mth component, 1 < m < d, 

of the d-component vector qi for i e {1, . . . ,K}. Now, on integrating by parts w.r.t. Xm and 
cancelling the boundary integral terms, with the justification given above, we have that 



V7 



/o Jn^D 



/o JnxD 

i-T 



i—l m,n—l 

K d r 

, , E (( 

=: V74 + ¥7,2- 
For the term Yt^i, we have that 

»T K d 



' i—l m.n—1 

K d r 



(Vgi(^)„) dq 
■ iy qi'p)n dg da; di 



da; dt 



9a;„ 



(Vg,(^)„ I dq Axdt 



IV 



7,1 



< 



JnxD 



E E k"^i^)" ~ "^^f^)" \i<l^)rn\ 



< 



I-T 

lo JnxD 



i—l m,n— 1 



K 



dXr^ 



dq dx dt 



,^ME E (; 



i—l m,n— 1 



9Ve 




dXr, 



dq dx dt 



M \uf^{^ - u,\ kllV.^el |V,(^| dg da;dt 



/o JnxD 
< Vb 



2Vb 







M \u 



nxD 



At.+ 

■e,L 



^x^el dq dx I II Vgiy9 IIloo(OxD) 

dt 



T r 



< 2%A 



At.H 



At.H 



MV^/'elV.V^eldg da; 



l|Vg(^ IIl^(oxD) dt 



MlV^V^^rdg da; 



\\^ q^\\L=°(nxD) dt, 



where in the transition to the last line we used the Cauchy-Schwarz inequality in conjunction with 
the upper bound (|6.47l) . Hence, 



|V7,i| <2V6| 



At,+ 



^e|U2(0,T;L2(O)) '^^\\L^(Q,T:Llj(nxD)) 

X II Vg'P IIl°°(0,T;L°°(OxD))- 



Thanks to (|6.39cp with r = 2 and (|6.40ap . ¥7,1 tends to as L ^ (and At ^- 0+). 

Let us now consider the term ¥7^2- Proceeding similarly as in the case of the term ¥7^1, using 
(|07|) . yields 



|V7^2| <Vb\ 



At,+ 



Noting (|6.39cl) with r = 2, we deduce that ¥7^2 converges to as L — > (and At — > 0+). Having 
shown that both ¥7,1 and ¥7,2 converge to as L ^ (and At — >■ 0+), it follows that the same 
is true of V7 = ¥7,1 + ¥7^2- We have already shown that V5 and Vg converge to as L — > (and 
At — > 0+). Since the sum of the first three terms on the left-hand side of (|6.53p converges to 0, it 
follows that the left-most expression in the chain ()6.53p converges to the right-most term, in the 
limit of L — >■ 00 (and At — >■ 0+). That completes Step 3.5. 
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Having dealt with (|4.4[) . we now turn to (|4.3|) . with the aim to pass to the hmit with L (and 
Ai). In Steps 3.6 and 3.7 below we shall choose as our test function 

w e C^{[0,T];C^{n)) with w{-,T) ^ 0, and V ^ ■ w ^ on fl for all t e [0,T]. 

Clearly, any such w belongs to i^(0,T; V) and is therefore a legitimate choice of test function in 
()4.3p . Furthermore, for any cr > 1, the set of such smooth test functions w is dense in the space 



of all functions in W^'^{0,T;Vcr) that vanish at t = T. As each term in (|4.3p has been shown 
before to be a continuous linear functional on L^{0,T; Va), o > ^d, a > 1 and W^'^{0,T; Va) is 
(continuously) embedded in L^{0,T;Va), cr > ^d, a > 1, the use of such smooth test functions 
for the purposes of the argument below is fully justified. 

Step 3.6. The terms on the left-hand side of (|4.3p are handled routinely, using (|6.35l) and, 
respectively, integration by parts in time in conjunction with (j6.39cl) with r = 2, (I6.39b[) and 
recalling that — > uq weakly in H . In particular, the second (nonlinear) term on the left-hand 

, y w) dt, and then 



side of (|4.3|) is quite simple to deal with on rewriting it as — /p (; 
considering the difference Jq {ue (S" 



At,+ 



(g) u 



At,- 



^'e'l^ ^ h'c'l ' '^^i which is bounded by 



ie.L 




I ^ Ai,+ ^ ^c,— II 



At, -I 



dt WV^wl 



L°°(0,T;L<=°(n))- 



By adding and subtracting uf'*£~ inside the first norm sign, using the triangle inequality, 
followed by the Cauchy-Schwarz inequality in each of the resulting terms, and then applying the 
first bound in (|6.35p . and (|6.39cp with r = 2, we deduce that the above expression converges to 
as L — > oo (and At — > 0+). The convergence of the first term on the right-hand side of (14. 3p to 
the correct limit, as L oo (and At 0+), is an immediate consequence of (|3.23p . We refer the 
reader for a similar argument to Ch. 3, Sec. 4 of Temam [35]. That completes Step 3.6. 

Step 3.7. The extra-stress tensor appearing on the right-hand side of (|4.3p is dealt with as 
follows. First, by using p.l6p and noting that w is, by assumption, divergence-free, and proceeding 



At,+ 



in exactly the same manner as in (|4.48p , but with ip^ ^ 
that 



now replaced by V'^l'^ 



ipe, we have 



Vs 



k 



T K 

E 

1=1 



V X w dxdt ^ k 



T K 

E 

i=l 



f I I ^^E(v. 

Jo Ju JD ,_i ~ 



w)qt 



■4^e)dq 



Ct{M ^e) ■■ Vxwdxdt 



dx dt 



We rewrite the second factor in the integrand of the last integral as follows: 



M^ Vq^V'/'. 



Hence we obtain the following inequality: 



Vs < 2kVb f f \Vxw\ f 
Jo Jil ~ ~ Ud 



At, 



+ 2k 



Jo Jn \_Jd ^^-^ » ~ 



M2 V,J^^ 



dq 



dx dt 



M5 V. 



dx dt 



Vs,! + Vg 



We emphasize at this point that in the term Vs,2 wc intentionally did not move the modulus sign 
under the integral: we shall be relying on the weak convergence result (I6.40bp to drive the term 
Vs,2 to in the limit of L — > oo (and At 0+), so it is essential that the modulus sign is kept 
outside the integral. 
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For Vg.i, we have, by using that | — ^Jc^\ < \/\c\ — C2I for any ci, C2 € IR.>o 

Vg.i <2fcV6 / ||V,ii;|Uoo(o)||^^fl'+-^;e|li 
Jo ~ ~ ' " 

< 2fc\/6||V3; U;||Loo(o,T;L==(f2)) 



(JlxD) 



'^^l'''IU2(o,T;L2^(OxD))- 



we deduce that the term 



By noting (j6.40dp with p ~ 1 and the bound on the seventh term in 
Ys,i converges to in the hmit of L -> 00 (and At — > 0+). 

Finally, for Vs. 2, we first define the (A'(i)-component column- vector function S := [Sj , . . . , 

where S, := 



"T IT 



, and note that 



Jn 



K 



2k 



i=l 



M5 



JnxD 



i/'e I • S dg dx dt 





dx dt 







The convergence of V8,2 to will directly follow from (I6.40bp once we have shown that S e 
L^(0,T;L2(f2 x D)). The latter is straightforward to verify, using (lOT)) : 



iSpdg da;dt = 



nxD 



E 

1=1 

E 



iSjpdq da; dt 



nxD 



Mips (V:r 



dq dx dt 



< 



K 

E 

1=1 



|V:,wp|(j,pMV£ dg do; dt 



nxD 



AI ips dq dx dt 



where in the transition to the last line we used (|6.47p . 

Thus we deduce from (j6.40b| that ¥§,2 converges to as L — )• cx) (and At — 0+). As both 
V8,i and V8,2 converge to 0, the same is true of Vs, which then implies (|6.40ep . thanks to the 
denseness of the set of divergence- free functions contained in C^([0,T]; Cjf (fi)) and vanishing at 
t = T in the function space L^(0,T; Va), u > ^d, a > 1. That completes Step 3.7, and the proof 
of dOOel) . 

Step 3.8. Steps 3.1-3.7 enable us to pass to the limits L — ^ cx), At 0+, with At = o{L^^) 
as L — 7- cx), to deduce the existence of a pair {ue.ipe) satisfying (|6.4ip . (|6.42p for smooth test 
functions ^ and w, as above. The denseness of the set of divergence-free functions contained 
in Ci([0,T];C'g°(f})) and vanishing at t = T in the set of ah functions in W^^^{id,T;V and 
vanishing at t = T, tr > ^d, cr > 1, and the denseness of the set of functions contained in 
Ci([0,T]; C^ip. X D)) and vanishing at t = T in the set of all functions in W'^'^{Q,T]H'{VL x D)) 
and vanishing at t = T, s > 1 + i(i\: + l)(t, yield (|OT|) and ((O^ . That completes Step 3.8. 

Step 3.9. Let X be a Banach space. We shall denote by Cuj([0,T];X) the set of all functions 
X G L°°(0,T;X) such that t e [0,T] ^ {x',x{t))x G M is continuous on [0,r] for aU x' G X' , 
the dual space of X. Whenever X has a predual, E, say, (viz. E' — X), we shall denote by 
C^*([0,T];X) the set of ah functions x G L°°(0,T;X) such that t G [0,r] ^ {x{t),e)E G E is 
continuous on [0, T] for all e ^ E. 

The next result will play an important role in what follows. 
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Lemma 6.3. Let X and Y be Banach spaces. 

(a) If the space X is reflexive and is continuously embedded in the space Y, then L°°{0,T; X)r] 
C^{[0,T];Y)^C^{[0,T];X). 

(b) // X has separable predual E and Y has predual F such that F is continuously embedded 
in E, then i°°(0, T; X) n C^* ([0, T]; F) = C^* ([0, T]; X). 

Part (a) is due to Strauss [48] (cf. Lions & Magenes [34], Lemma 8.1, Ch. 3, Sec. 8.4); part 
(b) is proved analogously, via the sequential Banach- Alaoglu theorem. That G Cw{[0,T]; H) 
then follows from € L°°{0,T]H) Ci H^{0,T;V'„) by Lemma [Oi;a), with X := H,Y ;= V'„, 
cr > id, cr > 1. That G C^i[0, T]]LIj{Q X D)) follows from Ti^!^) e L°°{0, T; L\;{n x D)) and 
V^e e H^O^T-^M-^H^n X D))') by LemmaHSJJb) on taking X := i* (17 x D), the Maxwellian 
weighted Orlicz space with Young's function <l>(r) = J"(l + |r|) (cf. Kufner, John & Fucik [31], Sec. 
3.18.2) whose separable predual E :— Efj{Q x D) has Young's function 5'(r) = exp \r\ — |r| — 1, and 
Y := Af~i(i7''(f7xi:)))' whose predual w.r.t. the duality pairing (Af-, ■)H^(nxD) is F := H'^{nxD), 
s > 1 + i(ii: + l)c;, and noting that ([0; T]; L* (f7 x D)) ^ C^,{[0,T]; Llj{n x D)). The last 
embedding and that F ^ E are proved by adapting Def. 3.6.1. and Thm. 3.2.3 in Kufner, John 
k Fucik [31] to the measure M{q) dq dx to show that L°°{fl x D) ^ x D) for any Young's 

function S, and then adapting Theorem 3.17.7 ibid, to deduce that F ^ L°°(r2 x D) ^ EfAil x 
D) = E. [The abstract framework in Temam [?5], Ch. 3, Sec. 4 then implies that and satisfy 
u,{-,0) = Mo(-) and V'e(T,0) = V^(-,-) in the sense of C^([0, T]; i/) and C^([0, T]; Li,^(17 x L>)), 
respectively.] 

The energy inequality (|6.43l) is a direct consequence of (|6.39al -c) and (j6.40a[ b.d'). 
on noting the (weak) lower-semicontinuity of the terms on the left-hand side of (|4.28l) and (|6.48p . 
That completes Step 3.10. 

Step 3.11. It remains to prove (|6.37p . The bounds on the first and third term on the left- 
hand side of ()6.36p imply that the sequences {/3^^''^''}l>i are bounded in L°°{0,T]L°°{VL)) n 
L2(0,T;iJi(fi)); the bound on the fourth term in ([Oel) yields that {p^l} i>i is bounded in 
ffi((0,r); (i7i(r2))'). In fact, by noting (ESj) and dHj), we have that {pf^t^^}L>i is a bounded 
sequence in L^{0,T; IC). Thus, there exist subsequences of {p^*^'^^}l>i (not indicated) with 
At = o{L~^) and, thanks to the uniform bound on the second term on the left-hand side of (j6.36p . 
a common limiting function G L°°(0, T; L°°{n)) n L'^{0, T; JC) n H^{0, T; (i/i(fi))') such that 

(6.54a) pf'j^''^'> p, weak* in L°°(0,r; L°°(17)), 

(6.54b) P^l'"^^ ^ Pe weakly in L^{Q,T; H^fl)), 

(6.54c) ^ ^ ^ weakly in ^^(o, T; (ff ^(r!))'), 

as L — > oo (and thereby A< 0+), It follows from (|6.54bl c) and Lemma 1.2 in Ch. 3 of Temam 
[IS" that belongs to C([0, T]; ^^(il)), in fact. Using (|6.54aL c) and (|6.39cp . we can pass to the 
limit as L ^ oo (and At 0+) in (H3)) to obtain that, for ah G L^{0, T; i?H^^)), 

(6.55) / dt+ / / [eVxPe-UePel • V:,¥>da;dt = 0. 



"'0 \ 9t /Hi{n) Jo Jn ^ ~ ~ 
We note also that Fubini's theorem, (|4.6p and (j6.40dp yield that 



r f pfi- ( M^.dq dxdt^ r 1 1 

Jo Jn JD ~ ~ Jo Jn JD 



M{i^^X-i^,)dq 



dx dt 



56) -11 1-0^1 -V'eldg da; di^O as L c» (and At ^ 0+). 

Jo JnxD ' ~ ~ 
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Thus, pf^^ M dq strongly in ^^(O, T; L^ifl)) as L cx) (and At 0+). Comparing this 

with (|6.54ap imphes that 



(6.57) Peix,t) ^ / M{q)ijj,{x,q,t)dq for a.e. (x, t) e 17 X (0, T). 

It foUows from Step 3.9 that, for s > 1 + ^{K + l)d, we have 



hm / M{q){^lj^{x,q,t) -ipo{x,q))^{x,q)dqdx = yip€H'{nxD). 
*-*0+JOxD ~ ~~ ~~ ---- 

Consequently, using (|6.57p and p.Sp we then deduce by selecting any (p — Lp ^ C°°{^) (= C°°{ri)(i 
1{D) C H%VL X D)) that 

lim / pg{xTt) ip{x) dx — lim / / M[q)'ij;g[x^q,t) dq] (p{x) dx 



= lim / M{q) ^pg{x,q,t) (p{x) dq dx — / M{q) 'ipo{x, q) ip{x) dq dx 
*^"+ JnxD ~ ~ ~ ~ ^ ~ JuxD ~ ~ ~ ~ ^ ~ 

(6.58) = / ( / M{q)Mx,q))dq]ip{x)dx^ f ip{x)dx. 



As pe G C{[0,T]-L^{n)), it follows from that Pe{x,0) = 1 for a.e. x e n. 

Clearly the linear parabolic problem (|6.55p with initial datum p^{x,0) = 1 for a.e. x & ^ has 
the unique solution p^ = 1 on x [0,T]. Using this in (|6.57p implies (|6.37p . and completes Step 
3.11 and the proof. □ 

7. Exponential decay to the equilibrium solution 

We shall show that, in the absence of a body force (i.e. with / = 0), weak solutions (j/e,V'e) 
to (Pe), whose existence we have proved via our limiting procedure in the previous section, decay 
exponentially in time to the trivial solution of the steady counterpart of problem (P^) at a rate 
that is independent of the specific choice of the initial data for the Navier-Stokes and Fokker- 
Planck equations. Our result is similar to the one derived by Jourdain, Lelievre, Le Bris & Otto 
|29| . except that the arguments there were partially formal in the sense that the existence of a 
unique global-in-time solution, which was required to be regular enough, was assumed; in fact, the 
probability density function was supposed there to be a classical solution to the Fokker-Planck 
equation; -ipo was required to belong to L°°{fl x D) and to be strictly positive, and u was assumed 
to be in L°°{0,oo;W'^'°^{n)) (cf. p.l05, (B.128), (B.129) therein; as weh as the recent paper 
of Arnold, Carrillo and Manzini IS] for refinements and extensions). In contrast, we require no 
additional regularity hypotheses here. 

Theorem 7.1. Suppose the assumptions of Theorem \6. 1\ hold and M satisfies the Bakry-Emery 
condition (cf. Remark \5.1\) with k > 0; then, for all T > 0, 

lhe(T)||^-fA||4(r)-l||i.^(^,^) 

(7.1) <e-''oTUuof + 2k f MJ^{^^o)dqdx]+- r imf^^.^^yds, 
where 70 := min , -^^j ■ In particular if f = 0, the following inequality holds: 

(7.2) ||«,(T)||2-K A||4(T)-1|12 <e'^°^Lof + 2fc / M J-(v7o) dq dx 

l"l M^ ' ^ ^ JnxD ~ ~. 

Proof. We take t — ti ~ At and write = to in (|4.23l) . and we replace the function T on the 
left-hand side of (|4.23l) by T^, noting that, prior to (|4.23p . in (|4.19l) we in fact had on the 
left-hand side of the inequality, and J-^ was subsequently bounded below by J-; thus we reinstate 
the J^^ we previously had. We recall that m° = y^L~(^i) and 13^ (tp^) — V'i^i (^1) ^nd adopt the 
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notational convention t_i := —oo (say), which allows us to write ^^^^(^o) instead of u^^' (ti) 
and ipf'^j^^ {Iq) instead of "0^^' (^i)- Hence we have that 

+ 2k f AfJ-^(^fi^+(ti) + a)dgdx+-A. / Y M{^fl+ -^f'j^-)'dqdxds 

JflxD ~ ^T^^Jto JnxD 

Jto JnxD £ + a ~ Jto JnxD 4'e l + 

(7.3) +2fc / Af J"^(^/;^^'+(to) +a)dgdx. 



1 



<IICr(Mr + - / lir '^(s)li(/^„Hn))'d5 



Closer inspection of the procedure that resulted in inequality (|4.23p reveals that (j4.23p could have 
been equivalently arrived at by repeating the argument that gave us (|7.3p on each time interval 



[t„_i,i„], n = 1, . . . ,iV; viz. 



I ^*.+ /'+ M|2 I / II At,+ Ai.-||2j 

2A6fc 



/" lly.yfi+(s)||'ds+ 2fc / M^^(V^f[^+(t„)+«)dgdx 
Jt„_i ' JnxD 

r I Mii^f'^ -^^f'^fdqdxds 
Jt„_i JnxD 



AtL 



< 



^t.+ f-^ M|2 I / II /•At,+ / M|2 



'2A;e / / M Z.^^ , dqdxds+ ^ / ^It + dgdgds 

t„-iJnxD Ve.L + Q; ~ ~ Jt^^iJnxD 4' ' 



<r(^-i)ir + - / iir'^(^)ii(H„Mn))'d^ 



(7.4) +2fc / Af J-^(v!^f^'+(t„_i) +a)dgdx, n = l,...,A^, 

summing these through n and then bounding on the left-hand side below by F. 

Here we proceed differently: we shall retain on both sides of (|7.4p . and omit the second, 
fifth and sixth term from the left-hand side of (|7.4p . Thus we have that 



2k [ M (V^f (t„ ) + a) dg dx 
JnxD ~ ~ 



< 



llyti'+(^n-l)ir + ^ / ■^||/^*'+(.)||^^.(^)yd. 



tn-l 



(7.5) +2k I MT^{ii^^'£+{tn-i) + a)dqdx, n^l,...,N. 

JnxD 

Thanks to Poincare's inequality, recall (|4.34p . there exists a positive constant Cp ~ Cp(r2), such 
that 

(7.6) ||<i'+(.,.)|| <Cp(f7)|lV.<i^+(.,.)|| 
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for s £ {t„-i,tn]] n = 1, . . . ,N. Also, by the logarithmic Sobolev inequality (|5.4p . we have for a.e. 
X £ fl that 



L 



M{q) [7Afl'+(x,q,s)+a] log 



dq 



< - / M(g) 
1^ Jd 



dq, 



for s € (t„_i,t„]; n ^ 1, . . . , N . Hence, for a.e. x £ il, 



M{q) [^)f^£^{x, g, s) + a] log[-0^['+(x, g, s) + a] dq 



< - / M{q) 



^q\l'4'fx'{x,q.s) + a 



dq 



(7.7) + ( / M{q) [i:f'j^+{x,q,s) +a\dq] log ( / M{q) [V'f 1'+ (x, q, s) + a] dq 



D 



D 



for s S (t„_i,i„], n ~ 1,...,N. Note that, thanks to (j4.31bp and the monotonicity of the 
mapping s E M>o logs G M, the second factor in the second term on the right-hand side of 
(17. 7p is < log(l + a). Since a G (0, 1), we have log(l + a) > 0; also, the first factor in the second 
term on the right-hand side of (|7.7p is positive thanks to (I4.31ap and by (j4.31bp it is bounded 
above by (1 -I- a). Hence the second term on the right-hand side of (17.71) is bounded above by the 
product {1 + a) log(l -I- a). We integrate the resulting inequality over 51 to deduce that 



M{q) [■ipf'*£^{x,q,s) + a] \og[iljf V^{x,q,s) + a] dq dx 



nxD 



2 

< - 

K 



nxD 



M(q) |V, g, s) + dq dx + \Q\ (1 + a) log(l + a). 



for s G (tn-i,tn], n — 1,...,N. Equivalently, on noting that slogs — J^(s) — (1 — s), we can 
rewrite the last inequality in the following form: 



nxD 



M{q) {x, q, s) -|- a) dq dx 



< 



^ JnxD 



M{q) \V,,hf'l+(x, q, s) -I- a|2 dq dx 



(7.8) +/ M{q){l~^P^l:+{x,q,s)-a)dqdx + \n\{l + a)\og{l + a), 

JnxD ~ ' ~ ^ ^ ~ 



for s G {tn-i, tn], n = 1, . . . , N . This then in turn implies, thanks to the fact that ipf'*£'^ {x, q, •) is 
constant on the interval (t„_i,t„] for all {x,q) G H. x D, that 



At 



Af(g)J-(Vfi'+(i„)+a)dq dx 



< 



2ao k 
A 

Kank 



nxD 



t„_i JnxD 



M lyqJipfl:'^ + dq dx ds 



M (1 - V^fi'^ - a) dg dx ds + At (1 + a) log(l + a) 



nxD 
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for n = 1, . . . ,N. Using this and (j7.6|) in (|7.5p then yields 



At f 2Xbk\\ II At.+ . N||2 



i+T^i'^-"— ) )iKr(ui 



' JnxD 

^ 2k f M [J'Hi^tL^itn) + «) - -^(^^^l'+(tn) + «)] dq dx 



< 



inxD 

||yfi'+(in-i)||' + 2fc / MJ-(^f[^+(t„_i)+a)dgds 

+ 2fc / Af [J-^(^f[^+(i„_i) + a) - J-(V^f['+(t„_i) + a)] dq dx 

^naok r Mil-$f'j:+-a)dqdxdt 

't„-i JnxD 



X 



(7.9) +^Ai|J7|(l + a)log(l + a) + i jj'*" 

for n = 1, . . . , iV. We now introduce, for n = 1, . . . , N, the following notation: 

. / 1 / 2Xbk\ Kao\ 
7(a) :=mm(^^ 

A„(a) ||ufi^+(t„)ir + 2fc / M J-(^fi^+(t„) + a)dgdx, 

nixD 



B^{a) := 2fc / M [-F^(V^fi'+(t„) + a) - J-(V!^fi'+(i„) + a)] dq dx, 

QxD 



Ka^k 



Cn{a):='-^ / / M{l^^f'£+{x,q,s))dqdxds 
^ Jt„-i JnxD ' ~ ~ ~ ~ 

+ 1^ At \n\ [(1 + a) log(l + a) - a] + i ||/^*'+||^^.(^)), d.. 

We shall assume henceforth that a is sufficiently small in the sense that (|4.22p holds. For all such 
a, 7(ck) > 0; further, trivially, An{a) is nonnegative; by ()4.12p . we have that i3„(a) is nonnegative, 
and by (|6.4p and since J-{1 + a) > for all a > 0, C„(a) is also nonnegative. In terms of this 
notation (17. 9p can be rewritten as follows: 

(1 + 7(a) At) An{a) + B,,{a) < A„_i(a) + B„_i(a) + C„(a), n = 1, . . . , iV. 

It then follows by induction that 

n 

An{a) < (l+7(a)Ai)""Ao(a)+^i^,(a), n= 1,...,7V. 

i=i 

That is, 



An{a) + Bnia) < (1 + 7(a) Ai)-"Ao(a) + \ Bo{a) + ^ (a) , n = l,...,N. 
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In particular, with n = iV, by omitting the nonnegative term B]^{a) from the left-hand side of 
the resulting inequality, and recalling that T ^ = NAt, we get 



At, 



i -(r)ll^ + 2fc / MT{i^f'£^{T) + a) dq dx 



< 1 



7(a) T 



N 



(0)||^ + 2fc / MJ-(^f['+(0) +a)dgda; 

flxD ' ~ ~. 



2k 



M 



nxD 



dq dx 



Kagk 



JnxD 



M (1 — ip^*£^{x,q, s)) dq dx ds 



(7.10) 



+ ^TN^(l + a) + i A/ 



At,+ ||2 1 



Using that uf^'+(0) = w° and = P^i-^^)-, we then obtain from ([7101) that 



||ufi^+(r)|r + 2fc / Af^(V.fr(T) + a)d<?dx 



< 1 



7(a) T 



SlxD 



\U 



0||2 



2/c 



MJ"(/3^«) + a) dg da; 



2k 



M 



QxD 



InxD 

J-^(/3^(^")+a)- J-(/3^(V^0) + a)' 



dq da; 



+ 



(7.11) 



Kttok 

A 

Kacjk 
A 



JnxD 



M (1 — ■0^^'^(a;, s)) dq dxds 



r|f7| J-(l + a) 



ll/'^*'^ll(Hi(o))'ds- 



Applying (I4.12p and (|4.25p in the second factor in the first term on the right-hand side of (|7.1ip 
and using (I4.24p in the square brackets in the second term on the right-hand side, we have that 



|<r(r)ir + 2fc / MT{^CL{T) + a)dqdx 
~ ' JnxD ~ ~ 

-N r 



< 1 



(7.12) 



N 

+ 3ak\n\ 



\u"\\^ +3ak\n\ +2k 



MJ'i^^+a) dq dx 



nxD 



JnxD 



M (1 - ipf'£+{x,q,s)) dq dx ds 



A 



T\n\T{l + a) + - f 11/ 
Jo ~ 



At,+ ||2 J 



We now pass to the limit a — > 0+, with L and At fixed, in much the same way as in Section |4TT 
Noting that limQ_5.o^ 7(a) =707 we thus obtain from (17.121) . p.lSp and (I3.20p . that 



\\ufl^{T)r + 2k I MJ-(^^-r(T))d9dx 



At, 



< 1 



10 T 
N 



nxD 

-N r 



2k 



M J"(V'o)dq da; 



nxD 



(7.13) 



KUok 



JnxD 



Af (1 - $^*£^)dq dxds-l- i 



11/ 



At.+ ||2 - 

\\(Hi(n)y as- 



In order to pass to the limits L — >■ 00 and At — >■ 0+ (with At — o{L "'^)) in the first two terms on 
the left-hand side of (|7.13p we require additional considerations. 
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Noting (I6.39cp for the sequence {u^i^}l>i of Theorem 16. 1[ and passing to a subsequence (not 
indicated), as L — > oo and At — > 0+ (with At ~ o{L^^)) we have that ||u^i(i) — ye(OII converges 
to for a.e. t E (0,T); let be one such t in (0, T). It then follows from (|6.39dp that, for any 



\{Ue{T) 



.At 



{T),v 



(7.14) 




.tX{t^)-u,{t^)\\ \\v\\ 



o{L- 



Since : [0, T] — !■ _ff is weakly continuous, we have that u^{T) G H. It follows from the bound on 



the first term in ((05l) . as t G [0, T] ^ 
that, for any vqEH and any v £ Va-, 



,At 
s.L 



{t)\\ G K>o is a continuous (piecewise linear) function, 



\{u,{T)~ufX{T),vo)\<\{u,{T) 



41 



\u,(T)\\+C:^ 



Recalling (|7.14p . it follows from the last inequality that 

limsupi_^|(y,(T) - utl{T),vo)\ <C\\vo-v\ 



As Va is dense in H, we thus deduce that {u'^l{T)}l>i converges to Ue{T) weakly in _ff as L — !• cx) 
and At — >■ 0+, with At = o{L~^). Hence, by the weak lower-semicontinuity of the Li^iU) norm 
and noting that uf\^ (T) 

(7.15) 



At,+ 



^ (T), we have 
|iue(T)|| < liminfL- 



At,+ 



mil. 



Analogously to ()7.14p . noting (j6.40dp for the sequence {V'^£}l>i of Theorem 16. H we have (on 
passing to a subsequence, not indicated,) that ^f'\^{T) converges weakly to ipe{T) in [H'^ {^l x 
D))' as oo and At 0+, with At = o{L-^). Thanks to Theorem [Ol ijj^ : [0,T] L\i{n x 
D) is weakly continuous; hence we have that ipeiT) € L\j{il x D). Similarly to the argument in 
the proof of O of Lemma l6.2[ it follows from the bound on the fourth term in (|6.35p . on noting that 
T{r)/r — >■ oo as r — >■ 00 , together with the de la Vallee-Poussin and Dunford-Pettis theorems, that, 
upon extraction of a further subsequence (not indicated), tpf'\^{T) converges weakly in L\f{^ x D) 
to some limit A, as i — > oo and At 0+, with At = o{L^^). The fact that A — ipeiT) follows 
from the weak convergence of ijj^X^T) to ■0e(T) in M-^{H'{p. x D))' (^ L\[{VL x D)). Finally, 
since r € [0, oo) i— > .F(r) G M>o is continuous and convex, on applying Tonelli's weak lower 
semicontinuity theorem in L\,i{Q, x D) (cf. Theorem 3.20 in Dacorogna !l6 ), 



(7.16) 



MF{i)e{T)) dq dx < liminfL- 



SlxD 



MF{^fX+{T))Aq dx. 



where we have noted that V'i^i(T) = ^pf^^^iT). 

We are now ready to pass to the limit in ([71^ . Using ([71^ and ([71^ . (|6.40dp . ((07)) 

and p.23p . and letting L — > oo (whereby At — >■ 0+ according to At — o{L^^) and therefore 
N = T/At ^ oo), we deduce from (fTTSl that 



2k 



(7.17) 



\\MT)r 



MT{ipe{T))dq dx 



2k 



Af J"(-0o)dq dx 



nxD 



The Csiszar-Kullback inequality (cf., for example, (1.1) and (1.2) in the work of Unterreiter et 
al. [50]) with respect to the Gibbs measure /i defined by d/i — M{q) dq yields, on noting (j6.37p . 
for a.e. x G fi, that 



IIV'e(S,-,r) - 1\\l\,(D) 



< 



MF{'^e{x,q,T))dq 



D 
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which, after integration over fl implies, by the Cauchy-Schwarz inequahty, that 

WMT) - lllii (^,,^) < m [ MFmT))Aqdx. 
" JnxD 

Combining this with (|7.17p yields (|7.ip . Taking / = 0, the stated exponential decay in time of 
{us,ijj^) to (0, 1) in the L'^{il) x Llj{n x D) norm follows from ((7?T|) . □ 

Remark 7.1. By introducing the free energy as the sum of the kinetic energy and the relative 
entropy: 

m :=i||y,(i)f + A:^^^Af J-(4(t))dgda;, 
we deduce from (j7.17p that, for any T > 0, 

€(T)<e-^oM0) + 2^/o'll/(s)ll^^.(^)),ds. 

Thus in particular when f — 0, the free energy decays to as a function of time from any initial 
datum {uo,tpo) with initial velocity uq € H and initial probability density function -00 that has 
finite relative entropy with respect to the log-concave Maxwellian M . 

It is interesting to note the dependence of 70 = min , , the rate at which the fluid 

relaxes to equilibrium, on the dimensionless viscosity coefficient v of the solvent, the minimum 
eigenvalue ag of the Rouse matrix A, the geometry of the flow domain encoded in the Poincare 
constant Cp{fl), the Weissenberg number X, and the Bakry-Emery constant k for the Maxwellian 
M of the model. We also observe that the right-hand side of the energy inequality (j6.43l) and 70 
are independent of the centre- of -mass diffusion coefficient e appearing in the equation (|1.9I) . o 
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Appendix A. Cartesian products of Lipschitz domains 

Let us suppose that D :— Di y. ■ ■ ■ x Dk, where Di, i = 1, . . . ,K, are bounded open balls in 
M'' centred at S M''. Then, D is a bounded open Lipschitz domain (cf. [Appendix C| below for a 
precise definition of a Lipschitz domain). We present two different proofs of this statement. 

The first argument is based on the observation that the Cartesian product of K bounded open 
Lipschitz domains Di C R"', i — 1, . . . ,K, with > 1, is a bounded open Lipschitz domain in 
j^niH hnjf _ First, note that as a Cartesian product of a finite collection of open sets, D itself is 
open. That D is a Lipschitz domain follows by combining Theorem 3.1 in the Ph.D. Thesis of 
Hochmuth [27j , which implies that the Cartesian product of a finite number of bounded domains 
Di C K"S i = I, . . . ,K, with Ui > 1, each satisfying the uniform cone property, is a bounded 
domain in ^"^ satisfying the uniform cone property; and Theorem 1.2.2.2 in the book of 

Grisvard [25) . which states that a bounded open set in R" has the uniform cone property if, and 
only if, its boundary is Lipschitz. 

In the special case of our domain D an alternative proof proceeds by noting that, as a Cartesian 
product of K bounded open convex sets, Di C R*^, i — 1, . . . , K, D is a bounded open convex 
set in R^'^ (cf. Hiriart-Urruty & Lemarechal [55], p. 23), and then applying Corollary 1.2.2.3 in 
Grisvard [3S] , which states that a bounded open convex set in R^'' has a Lipschitz boundary. 

Appendix B. Completeness and separability of L\i{D) and HIj{D) 

The completeness of the spaces H\i{D) and Lfj{D) follows from Theorem 8.10.2 on p. 418 in 
the monograph of Kufner, John & Fucik [31j. 

By Theorem 19 in Section IV.8.19 of Dunford & Schwartz C{D) is dense in LIj{D) = 
L'^{D;ii), where the Radon measure ft is defined on the a algebra of Borel subsets of D by 
Ijl{B) := /g M{q) dq; fj, is absolutely continuous with respect to the Lebesgue measure on R^'' 
and d/x = M{q) dq. Thus, given v G L'\^{D) and any e > 0, there exists G C{D) such that 

\\v-v\\lI,{d) < l^- 

By the Stone- Weierstrass theorem (cf., for example, Pinkus [43 ), C{D) is separable: there exists a 
countable dense set V C C{D), where V is the set of restrictions from R^"* to D of all polynomials 
with rational coefficients (cf. Theorem 1.4.5 on p. 30 of Kufner, John & Fucik [31]); hence, given 
e > there exists p € V such that 

y-p\\c(D) < ¥■ 

Clearly, C(D) C L\i(p) and therefore V C L\j{D); hence, and since fi{D) — 1, 

This shows that the countable set V C L\j{D) is dense in L\.j{D). Therefore L\j{D) is separable. 
By an argument analogous to the one in Section 3.5 on p. 61 of Adams & Fournier [T], H\j{D) is 
separable. 
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Appendix C. Density of C°°{D) in L\i{D) and HIj{D), and of C°°{n x D) in Ll^(n x D) 

AND Hljln X D) 

Since the set P defined in the previous subsection is dense in L\j{D) andV C C°°{D) C L\j{D), 
we deduce that C°°{D) is dense in L\j{D). By an identical argument, C°°(ri x 13) is dense in 
L\j{n X D). 

Next we consider the density of C°°(D) in Hlf{D), and of C°°(17 x D) in x £>). The 

argument below closely follows the proof of Theorem 1.1 on p. 307 in the paper of Necas [39] . 
We suppose that a fixed Cartesian co-ordinate system, referred to henceforth as the canonical 
co-ordinate system, is introduced in R", so that a vector q £ M" is represented in terms its 
co-ordinates as 

q = {qi,...,qn). 

Suppose that A is an n x rt orthogonal matrix with determinant equal to +1 and C is an n- 
component column vector. The affine transformation q^ cf^ := Aq^ + C defines a new co- 
ordinate system, which we shall say is equivalent to the canonical co-ordinate system, and we shall 
express q in terms of its co-ordinates in this second, equivalent, co-ordinate system as 

q = (gi,...,g„). 

We recall the following definition from Kufner, John & Fucik [31] . 

Definition Appendix C.l. We say that a bounded open domain V C M" is a Lipschitz domain 
if, and only if, there exist: 

(i) a positive integer m and m (different) Cartesian co-ordinate systems, referred to as local 
co-ordinate systems, each of which is equivalent to a fixed canonical co-ordinate system in 
M". When an n-component vector is expressed in terms of its co-ordinates in the rth local 
co-ordinate system, we shall write 

Qr = {qrl, ■ ■ ■,qrn) = {q[r) , Q^r)) ' 

where q'^^^ := {qri, . . . ,q'„i-i) and g^^^ := qm- 

(ii) a number a > and m functions 

Or e C°''^(Kr), r = l,...,m, 

where 

(iii) a number /3 > such that 

(111.1) the sets 

Ar := A-\{Qj = {q[r),q(^r)f ■ Q'(r) € Ar and q^;) = ar(g(^))}) 
are subsets of dV for r — 1, . . . , m, 

m 

and Ar : — > QjF is the affine transformation of co-ordinates Qj :— Ar{Q^) = 
ArQ"^ -\- Cr, with Ar G R"^" orthogonal and determinant equal to and Cr an n- 
component column vector, that maps the canonical co-ordinate system in R" to the, 
equivalent, rth local co-ordinate system. 

(111. 2) for every r = 1, . . . , m, the sets and C/+ , defined by 

U~ A-\{Qj = iQ[r)^1{r)f <l[r) ^ ^r, ar{q[r)) - /? < < Ar^^))}) 

U+ := Ar\{Qj = {q[r),q{r)f ■ <l[r) ^ A^, a^(g[^)) < qfr) < Or (g(^) ) + /3}) 

are such that U^T C T> and C V (— the complement of V w.r.t. R"j. 
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We shall write 

Ur A-\{Qj = iq[r),q(r)f ■ <l[r) ^ and a,(q'(,)) ~(3< g^,) < a,(g[,)) + /?}). 
Theorem Appendix C.l. C°°{D) is dense in H\i{D). 

Proof. Let q — {qi, . . . , qk) E D = Di x ■ ■ ■ x Dk, where Di C M'' is a bounded open ball 
in R'' centred at the origin, and let n :— Kd; thus, after relabelling the co-ordinates oi q, q — 
{qi, . . . , qn) e D, and is a bounded Lipschitz domain in M" by [Appendix A In particular, by 



using the notation introduced in the previous definition, there exist m local co-ordinate systems in 
M" and real- valued functions a^, r ~ 1, . . . ,m, such that for each q G dD there is an r G {1, . . . , m} 
for which 

q^ ^ A;^{{q[^yar{q[^)))'^), 

where g^^^ G A^, a,. €E C°'^(Ar) and Ar is as in the definition above. Since the Jacobian of each 
mapping A^. is equal to +1, the specific choice of the matrix and the vector in the definition 
of the local co-ordinate transformation A^. does not affect the argument below. We shall therefore 
assume for ease of exposition that A,. — I and = 0, so that A^. is the identity mapping. Thus, 
for example, we shall write u{q[^y q"^^-^) instead of u((A;:^((g'^^j, qj^^^))'^)'^). 

There exist functions tpr S C^(Ur), with < Lpr < \, r — \, . . . and a function ipm+i G 
C^^{D), with < (fm+i < 1, such that we have the partition of unity 

m+l 

(fij, = 1 on D. 

i=l 

Now, given u G H\,j{D), let Ur '■— uipr for r = 1,...,to + 1; clearly, u,. e H\j{D) and 
Um+i S H^{D) C H\,i{D). Having thus decomposed u G H\,i{D) as 

m+l 

(C.l) u = Ur on D, 

r=l 

we give a brief overview of the rest of the proof. Let us first note that on any compact subdomain 
B oi D the Maxwellian M is bounded above and below by positive constants; on such subdomains 
B we have that H\j{B) = H^{B). On D itself the situation is quite different: of course, H^{D) C 
HIj{D); however the converse inclusion is false as there exist elements in H\j{D), which, due 
to the decay of the weight function M to on dD, tend to infinity at a faster rate than can be 
accommodated within H^{D). The idea of the proof is therefore to translate each of the functions 
Ur, r = 1, . . . , m, in the direction of the boundary in order to shift the locations of potential 'bad' 
behaviour in Ur G Hli{D) from dD to D ; this then ensures that, for any (sufficiently small) 
positive value of the shift parameter A, the shifted function, Ur\, belongs to H^{D), and can 
be therefore well approximated in H^{D) by functions (denoted Ur\h below) that lie in C°°{D), 
thanks to the density of G°°{D) in H^{D). As (pm+i G C^{D), the function Um+i = ucpm+i has 
compact support in D, and can be therefore well approximated in H^{D) (and thereby also in 
H\j{D)) by smooth functions, without the need to form its shifted counterpart Um+i\ beforehand 
and approximate that instead; however, Um+i will be shifted in any case, by the same amount 
A as the functions Ur, r = 1, . . . ,rn, simply to ensure that a shifted counterpart of the partition 
(ICH holds. 

With this motivation in mind, we define shifted counterparts, Ur\, of the Ur by 

Ur\{q[r),q\r)) — (^(r) > 9(r) ~ ' for r 1 

Step 1. The first step in the argument amounts to showing that 
(C.2) lim Ur\ ~ Ur in HIj(D), r = l,...,m + l 



m - 



+ 



This clearly holds for r = m + 1. Let us therefore assume that r G {!,..., m}, let gr signify 
the function Ur or any of its first partial derivatives, and define Vr :— Ur ■ Clearly, Vr C Ur, 
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r = 1, . . . ,m. By the triangle inequality, with q — {q[r)'lfr))^ 

1 

l5r(g(,),(?f,))-.9r(9M,g(;)-A)pM(g)dg ' 

^ ' \9r{q[r),q^r) " A) P \[M {q[^) , q^^^ - X)]i - [M{q[^),q^^^t^\' dq 

IJ V,. 

(C.3) 



< 



Uv,. 
T1+T2. 



We begin by considering T2. Let A G (0,/3/2]; then, for /3 > sufficiently small, there exists 
Co G (0, 13] sufficiently small such that 



\[Miq[r),qh - ^)]* - [M{q[r),q'^r))]^\' < M{q[^),qf,^ - A), 

for g^^j e Ar and g"^^ e (ar('?(r)) ^ cq, ar(g(^j)). This follows on noting that for all A G {0,(3/2] 
and Co sufficiently small. 



M(<?[,),gf,)-A) >M((7[,),q[;)) for 



q{r) e {ariq[r)) ~ Co , fl^ ((?(^) )) . 



Hence, and by the absolute continuity of the Lebesgue integral, for any e > there exists c\ € (0, cq] 
small (independent of A e (0, /3/2]) such that 



(r) 



(C.4) < / dg[,,) 



l5r(9(r),'?(r) - A)|^M(g(^),g['^) - A) dg"^) < - 



'ar(9|^))-ci 

Now for ci > fixed, there exists A > sufficiently small such that 

rar(9j^))-ci 



r) 



\9Mryql.) - X)? |[M(g[.),<ZM - A)]5 - [M(g[,),qf,))]5pd< 



(r) 



< max 



[M{q[ryql)->^)]^ 



(C.5) 



X / dgf 



(r) 



ar(q[r))-ci 



\9r{q[r),qir) " A) ^ M(g[,) , qf,, - A) dgf,, < - 



e . 



Summing (|C.4I) and (jC.5|) and taking the square root of both sides of the resulting inequality, we 
deduce that for any e > there exists A > such that 



\9r{q[r),q'^r) " A)!' IIMK.),?^ - A)]t - [M(g[,),gf,))]^|2dg 



< e. 



Hence, 
(C.6) 



l5r(g[.),gr.) - X)\''\[M{q[ryq^r) ' ~ [M {q[,y q^,)P \' dq 



converges to as A — >■ 0+. That concludes the analysis of term T2. 
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Concerning Ti , continuity in the L'^-norm of the translation operator implies that 



(C.7) 



converges to as A 0+. That concludes the analysis of term Ti. 
Finally, I^M and dSIZl) imply jUJ). 

Step 2. Having shown (jC.2p . it now suffices to prove that each of the functions Ur\, r = 
l,...,r7i + l, is a limit in H\j{D) of C°°{D) functions. To that end, for r = 1, . . . , m + 1 we define 
the functions 

If ( \(1-P? \ 

(for h sufficiently small, e.g., < ft- < ^A, to ensure that the integral is correctly defined), with 
Co chosen so that 

|2 \ 



19 -P 

Co/i" J\q-p\<h I I? - Pp - 



Clearly, Ur\h € C°°{D), r = f , . . . , m, and for A e (0, /3/2] fixed, we have that 

lim UrXh — in H^{D). 

Thus, a fortiori (noting that M e L°°{D)), for A G (0,/3/2] fixed, 

YwTL UrXh^Urx in Hli{D), r = l,...,m. 

/i->0+ 

Hence, given e > and A e (0, /3/2] fixed, there exists h > such that 
(C.8) \\urXh- Urx\\Hi,iD) < — T , r = l,...,m. 

Further, using that M S L°°{D) and the closeness of Um+iA/i and in H^{D), and thereby 

in H\.j{D), we also have 

(C.9) |!Wm+lA?i - "m+1 a|1_H-i (D) < — T 

for /i > sufficiently small. 
We define 

m+l m+1 

"A/i := ^ WrAh and wa := ^ "rA- 

r— 1 r=l 

The inequalities (jC.8|) and (jC.9[) then imply that 

m+l 

(C.IO) ||uA/« - Wa|Ui^(D) < ^ ||UrA;« - ■"rA|Ui,(n) < £• 

Since the functions Urxin r = 1, . . . , m + 1, all belong to C°°{D) the same is true of uxh- Finally, 

\\u - Uxh\\Hl^(D) < \W - Ux\\h\,{D) + \\UX - Uxh\\Hl^{D). 

and therefore the stated density result immediately follows, on recalling (|C.ip . (IC.2p . the definition 
of ux, ([CTOl) and that uxh e C°°(i5). □ 



By an identical argument, C°°{^ x D) is dense in H\^{p. x D) 
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Appendix D. Compact embeddings in Maxwellian weighted spaces 

D.l. Step 1: Compact embedding of H\,j{D) into L\j{D). We are grateful to Leonardo 
Figueroa (University of Oxford) for suggesting the proof presented in Section [D.l\ 

Let D := Di X ■ ■ ■ X Dk, where Di = B{0,\/bi), hi > 0, i = 1,...,K, and suppose that 
M{q) := Mi(gi) • • • Mk^Qk)- We shall prove that 

(D.l) Hl,{D) ^ Ll,{D). 

We begin by recalling from the Appendix of Barrett & Siili [TD] that 

(D.2) iJ],^(A)^il/,(A) 

for i = 1, . . . , which was proved there using a compactness result due to Antoci We shall 
prove (|m|) for the case oiK = 2, with D = Di x D2 and M{q) = Mi(qi)M2((72)- For > 2 the 
proof is completely analogous. 

Let u S H\j[D). As M = Mi x M2, it follows from Fubini's theorem that, for almost all 

qi e Di, 

u{qi,-)eLl,{D2) and •) € ^m, (1)2), 

where a is any d-component multi-index with < |a| < 1. Fubini's theorem also implies that, 
given G C^(I?2) and a d-component multi-index ck2, with < |q!2| < 1, we have 



(-1) / u{quW^2M7 



(fi dqi = 



Di 



ipi dqi, 



for all ipi £ C^{Di). Therefore, d°"'[u{qi, ■)] = d'-°'°''-^u{qi, ■) in the sense of weak derivatives on 

D2 for almost all qi £ Di. As d''^'°'^^u{qi, •) belongs to L|j^(I?2) for almost aU qi e Di we have 
that 

(D.3) u{qi, •) e HIj^{D2) for almost ah qi e Di. 

Analogously, 

u{-,q2) £ H\,j^{Di) for almost all 92 € -D2. 

As each of the partial Maxwellians, Mi and M2, is bounded from above and below by positive 
constants on compact subsets of their respective domains, there exists a sequence (1?^ („) : n G N) 
of open proper Lipschitz subsets oi Di, i — 1,2, such that 

00 

A,(n) c A,(n+i), ?^ e N, |jA,(n)=A and 

n=l 

e.g., Di_(^n) ~ B(0, ^{^i )', the compact embeddings stated here follow by the Rellich-Kondrachov 
theorem (cf. Adams & Fournier [1], p. 168, Theorem 6.3, Part I, eq. (3)) applied on I?i^(„), i — 1,2, 
n e N. Letting, for n S N, ^(n) •= ^i=i Di,{n) Q D, and noting that, by Appendix A, i?(„) is a 
Lipschitz domain, the above properties get inherited by -D(ri) from _Dj 

00 

!)(„) C I?(„+i), n e N, y and (£>(„)) '-^ L^,^ (£>(„)). 

71=1 

Let D,'"^ := D, \ and \ £>(„)• It follows from Opic gT], Theorem 2.4, that the 

above compact embeddings on members of a nested covering imply the following characterizations 
(the first, for i e {1,2}): 

(D.4) ^iL.(A) ^ lim sup / u^dM,/\\u\\l, ,^ .=0, 

(D.5) Hl,{D) ^ Ll,{D) ^ lim sup / AM /\\u\\]j, ^ Q, 
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where dMi := Mi{qi) dqi, i = 1,2, and dAI := M{q) dq. By virtue of (ID.2p . the left-hand side of 
(ID.4[) holds; hence, its right-hand side also holds. Using (|D.3p and (jD.4[) with i — 2, we deduce 
that for any e > there exists n — n{e) G N such that 



dM 



Di 



u^{qi,-)dM2 



I3i 



dMi<e^ ||w(<Zi,-)||?^i^^^(B,)dAfi 



/ [/ u^q,,-)dM2+ f |Vq,u(gi,-)pdAf2 

JDi Ud2 ~ JD2 



dMi 



< e u 



and similarly for /^(n) 



dM. Then, as L>(") = (L>i x D^"^) U x L>2), the right-hand side 



X_D2 



of (jP.sp holds; therefore, so does its left-hand side; hence (jD.l[) . 

D.2. Step 2: Isometric isomorphisms. Let be a bounded open Lipschitz domain in W^. We 
now show the isometric isomorphism of the following pairs of spaces, respectively: L\.j{Vl x D) and 
L^{n;L\j{D)); Hlf(VLxD) and L^ifl; Hli{D)y, Hlf{nxD) smd H^{n;Llj{D)). The definitions 
of Hlf{n X D) and H];j{VL x D) are given below. 



D.2.1. Isometric isomorphism of L\.j( fl x D) and L^{Q,; L\j{D)). Let 



L\n;LlAD)):^{veMU^,Lli{D)) : f \\v{x)\\ 

Jn 



dx < 00}, 



where 



L\,j{D)) :— {v : ^ L\j{D) : v is weakly measurable on f2}. 



Let {<y5j}j^i be a complete orthonormal system in the (separable) Hilbert space L\j{D) with 
respect to the inner product (•, •) of L\j{D). For v G i^(f^; L\j{D)), we define the function 

N 

VN{x,q) ■.= ^{v{x),Lpj)Lpj{q). 

As V is weakly measurable on $7, each of the functions x 1— >■ (pj), = 1, 2, . . . , is measurable 

on r2; therefore {x,q) 1— ?> (^(a;;), V'j) is measurable on $7 x D for all j = 1,2,.... Similarly, q 1— > ^Pj{q) 
is measurable on D for each j = 1,2,..., and therefore (x, q) (pj{q) is measurable on J7 x _D. 
Hence, also Vat is a measurable function on fi x _D. Now, 

N N 

m— 1 

By the Cauchy-Schwarz inequality Mipj (pk = M^(pj ■ M^ipk G L^{D) for all j, k> \; hence also 
M(-) 11^(3;, OP € L'^{D) for a.e. x & VL. Thus, by the orthonormality of the ipj, j = 1, 2, . . . , in 

A' 

/ A/(q)|VW(x,(?)|2dq= V|(t;(3;),^,)|2, a.e. 3; G 

By Bessel's inequality in L\j{D), the right-hand side of this last equality is bounded by ||w(a;)|||2 

for a.e. x and, by hypothesis, x >->■ v{x) € ^^(0); therefore, by Fubini's theorem, M \Vm\^ & 
L^(f2 X D). Upon integrating both sides over $7, and using Fubini's theorem on the left-hand side 
to write the multiple integral over 51 and D as an integral over J7 x D, we have 



(D.6) 



[ M{q)\VNix,q)\^dqdx = y^ [ |(«fe),^,)pdx 
JnxD ~ ~ ~ "'o 
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Now, let 

N 

VnIx) :^^\{v{x),ipj)\^, xgQ. 
i=i 

The sequence {yN{x)}N^i is monotonic increasing for almost all x £ Q; also, according to Bessel's 
inequality in L\j{D) we have that 

< 2/iv(s) < Ik(s)|li2^(£,) ViV>l, a.e. jGfi. 

Thus {yN{x)}'^^i is a bounded sequence of real numbers, for a.e. x&Q,. Therefore, the sequence 
{j/jv(2)}?^=i converges in M for a.e. g € K, with 

oo 

y{x) = Ymi yN{x) \{v{x),Lpj)\^ , a.e. x G 1^. 

N—^oo ^ — ^ 

By the monotone convergence theorem. 



J™ / l(^(s)'V'i)PdS = lim / yN{x)dx 

p p OO 

(D-7) = yix)dx= Y.\iv{x),^j)\'dx. 

Jq Jn --^ 



This implies that 

N 



^ Jjiv{x),^j)\'dx' 



N=l 



is a convergent sequence of real numbers. Hence, it is also a Cauchy sequence in M. 
Since, for any N > L> 1, 

/ \VN{x,q)-VL{x,q)\^dqdx= V / \{v{x),^j)\^ dx, 
Jqxd ~ ~ ~ j^j^j^-^Jd 

it follows that {Vjv}^^^ is a Cauchy sequence in L\j{Q, x D). Since L\^{Q, x D) is a Hilbert space. 



there exists a unique V e L\i{^ x D) such that 



(D.8) V= lim Vn mL\i{VlxD). 

Thus we have shown that the mapping 

oo 

I:v& L\n,LUD)) ^V:= J2{v{-),^j) € 1.^(11 x D) 

is correctly defined. Next, we prove that X is a bijective isometry, and this will imply that the 
spaces L'^{0,; L\j{D)) and Ll^{fl x D) are isometrically isomorphic. 

We begin by showing that X is injective. As I is linear it suffices to prove that if I{v) = then 
V = 0. Indeed, if I{v) = 0, then 

oo 

''^{v{x), ipj) <Pj{q) = for a.e. {x, q) G x D. 

Since {<fij}jLi is an orthonormal system in L\f{D), it follows that {v{x),(fij) = for a.e. x € 

and all j = 1,2, The completeness of the orthonormal system {(pj}j±-^ in L\j{D) now implies 

that v{x) = in Ll_j{D) for a.e. x e n, i.e. w = in L'^{n; L\j{D)). 

Next we show that X is surjective. Suppose that V S L\i(Q. x D). Then, by Fubini's theorem, 
V{x,-) e L\j{D) for a.e. x & Q,. Since {<^j}^j is a complete orthonormal system in L\j{D), it 
follows that 

•oo 

V{x,-) = Y^{V{x,-),^i)^i{-). 
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On defining v{x) := V{x, ■) E L\j{D), we have that I{v) — V . Hence X is surjective. 
Finally, we show that T is an isometry. Clearly 



ED 
EH 



N 



oo 



Applying Parseval's identity in L\j{D) to the infinite series under the last integral sign, we deduce 
that 

Thus we have shown that ||Iw||L|,j(oxr>) — II^IIl2(0;L2(£))), whereby T is an isometry. 

D.2.2. Isometric isomorphism of H*^^ (ft X D) and L"^ H\^(D)) . Let us begin by observing that 
L\j{Vl X D) C Ll^^{il. X D), and therefore any V in L\,f{rt x D) can be considered to be an element 
of 2?'(il X D), the space of R- valued distributions on Vl x D. Let Vg denote the distributional 
gradient with respect to g, defined on 2?'(il x D). We define 

H^t^jin X D) := {V e Llj{n X D) : V qV G L\j(9. x D)}. 

A completely identical argument to the one above shows that x D) is isometrically iso- 

morphic to L'^{V,; HIj{D)); the only change that is required is to replace L\j{D) by H\,i{D) 
throughout and to take {'/'jj^i to be a complete orthonormal system in the inner product (•, •) 
of the (separable) Hilbert space H\j{D), instead of L\j{D). 

D.2.3. Isometric isomorphism of {Q. X D) and L\j{D)) . Let 

Hlf{n X D) := {V e Llj{n X D) : V .^V G L|^(17 x D)}. 

Concerning the isometric isomorphism of H]^j{Vt x D) and H^{Q,; L\j{D)) we proceed as follows. 
Given v £ H^i^; Lfj{D)) C L^{^', L\i{D)), we define, as in the proof of the isometric isomorphism 
of L'^{VL;L\j{D)) and L\i{p. x D) above, the function 

oo 

V : {x,q)enx D>^V{x,q) := ^(^(a:), (g) G R, 

i=i 

where {Vjlj^i is a complete orthonormal system in L\j{D). We showed above that V E x 
D), and ||l^||L|,,(oxr>) = \ML^(n;Ll,{D))- 

Now, let Yx denote the distributional gradient with respect to x, defined on T)'{il x D), and 
let Px denote the distributional gradient, defined on V {V,; L\^{D)), the space of L|^(-D)-valued 
distributions on fl. Applying V^; to 

CO oo 

V = ^(w, Vj) Vj in VlVl X D) and noting that V^V^ = ^(-PxW, Vj) Vj^ 
i=i i=i 

it follows from the isometric isomorphism of L\i{Vl x D) and L"^ {Vt\ L\^{D)) that 

= M\ll,{n;Ll,[D)) + \\DM\l2(^n;Ll,[D)) = M\]mn;Ll,{D)y 
which shows that Il\f{Q. x D) and L\j{D)) are isometrically isomorphic. 
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D.3. Step 3: Compact embedding of Hl^iifl x D) into ^^^(^1 x D). We use the results 
of Step 2 to identify ttie space L\j{n x D) witli L"^ [Vl] L\[{D)) and the space H\^{VL x D) — 
Hlf{n X D)n H^j^n x D) with H^{n;Llj{D)) n L'^{n;Hlj{D)). Upon doing so, the compact 
embedding of Hlj{il x D) into Lj^^^ x D) directly follows from the compact embedding of 
H^{n;Llj{D))nL'^{n;Hli{D)) into L2(f]; L|^(D)), implied by Theorem 2 on p. 1499 in the paper 
of Shakhmurov [46] . thanks to the compact embedding from Step 1. 
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